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ABSTRACT 


PHASE  1 :  FOAMS 

Results  of  optical  measurements  of  peak  pressures  in  the  0.1-  to 
6-kbar  range  transmitted  through  3-  to  12-mm-thick  layers  of  certain 
foams,  namely,  0.7  to  1.4  g/cc  aluminum,  0.67  g/cc  polyurethane, 

1.1  g/cc  beryllium,  1.1  and  1.7  g/cc  graphite,  and  1.0  g/cc  silica,  are 
presented.  Superficial  impact  momentum  densities  are  in  the  range  of 
0.6  to  3  x  104  taps  and  are  delivered  by  explosive  propulsion  of  a  thin 
sheet  of  aluminum. 

Within  the  range  of  study,  peak  pressures  scale  approximately 
linearly  with  the  ratio  of  superficial  momentum  density  to  superficial 
mass  density  of  foam;  the  range  of  validity  of  the  scaling  and  the  scale 
factor  depend  strongly  on  the  kind  of  foam.  The  effect  of  degree  of  dis¬ 
tention  on  scaling  relations  in  aluminum  is  slight  but  is  somewhat  greater 
in  graphite.  Of  the  foams  considered,  polyurethane  transmits  the  least 
peak  pressure  per  unit  of  added  weight  but  shows  the  greatest  tendency  to 
increase  the  delivered  momentum  through  rebound. 

Through  the  use  of  parameter  values  measured  in  shock  experiments  an 
elementary  theory  predicts  fairly  reliably  the  conditions  under  which  the 
transmitted  shock  structure  consists  of  an  elastic  forerunner  alone  with¬ 
out  a  following  locking  wave.  Presently  available  theory  fails  to  predict 
accurately  either  observed  peak  pressures  or  wave  shapes. 

PHASE  2:  SOLIDS 

Attenuation  of  shock  waves  is  studied  in  specimens  made  of  1060  and 
2024  aluminum,  OFHC  copper,  gold,  and  Armstrong  C-7  epoxy.  The  strains 
induced  in  the  specimens  are  virtually  one- dimens iona 1 .  Comparison  of 
experimental  and  calculated  results  indicate  than  an  e  1  as  top  1  as t i c  model 
for  the  relation  between  stress  and  strain  should  be  used  for  aluminum. 


•  *  • 
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copper,  and  the  epoxy.  Such  a  relation  permits  the  prediction  of  the 
observed  attenuation  more  exactly  than  does  a  relation  in  which  rigidity 
is  neglected. 

Results  for  gold  show  that  the  shock  is  also  attenuated  more  rapidly 
than  predicted  by  the  use  of  a  s t r ess- s t ra i n  relation  which  neglects 
rigidity.  More  data  are  needed  before  a  more  appropriate  s t ress- s t ra i n 
relation  can  be  obtained  for  gold. 
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1.  I  INTRODUCTION 


An  interposed  layer  of  porous  or  foamed  material  can  often  sharply 
reduce  the  peak  stress  reached  in  a  hard  body  as  a  result  of  deposition 
of  X-ray  energy.  The  foam  may  then  be  a  means  of  protecting  the  body 
against  pressure  damage.  We  have  studied  this  property  in  several  widely 
different  foams  useful  in  lowering  impact  stresses  due  to  collisions  last¬ 
ing  10  to  100  /xsec  in  which  momentum  densities  in  the  neighborhood  of 
1  x  104  dyne  sec  cm  2  (taps)  are  exchanged.  The  first  goal  has  been  ex¬ 
ploration  of  the  effect;  the  second,  semi qu an ti tati ve  prediction  using 
easily  accessible  physical  properties  of  foams;  and  the  third,  under¬ 
standing  of  the  detailed  flow  within  the  foam  during  the  impact.  Methods 
of  investigation  have  included  measurements  of  static  properties  of  foams, 
observations  of  foam  behavior  under  shock,  and  calculations. 

The  work  on  this  and  on  previous  contracts1,2  has  yielded  consider¬ 
able  exploratory  information  ondistended  plastics,  silica,  graphite, 
beryllium,  and  several  kinds  of  aluminum.  In  all  these  materials  it  is 
established  that  there  is  a  protective  effect  and  that  there  are  two 
mechanisms  responsible:  the  two-wave  nature  of  the  flow  in  certain  dis¬ 

tention  ranges  of  these  foams  and  the  relatively  large  disparity  between 
crushing  wave  speed  and  speed  of  overtaking  rarefactions.  A  heavy  enough 
impact  upon  a  foam  layer  gives  rise  at  the  struck  surface  to  a  wave  of 
pronounced  change  moving  into  the  foam;  the  voids  are  removed,  and  the 
density  greatly  increased.  If  high  stress  is  not  kept  on  the  struck 
surface,  release  waves  start  from  that  surface  toward  the  shock  front  and, 
because  of  the  large  changes  in  the  nature  of  the  material  brought  about 
by  the  crushing  or  main  wave,  these  release  waves  quickly  overtake  the 
impact  wave  and  steadily  reduce  its  speed  and  pressure.  This  tends  to 
lower  the  pressure  finally  delivered  to  the  protected  structure.  Some¬ 
times  splitting  off  from  and  moving  ahead  of  the  main  wave,  another  wave 
of  much  less  density  change  may  carry  off  at  high  speed  an  important 
part  of  the  total  impact  momentum  at  the  first  or  outer  foam  interface 
and  deliver  it  to  the  structure  to  be  protected  at  a  relatively  low,  con¬ 
stant  pressure  applied  for  a  relatively  long  time,  this  wave  is  related 
to  the  linear  or  proportional  region  in  a  statically  measured  stress-strain 
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diagram.  Since,  when  it  is  present,  it  comes  to  the  structure  ahead  of 
the  main  or  crushing  wave  we  call  it  the  first  or  forerunning  wave. 

As  is  plain  from  the  summary  in  Section  2  of  this  report,  much  re¬ 
mains  to  be  learned  of  the  generally  true  details  of  flow  in  foams.  The 
further  research  outlined  in  the  fifth  and  last  section  would  help  to 
throw  light  on  these  details.  Sections  3  and  4  contain  descriptions  of 
this  year's  experiments  and  calculations,  respectively. 
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2.  SUMMARY 


A.  EXPERIMENTAL  METHODS  AND  RESULTS 

He  explosively  put  a  certain  amount  of  momentum  per  unit  area  in  the 
range  0.5  to  5  x  104  taps  into  one  side  of  a  foam  layer  3  to  12  mm  thick 
and  by  high-speed  photography  of  light  images  formed  in  a  polished  sur¬ 
face  measured  the  shapes,  pressures,  and  speeds  of  stress  waves  in  a  steel 
anvil  in  contact  with  the  opposite  side  of  the  foam.  The  basic  method  is 
the  optical  lever  arm  described  by  Fowles.3  Except  for  a  few  cases  in  which 
the  explosive  lay  directly  against  the  foam,  momentum  was  delivered  by  an 
aluminum  sheet.  0.012  to  0.040  inch  thick.  For  one -d i mens i ona 1  (1-D) 
symmetry  the  sheet  was  driven  by  an  explosive  lens  across  a  vacuum  gap 
onto  the  foam;  in  this  arrangement  the  anvil  is  wedge-shaped.2  The  more 
usual  arrangement  allowed  a  small  amount  of  obliquity  and  involved  a  run¬ 
ning  detonation  front  in  sheet  explosive  lying  against  the  aluminum  sheet; 
these  are  called  two-dimensional  ( 2 - D)  experiments.2  For  the  1-D  experi-* 
ments  we  know  accurately  the  value  of  momentum  density  applied  to  the 
foam;  in  the  2-D  the  relation  between  the  momentum  densi tv  and  the  explo¬ 
sive  thickness  is  known  only  to  an  accuracy  of  ±10  percent.4  Because  of 
its  ease  of  performance  and  because  of  a  wish  to  survey  behavior  in  a 
wide  variety  of  foams,  we  used  the  2-D  arrangement  this  year  in  all  exper¬ 
iments.  Figure  l  is  an  illustration  of  the  arrangement  including  all  its 
variations;  Fig.  2  is  a  record  made  with  the  basic  arrangement,  i.e., 
flyer  in  contact  with  a  plane  parallel  slab  of  foam  lying  on  a  plane  par¬ 
allel  steel  anvil.  The  light  source  responsible  for  the  images  seen  in 
the  figure  is  an  explosive  argor.  bomb  with  a  ruled  translucent  grid  de¬ 
scribed  bv  Pressman.'  The  small  amount  of  obliquity  must  always  be  kept 
in  mind  when  comparing  peak  pressures  and  wave  shapes  in  1-D  and  2-D  ex¬ 
periments  in  the  same  thicknesses  of  the  same  material,  however,  in  the 
survev  comparing  results  in  widely  different  materials,  obliquity  is  not 
critically  important. 

The  kinds  of  foam  studied  this  year  have  been  polyurethane,  open- 
and  closed-cell  aluminum,  silica,  beryllium,  and  graphite. 
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FIG.  1  CROSS-SECTIONAL  VIEW  OF  GENERAL  EXPERIMENTAL  ARRANGEMENT 
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FIG.  2  TYPICAL  SMEAR  CAMERA  RECORD 


For  each  shot  Table  1  lists  foam  kind,  densi'y,  and  thickness,  then 
explosive  thickness  (approximately  0,016  inch  provides  1  *  104  taps,  except 
when  there  is  an  air  gap  between  flyer  and  foam  0.0243  inch  ’  1  *  104  taps), 
and  anvil  thickness.  Sometimes  we  measured  detonation  speed  with  time- 
of-arrival  switches  or  “pins"  on  the  explosive,  always  the  smear  camera 
photographed  the  motion  of  images  in  the  anvil  free  surface.  When  the 
aluminum  sheet  (called"  flyer"  or  ‘‘barrier’’)  was  separated  by  an  air  gap 
from  the  foam  or  when  the  foam  was  preheated,  “standoff”  or  " ho  t ,  "  r espec - 
tively,  was  entered  in  the  column  headed  “Type  of  Experiment.”  Two  anvils 
with  nonparallel  surfaces  are  also  noted  in  this  column. 

In  Table  1  the  emergent  waves  are  identified  by  the  letter  A  for 
elastic  forerunner  and  B  for  main  or  locking  wave.  When  the  first  wave 
is  double,  its  parts  may  be  signaled  by  A- 1  or  A-2.  Pressure  jumps  (not 
peak  pressures)  in  each  wave  are  entered  after  each  wave  designation  and 
are  presumed  to  exist  at  the  foam-anvil  interface  although  they  are  mea¬ 
sured  at  the  anvil  free  surface.  Apparent  wave  speed  means  the  speed  of 
the  front  along  this  free  surface  and  within  experimental  error  should 
equal  detonation  speed  (7.35  t  0.07  mm/^sec  at  ordinary  temperatures) 
uniess  a  beveled  anvil,  tapered  slab,  or  tilted  fiver  is  used.  In  each 


experiment  involving  a  preheated  specimen  detonation  speed  was  measured 
again,  and  these  results  are  noted  in  Table  1.  Average  wave  speed  (in 
the  foam)  is  inferred  from  the  separation  [D  in  Fig.  1(b)]  along  the 
anvil  mirror  surface  between  detonation  and  point  of  wave  emergence  by 
assuming  a  planar  wave  front  in  the  foam  as  well  as  in  the  anvil. 

Because  this  deduced  value  depends  strong'y  on  the  value  taken  for 
the  apparent  wave  speed,  7.35  mm'/xsec  is  used  whenever  the  experimental 
measurement  of  apparent  speed  agrees  with  that  value  within  the  uncer¬ 
tainty  of  the  measurements.  In  computing  the  pressure  jumps  from  the 
camera  record  the  same  rule  is  followed,  although  the  computation  is 
much  less  sensitive  to  apparent  speed.  It  should  be  noted  that  the  aver¬ 
age  speed  of  a  forerunner  as  read  in  Table  1  will  be  higher  than  its 
characteristic  speed  by  an  amount  increasing  with  the  delay  between  ini¬ 
tial  impact  and  forerunner  emergence  ahead  of  the  main  shock. 

Figure  3  shows  approximate  pressure  histories  at  the  foam-anvil 
interface  sketched  from  the  smear  camera  photographs.  Written  beside 
each  outline  is  the  total  impulse  showing  in  it  up  to  the  point  on  the 
time-axis  marked  by  the  arrow,  when  the  observation  was  interrupted  bv 
reverberations  or  the  pressure  returned  finally  to  zero. 

B.  STATE  OF  OUR  UNDERSTANDING  OF  FOAM  BEHAVIOR 
UNDER  SHOCK  LOADING 

1.  Tin  Fohf.ri'nnkh 

In  some  solids  a  single  impact  may  give  rise  to  a  train  of  two  shock 
fronts  of  which  the  speed  and  pressure  of  the  first  are  characteristic  of 
the  material. h  El  a s 1 i r - r i gi d  foams  often  support  two  successive  waves  stem 
ming  from  the  same  impact.  Noticing  that  strain  hi  quasi-stat  u  ,  l-D  com¬ 
pression  of  elastic*  rigid  foam  is  proportional  to  stress  up  to  a  certain 
limit  (defined  as  the  stress  where  strain  falls  out  of  proportionality  h\ 
0.2  percent)  or  yield  stress,  we  have  sought  to  identify  this  elastic  or 
proportional  limit  state  with  the  state  behind  the  first  or  forerunning 
wove;  hut  in  the  one  e  1  ast  i  c  -  r  l  gi  d  foam  examined  accurately-'  fm  the  cor¬ 
relation  (0.67  g  cc  polyurethane)  we  have  forecast  a  forerunner  speed 
from  the  mouulus  of  linear  compression  whirh  is  20  percent  lower  than 
measured.  This  discrepancy  may  be  due  to  1  ark  of  uniformity  in  the  phys¬ 
ical  properties  among  specimens  of  the  same  kind  and  similar  density.  In 
none  of  the  several  foams  studied  has  the  proportional  limit  or  yield 
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stress  amounted  to  more  than  half  the  stress  seen  in  a  steel  anvil  struck 
by  the  forerunner.  This  can  only  be  due  to  the  effect  of  the  rate  of 
strain  although  the  magnitude  of  the  effect  may  be  less  than  appears  in 
the  disparity  between  the  pressure  induced  by  the  forerunner  in  the  anvil 
and  the  yield  stress,  since  the  free-running  forerunner  stress  is  enhanced 
upon  reflection  from  the  wall.  Because  the  anvil  is  relatively  motionless 
under  the  impact,  the  product  of  pressure  and  volume  increments  across 
the  reflected  shock  must  equal  the  same  product  across  the  forerunner; 
hence  for  any  likely  equation  of  state  the  first  wave  stress  cannot  more 
than  double  itself  at  the  wall. 

Although  we  have  observed  forerunner  speed  in  only  one  foam  mater¬ 
ial,  e.g.,  polyurethane,  we  have  measured  average  speeds  (Table  1  of  first 
waves  in  several  foams  by  an  addition  to  the  experiments  by  which  we  ob¬ 
served  press" re  histories.  Since  the  forerunner  may  not  break  out  ahead 
of  the  main  wave  until  the  main  wave  has  slowed  to  the  characteristic 
forerunner  speed,  observed  average  first  wave  speeds  are  upper  limits 
upon  this  characteristic  speed.  In  all  experiments  except  one  involving 
0.9  g/cc  open-cell  aluminum  and  another  in  polyurethane  (C-P-R),  observed 
average  wave  speeds  agree  within  experimental  uncertainty  with  modulus 
speed  or  are  higher  than  modulus  speed.  If  strain-rate  affects  forerunner 
speed  the  influence  must  be  less  than  20  percent. 

As  the  forerunner  front  moves  through  polyurethane  it  tends  to  lose 
its  abruptness  or  discontinuous  character.'  Si 'tea  foam  has  never  shown 
o  steep  pressure  rise  in  the  forerunner  but  r<nularly  transmits  first 
waves  reaching  end  stress  in  two  distinct,  gradual  stages  usually  more 
than  a  microsecond  apart.8  Such  behavior  suggests  a  compound  nature  of 
the  first  wave  a  purely  elastic  wave  mav  be  only  a  part  of  the  forerunner. 
Or  the  loss  of  sharpness  mav  be  u  feature  of  stress  relaxation  in  the 
fo  re  runt- e  r of  which  we  also  have  evidence  in  both  polvurethnne  and  alu¬ 
minum. w  Closed-cell  aluminum  on  the  other  hand  has  always  given  a  sharp 
forerunner,  open-cell  aluminum  and  graphite  have  shown  both. 

Evidence  for  the  compound  nature  of  the  forerunner  ;s  our  observa¬ 
tion  that  measured  sound  speed  has  invariably  been  significantly  higher 
than  the  speed  predicted  fiom  the  modulus  of  linear  compression.  The 
ratio  of  sound  speed  to  modulus  speed  ranges  from  3  in  0.9  g  rr  open-cell 


**  [>•*••(•  of  *  b«r*l  of  1  Ik  »<jp*nno)<  «i*«r(i»  ihroufK  a  iluclnaaa  ■>(  (<>»*. 
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aluminum  and  1.2  g/cc  beryllium  to  about  1  1/3  in  0  67  g  rr  polyurethane. 
One  measurement  of  modulus  in  si  lira  yields  u  speed  agreeing  with  sound 
speed. 

Most  of  the  average  wave  speeds  are,  however,  lower  than  sound  speed 
in  the  same  kind  and  density  of  material,  which  is  further  evidence  for 
the  complex  nature  of  the  forerunner. 

The  foregoing  comments  on  the  forerunner  are  summarized  in  Tables  2 
and  3.  In  Table  2  for  each  kind  and  density  of  foam  studied  appears  a 
sound  speed  (column  2),  one  to  thre-  observed  shock  speeds  (columns  3,  4, 
and  5)  and  the  forerunner  speed  forecast  from  the  modulus.  Of  the  ob¬ 
served  shock  speeds  that  under  Method  A  (column  3)  is  from  a  flash  X-ray 
shadowgram  reported  by  Fowles  and  Curran,1  that  under  Method  B  (column  4) 
stems  from  simultaneous  optical  measurements  of  stress  and  speed  in  a 
foam  wedge  reported  by  Rempel,*  and  those  under  Method  C  (column  4)  are 
the  average  wave  speeds  in  Table  i. 

Table  3  lists  for  each  kind  and  density  of  foam  a  value  of  static 
yield  stress  (column  3).  The  tangent  to  the  quasi-static  st ress- st rain 
curve  at  the  point  where  the  second  derivative  of  stress  with  respert  to 
strain  vanishes  meets  the  extension  of  the  linear  portion  at  the  stress 
value  shown  in  column  4.  Column  5  reports  the  observed  stresses  in  a 
steel  anvil  resulting  t rom  forerunner  impact;  and  in  column  6  are  the 
values  of  stress  induced  in  a  rigid  wall  if  the  foam  Hugoniot  followed 
the  linear  portion  to  the  intersection  in  column  4  and  then  fell  onto 
the  tangent. 

Except  for  the  peak  pressure  behind  it  we  know  nothing  of  the  nature 
of  the  wave  reflected  bv  the  forerunner  off  the  anvil.  It  would  be  a 
broad  compression  fan  if  the  shock  equation  of  state  were  like  the  quasi- 
static  in  shape  and,  if  it  were,  we  could  understand  why  we  hove  not  seen 
the  shadow  of  the  w«  1  1  -  re  f  1  ec  t  ed  wave  in  X-ray  shadowgrams.11 

2.  The  Main  Wave 

Two  aspects  of  the  main  wave  are  important  for  our  protective  pur¬ 
poses:  what  is  the  peak  pressure  and  momentum  it  brings  to  the  foam- 

covered  structure  and  how  much  foam  is  needed  to  stop  it  completely7 
Our  experimental  answers  to  these  questions  in  specific  cases  have  ap¬ 
peared  in  Table  1  and  Fig.  3.  We  have  generalized  these  results 
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Table  2 


COMPARISON  of  measurements  of  forerunner  speed 
IN  FOAMS  WITH  PREDICTIONS 


MATERIAL 

SHOCK  WAVE  SPEED  ( mm/fittc ) 

AND 

Meaaured 

Forecaat 

froa 

Modu 1 u  a 

(3) 

Utn Oil! 

(g/cc) 

Method  A* 

Method  B* 

Method  C§ 

( average ) 

P-C  -0.6 
(Poly-  n  L(. 

iron)  °'66 
0.644 

1.73  (1)A 

1.74 

1.65  t  0.05 

1.8  t  0.3 

1.40  i  0.  10 

0.645 

1.65  t  0. 10 

P-C  0.714 

(CPR)  Q  79g 

2.11  (3) 

1.2  ±  0.10 

1.50  t  0.08 

Si  -1.1 

1.35  (1) 

1.02 

1 

1.0  i  0.05 

0.931 

1.3  ±  0.065 

1.048 

!  1.4  ±  0.3 

1 . 053 

1.8  ±  0.4 

MD-AK  1.41 

4.6  (3) 

2.06  ±  0. 10 

1.445 

2.  14  ±  0.27 

1.458 

2.45  i  0.30 

MD-AO  0.948 

2.54  (3) 

0.83  i  0.04 

0.925 

0.57  i  0.01 

1.30 

2.08  i  0. 10 

1.279 

1.99  i  9.45 

1.41 

2.69  t  0.22 

1.38 

4.75  (3) 

C-ATJ  1.78 

2.32  (3) 

1.32  i  0.  13 

-1.73 

1.0  i  0.  15 

C-PT-  1.18 
0114  -1.15 

1.8  (3) 

l.o  t  0. 15 

1.02  t  0.09 

Hr  1.195 

6.71  (3) 

2.175  *  0.36 

*  Foalaa  and  Curran. * 

*  Raapal.2 

^  Froa  work  dona  t h 4 •  yaar.  The  average  aave  tpaada  ara  froa  Tabla  I, 

•nd  ara  thoae  of  the  leading  eleaenta  in  tha  disturbance. 

Nuabere  in  perentheeee  indicate  the  firat  report  of  tba  aaaauraaanl. 

All  entriea  ara  baaad  on  conaidarat ion  of  a  aingla  aaaauraaanl  or  aipariaant 
eicept  tha  firat  antry  in  tha  laat  coluan.  Thia  foracaat  ta  enaputed  froa 
tha  average  ratio  cf  danaity  to  aodulua  apaad  aaong  S  aaaauraaanta  (Tabla  h). 
Tha  uncartamty  liaita  corraapond  to  tha  atandard  donation  of  tha  5  ratioa 
froa  tha  aaan. 
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Table  3 


FORERUNNER  STRENGTH  OF  FOAMS  COMPARED 
WITH  QUASI -STATIC  Y I El.D  STRENGTH 


MATERIAL 

DENSITY 
( g/cc ) 

STATIC 

YIELD  STRESS 
(0.2*  0 f  f » e t ) 
(±0.01  kbar ) 

INTERSECTION 

OF  TANGENTS 

TO  STATIC  CURVE 
(kbar) 

— 

PRESSURE 
INDUCED  IN 
STEEL  BY 
FORERUNNER 
(kbar) 

MAXIMUM 

PRESSURE 

IN  RIGID 

WALL  EXPECTED 
FROM  STATIC 
COMPRESSION 
(kbar) 

P-C 

(Poly  Iron  ) 

0.660 

^0.67 

0.644 

0.668 

0.33 

0.38 

0.91  ±  0.  15 

1.0  ±  0.20 

0.48 

P-C.(CPR) 

^0.714 

0.798 

0.697 

0.55 

0.70 

1.9  ±0.32 

0.96 

Si 

11 

0.  165 

0.097 

0.152 

0.66  ±  0. 15 

0.  19 

MD-AK 

1.41 

1.46 

1.445 

1.405 

0.76 

0.945 

2.6  ±  0.3 

1.49  i  0.05 

1.  11 

MD- AD 

m 

0.069 

0.33  ±  0.03 
0.37  ±  0.05 

0.  10 

C-ATJ 

wm 

0.50 

0.54 

0.665  ±  0.04 

0.80 

C-PT-01 14 

1.18 

M.15 

0.055 

0.06 

0.72  ±  0.07 

0.10 

MD-AD 

0.345 

0.61 

0.60  i  0.05 
0.90  ±  0.04 
1.49  ±  0.f5 

0.73 

Be 

1.  195 

0.28 

0.46 

1.5  ±0.  15 

0.62 

somewhat  in  Fig.  4,  where  all  our  experiments  in  aluminum  foam  we  plot 
the  peak  observed  pressure  in  the  anvil  against  the  quantity 


where  / Q  is  the  thickness  of  sheet  explosive  converted  to  momentum  den¬ 
sity  according  to  the  relation  0.016  inch  of  explosive  ( EL -  50  6  D )  equals 
1  *  10*  taps,  l  is  the  original  foam  thickness,  and  r  Q  original  foam  den 
sity.  In  Fig.  4  we  have  tried  to  allow  for  errors  in  effective  foam 
thickness  arising  from  the  obliquity  of  impact  and  for  the  uncertaints 


in  the  conversion  factor  between  explosive  thickness  and  input  momentum 
density  by  the  addition  of  horizontal  uncertainty  bars  through  the  data 
points. 

All  data  points  in  Fig.  4  stemming  from  open -cell  aluminum  lie  neer 
a  straight  line;  that  is,  within  the  ranges  of  parameter  variations 
studied  ,  peak  pressure  passing  into  a  wall  depends  linearly  upon  / n/ >  Ql . 
It  is  probably  oversimplifying  to  draw  a  single  straight  line  through  data 
points  for  such  a  wide  variety  of  densities  of  foams;  but  certainly  dis¬ 
tention  does  not  sharply  affect  the  relation  between  peak  pressure  and 
Iq/PqI'  Figure  4  may  show  a  small  difference  due  to  the  closed-  or  open¬ 
cell  structure.  Furthermore,  the  value  of  I^/p^l  at  which  the  main  wave 
stops  just  short  of  the  wall  is  near  the  intersection  of  the  straight 
line  mentioned  above  and  the  pressure  level  induced  in  the  anvil  by  the 
forerunner.  In  the  figure  this  intersection  lies  between  W  ■  ‘-95 
and  2.7  *  104  cm/sec  for  1.4  g/cc  closed-cell  aluminum  and  below 
1.6  x  104  cm/sec  in  0.9  g/cc  open  cell  foam. 

In  Fig.  5  we  have  similarly  organized  experimental  results  in  poly¬ 
urethane,  beryllium,  silica,  and  two  kinds  of  graphite.  All  experiments 
providing  information  for  Fig.  4  were  2-D,  in  Fig.  5,  despite  the  fact 
that  values  of  momentum  and  foam  thickness  are  not  strictly  comparable 
between  1-D  and  2-D  experiments,  we  have  included  data  from  1-D  experi¬ 
ments  desc r i bed  earl i er 2  (Shots  9155,  9180,  9216,  9228,  and  9217  in  poly¬ 
urethane,  and  Shots  9325  and  9345  in  silica)  as  well  as  data  from  2-D  shots. 
Values  of  input  momentum  for  2-D  experiments  were  calculated  bv  means  of 
the  theory  described  in  Section  3E  below.  For  1-D  experiments  input 
momentum  was  computed  from  measured  flyer  speed  and  mass. 

Polvurethane  (0.67  g/cc)  data  from  the  three  1-D  experiments  can 
easily  be  represented  by  a  str.iight  line  which  crosses  the  ordinate  for 
forerunner  strength  in  the  anvil  at  /Q  *  1.98  x  I04  cm. see.  Except 
for  Shots  8775  and  8863*  the  2-D  shots  m  this  material  give  results  con¬ 
sistent  with  this  straight  line.  We  do  not  know  the  reason  for  the  in¬ 
consistent  results  from  Shots  8775  and  8863;  they  may  indicate  breakdown 
of  the  scaling  law  in  the  region  of  high  thicknesses  and  momenta. 


Since  the  eecursion  of  the  enrtl  tirror  due  to  the  nein  were  in  Shot  •  8**1  end  enhenced  b* 

coincidence  eilh  the  f t r • t  reeerberst ion  in  the  en»il  »f  the  forerunner,  the  *slue»  for  wein  •••» 
peek  strength  reported  eerlier-  he*e  been  reduced  by  forerunner  strength  (1  kbsr)  in  Fig. 
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FIG.  4  EFFECTIVENESS  OF  ALUMINUM  FOAM 
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Judging  from  the  one  experiment  in  0.335  g/cc  po 1 yu reth ane  (Shot  9228 ) , 
in  which  there  was  no  main  wave,  the  line  summarizing  peak  pressures  may 
lie  signi ficantly  to  the  right  of  that  for  0.67  g/cc  material. 

* 

A  strong  anomaly  appears  in  the  silica  data  represented  in  Fig.  5. 
Shots  9325,  9893,  and  9891  are  consistent  in  the  sense  established  above 
for  polyurethane  and  aluminum;  there  is  a  linear  relation  (Silica  1  in 
Fig.  5)  between  peak  anvil  pressure  and  the  quantity  /0/P0l  that  indicates 
the  main  wave  is  just  held  6ff  the  wall  near  I^/p^l  *  1.24  *  104  cm/sec. 

But  data  from  Shots  9931,  9893,  and  9345  seem  to  fall  along  another  whol  ly 
different  line  (Silica  2  in  Fig.  5)  which  does  not  cross  the  ordinate  for 
forerunner  strength  at  any  point  in  the  first  quadrant.  A  degree  of  con¬ 
sistency  would  arise  from  discarding  the  results  of  Shot  9931  but  we  know 
of  no  reason  that  experiment  should  not  be  considered.  We  propose  as  a 
tentative  explanation  that  the  value  of  the  quantity  I fl/ p Q1  corresponding 
to  emergence  of  the  first  wave  from  main  wave  just  at  the  wall  may  sepa¬ 
rate  regions  of  different  rates  of  change  of  peak  pressure  with  J0/p0Z. 

Shot  9325  had  no  forerunner;  Shot  9931  did.  Since  the  forerunner  draws 
off  momentum  from  the  main  wave,  it  is  reasonable  that  peak  pressure 
falls  faster  with  increasing  foam  thickness  in  the  neighborhood  of  the 
first  appearance  of  the  forerunner  than  elsewhere. 

Results  in  both  kinds  of  graphite  are  consistent  in  the  sense  de¬ 
scribed  above.  ATJ  graphite  (1.73  g/cc)  stops  the  main  wave  at 
I0/P0l  ="  0.7  x  10 4  cm/sec;  PT-0114  graphite  (1.1  g/cc)  at 
I0/p0l  =  0.4  x  104  cm/sec.  Neither  of  these  crossings  is  confirmed  by 
experimental  data  from  shots  in  which  the  main  wave  is  completely  held 
off  the  wall. 

A  straight  line  through  the  two  data  points  in  beryllium  crosses 
the  ordinate  corresponding  to  forerunner  strength  at  /0//0  l  **  0. 2  x  104  cm/sec. 
As  in  graphite  this  generalization  of  foam  "  stopping  power’*  i  s  tentati  ve. 
Since  the  forerunner  broke  away  from  the  main  wave  very  close  to  the 
wall  in  Shot  10,446,  the  line  summarizing  the  beryllium  information  in 
Fig.  5  may  be  less  steep  than  it  would  be  i f  there  were  data  for  higher 
values  of  the  abscissa. 

Experimentally  observed  momentum  information,  extracted  from  Table  1 
and  Fig.  3  and  Rempel,2is  collected  in  Table  4.  Reports  by  Abrahamson 
that  explosive  impulse  delivered  to  a  slab  lying  against  sheet  explosive 
depends  on  the  density  of  the  slab  have  led  us  to  convert  thickness  of 
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'fable  l 

EXPERIMENTALLY  OBSERVED  IMPULSE  GAINS 


MATERIAL,  NOMINAL 
DENSITY,  SHOT  NO. 

INPUT 
MOMENTUM 
DENSITY,  I0 

(  10*  t spa  ) 

IMPULSE  IN  SMEAR 
CAMERA  RECORD,  IR 
(10*  taps) 

IMPULSE 

GAIN 

!R 

io 

Polyurethane, 
(Poiytron ) 

( 0 . 6 i  g/cc) 

incomplete^ 

9832' 

3.1 

2. A. 

9833 

3.1 

2.35, 

incomplete 

ifipH 

8863 

3.1 

2.5. 

incomplete 

9135  (1-D)$ 

1.1 

1.6, 

incomplete 

9 

9180  (1-D) 

1.1 

0.9, 

incomplete 

9216  (1-D) 

2.2 

1.1, 

incomplete 

9217  (1-D) 

i.i 

1.2, 

incomplete 

-- 

Polyurethane 

(C-P-R) 

(0.71  g/cc) 

10261  (standoff 

0.95 

1.15 

1.3 

10331  (standoff) 

( 109°C) 

0.<»4 

1.3 

1.6 

10332  (standoff) 

0.94 

1.3 

1.  3 

Open-cel 1 

Average 

1.3 

Aluminum 

(0.7  g/cc) 

9809 

1.0 

1.2 

1.2 

Open-cell 

Aluminum 

(0.9  g/cc ) 

9800 

1.6 

2.3 

1.4 

9801 

1.6 

1.5, 

incomplete 

9802 

1.6 

0.84, 

incomplete 

-- 

9803 

0.94 

1.1 

1.2 

Open -cell 

Aluminum 
(1.4  g/cc ) 

9930 

1.1 

1.1 

1.0 

9932 

1.6 

1.2, 

incomplete 

-- 

C.l  osed-cel  1 

Aluminum 
(1.4  g/cr. ) 

9892 

1.1 

1.05 

1.0 

9894 

1.6 

1.9 

1.2 

Si  1 ica 
(1.0  g/cc) 

9891 

0.72 

0.1, 

incomplete 

-- 

9893 

1.1 

0. 6i . 

incomplete 

-- 

9931 

1.1 

1.2 

1.1 

9345  (  1-D) 

0.68 

0.52, 

incomplete 

-- 
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Table  4  cone  I  tided 


MATERIAL,  NOMINAL 
DENSITY.  SHOT  NO. 

INPUT 
MOMENTUM 
DENSITY.  Iq 

(  1  fl4  taps) 

IMPULSE  IN  SMEAR 
CAMERA  RECORD,  IR 
(104  taps) 

IMPULSE 

GAIN 

!r 

io 

Graphite. 

(ATJ) 

(1.7  g  cc ) 

10082 

1.0 

1.2 

1.2 

10084 

1.0 

1.3 

1.3 

Graphi te, 

(IT  0114) 

(1.0  g/ cc ) 

10035 

1.0 

1.2 

1.2 

10036 

1.0 

1.2 

1.2 

10037 

0.66 

0.58,  incomplete 

-- 

Beryl  1 ium. 

(  I  .  1  g,  cc 

10446  (standoff) 

1.12 

1.4 

1.3 

10453  (standoff) 

0.64 

0.56 

0.87 

Average 

1.0 

Notes : 

* 

Iq  =  (t/0.016)  104  taps  where  t  is  explosive  thickness  in 
inches  for  all  shots  except  those  marked  standoff  and  1-D. 

^  When  the  impulse  in  smear  camera  record,  Ip,  is  marked 
“incomplete"  the  record  was  broken  off  by  reverberations  in 
2-D  shots  or  by  driver  arrival  in  1-D  experiments  before 
pressure  returned  to  zero  at  the  foam-anvil  interface. 

§  For  1-D  shots  Iq  =  (t  / 0 . 02  0 )  1.13  X  104  taps  where  t  is 

flyer  thickness  in  irches.  Descriptions  of  these  experiments 
and  of  Shot  8863  appear  in  Rempe  1  . ‘ 

^  When  standoff  used  Iq  is  calculated  from  the  Gurney  formula 
bv  treating  flyer  as  a  free  plate. 


EL-506D  sheet  explosive  to  input  momentum  density  JQ  in  two  di f fe rent  way s, 
depending  upon  whether  the  explosive  impact  drives  the  aluminum  flyer  alone 
and  the  impact  on  the  foam  is  secondary,  or  whether  barrier  and  foam  are 
arranged  in  contact  with  each  other  before  detonation.  Unfortunately, 
the  camera  records  of  the  pressure  histories  in  the  1-D  experiments,  in 
which  input  momentum  is  more  accurately  known  than  in  the  2-D  shots,  were 
interrupted  before  pressure  returned  to  zero  at  the  foam-anvil  interface; 
so  all  our  observations  of  impulse  gain  depend  on  uncertain  values  for 
input  momentum  The  table  does  furnish  evidence,  however,  that  the 

momentum  gain,  the  ratio  of  observed  impulse  in  the  anvil  / fl  to  7Q,  is 
lower  for  1.4  g/cc  aluminum,  both  open-  and  closed-cell,  than  for  any 
other  material  except  possibly  silica,  and  that  the  gain  for  polyurethane 
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is  higher.  If  the  foamed  aluminum  and  beryllium  lock  to  crystal  density, 
it  is  certainly  reasonable  that  their  shock  impedances  are  the  nearest  to 
that  of  steel  of  the  materials  tested.  Locked  silica  also  probably  has 
a  relatively  high  shock  impedance.  Taking  shock  impedances  and  m2of 
impacting  plates  as  constants,  the  momentum  gain  in  the  forward  direction 
is  2m2/'(mj  +  m2)  if  m  2  is  the  impedance  of  the  plate  initially  at  rest. 

It  is  clear  that  our  results  for  gain  generally  vary  among  themselves  in 
the  way  expected  from  this  formula. 

Just  as  we  related  certain  features  of  the  forerunner  to  quasi-static 
compression  data  by  a  theory  of  double  shock  structure,  we  seek  bases  for 
forecasting  such  important  characteristics  of  the  main  wave  as  peak  pres¬ 
sure  and  total  impulse  or  wave  shape.  We  have  developed  three  theoretical 
tools  for  this  task:  a  simplified  characteristics  method,1,2  a  more  com¬ 
plete  characteristics  method,  and  the  (Q)  method  of  artificial  viscosity.  12,13 
As  data  for  use  with  these  tools  we  have  the  value  of  the  unfoamed  den¬ 
sity  of  the  material  constituting  the  foam  and  the  values  of  final  strain 
in  quasi-static  compression.  As  a  check  on  these  values  foam  samples 
under  shock  compression  were  flash  X-rayed,  but  due  to  edge  effects  both 
in  the  compression  and  in  the  X-ray  beam  passing  through  the  locked  mass, 
the  results  have  so  far  been  difficult  to  interpret  exactly,  although 
they  clearly  show  compression  to  near  undistended  density  behind  the 
main  wave.1,2 

For  several  foams  the  simplified  characteristics  theory  predicts 
approximately  the  greatest  value  of  the  quantity  IQ/pQl  cor respondi ng  to 
the  transmission  to  a  wall  of  the  forerunner  without  the  main  wave.  We 

calculate  the  particular  value  5  of  lQ/p0l  corresponding  to  collision  of 

* 

the  locked  wall  reflection  and  the  oncoming  locked  wave  at  the  time  t 
when  pressure  and  particle  speed  jump  in  the  main  wave  disappears:1* 


In  this  formula  V'Q  =  initial  specific  volume,  ~  locked  specific  volume, 
V't  *  specific  volume  in  the  forerunner,  and  l't  =  forerunner  speed.  In 
theory  when  Iq/^qI  =  S,  there  is  a  second  pressure  pulse  delivered  to  the 
wall,  but  the  rate  of  reduction  of  peak  pressure  with  foam  thickness  in 
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the  neighborhood  of  i  Q  /  /Q  Q  Z  =  5  is  high  and  experimentally  it  is  a  good 
estimate  of  conditions  needed  to  avoid  the  main  wave. 

Except  for  premature  but  unimportant  interruption  of  the  record, 

1-D  Shot  9 1 5 5 2  produced  an  excellent  pressure  history  (Fig.  6)  and  we 
can  use  our  study  of  it  to  illustrate  many  general  features  .of  our  find¬ 
ings.  First,  we  chose  an  elastic  pressure-volume  locus  of  the  slope  and 
peak  pressure  to  give  the  observed  forerunner  stress  in  the  anvil  (1  kbar) 
and  the  observed  forerunner  speed  (1.65  mm//xsec),  and  to  give  the  original 
foam  specific  volume  1.515  cm3/g.  Next  by  calculations  with  the  simple 
theory  we  found  a  locked  specific  volume  which  predicted  the  observed 
time  interval  between  first  and  second  wavefronts  (i.e.,  1.7  to  2.6mRP() 
a  good  compromise  value  is  1.0  cmVg.15  Next  we  chose  a  value  for  wave 
speed  in  locked  foam  which  would  give  the  observed  value  of  the  time  in¬ 
terval  between  the  main  (B)  and  reverberation  (C)  wave;  this  value  can 
be  2.11  mm/^sec.  By  assuming  the  locked  equation  of  state  is  a  straight 
line  passing  through  the  point  P  =  1  kbar  and  V  ~  1.00  cc/g  with  slope 
related  to  assumed  shock  speed  in  the  locked  mass,  we  then  had  an  energy- 
independent  equation  of  state  (Fig.  7).  (Actually  for  umimportant  reasons 
the  slightly  curved,  solid  line  in  Fig.  7  was  used  in  the  Q-method  calcu¬ 
lations;  the  straight,  dotted  line  described  in  the  text  was  used  in  the 
computations  by  characteristics.)  Using  this  eqiation  of  state  for  the 
foam  and  a  simple  straight  line  P- 1  equation  of  state  for  the  aluminum 
in  the  flyer,  we  carried  out  calculations  graphically  by  the  character¬ 
istics  method  and  electronically  by  the  Q-method,  Pressure  and  particle 
speed  distributions  resulting  from  these  two  methods  at  time  2.47  m sec 
after  flyer  impact  are  plotted  in  Figs.  8  and  9.  Since  the  characteristics 
and  artificial  viscosity  solutions  agreed  closely,  we  used  the  artificial 
viscosity  method  to  compute  the  pressure  history  in  a  rigid  wall  behind 
5  mm  of  the  foam  whose  equation  of  state  is  represented  in  Fig.  7.  This 
history  is  shown  in  Fig.  6  alongside  the  observed  history. 

There  are  a  number  of  striking  and  obvious  discrepancies  between  the 
two  curves  of  Fig.  6.  About  10  percent  of  the  peak  pressure  by  the 
Q-method  can  be  attributed  to  the  false  assumption  of  rigidity  in  the 
wall.  To  get  better  agreement  shock  speed  in  the  locked  material  would 
have  to  be  lowered  to  about  0.9  mm/^sec  which  would,  of  course,  separate 
the  B  and  C  fronts  far  beyond  the  observed  relation  unless  we  attribute 
to  the  reverberation  front  a  much  higher  speed  during  its  return  from  the 
flyer  than  during  its  trip  from  wall  to  flyer.  This  is  quite  reasonable 
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"SSURE,  P  - -  k  bar 


FIG.  8  PRESSURE  PROFILE  IN  A  FOAM  ACCORDING  TO  METHODS  OF 
CHARACTERISTICS  AND  OF  ARTIFICIAL  VISCOSITY 
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FIG.  9  DISTRIBUTION  OF  PARTICLE  SPEED  IN  A  FOAM  ACCORDING  TO 
METHODS  OF  CHARACTERISTICS  AND  ARTIFICIAL  VISCOSITY 
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but  would  involve  a  serious  modification  of  the  equation  of  state.  The 
relaxation  behind  the  C  iront  seen  in  Shot  9155  is  much  slower  than  that 
calculated,  but  this  discrepancy  could  be  removed  by  alteration  of  the 
X  equation  of  state  to  provide  gradual  isentropic  relaxation  from  the  locked 

f  state  (P  =  1  kbar,  V  =  1.0  cc/g)  to  some  state  near  (P  =  0,  V  =  1.515  cc/g). 

Neither  the  characteristics  method  nor  the  Q-method  is  expected  to  give 
information  on  shock  width,  of  course,  but  the  rise  time  in  the  main  wave 
is  unusually  long  compared  to  ordinary  materials,  and  the  pressure  plateau 
behind  the  main  wave  is  completely  unexpected.  Since  waves  in  the  flyer 
make  between  three  and  ten  round  trips  each  microsecond,  there  should  be 
a  pressure  gradient  across  the  locked  foam  layer,  and  we  expect  as  a  re¬ 
sult  to  see  a  gradient  behind  the  main  wave  in  the  anvil. 

The  pressure  plateau  is  frequently  but  not  always  found  in  experi¬ 
ments  involving  polyurethane,  silica,  and  aluminum  but  never  in  graphite 
where  steady  pressure  relaxation  is  the  rule,  The  phenomenon  is  partic¬ 
ularly  striking  in  Shot  9931  in  silica.  Shot  9803  in  aluminum.  Shot  10,446 
in  beryllium,  and  Shot  9832  in  polyurethane  (Fig.  3).  Both  the  unusually 
long  rise  time  in  the  main  front  and  the  plateau  must  be  connected  with 
the  porous  nature  of  the  material,  although  there  is  the  possibility  in 
some  experiments  of  an  interaction  between  main  wave  and  a  simple  com¬ 
pression  wave  moving  back  into  the  foam  from  the  wall  after  forerunner 
arrival.  This  would  give  the  appearance  of  a  thick  shock  front.  (Shots 
9800,  9801,  and  9809  in  aluminum  showed  slow  rises  to  main  wave  pressure 
in  the  absence  of  a  forerunner.)  Since  graphite  is  porous  but  shows  no 
main  wave  plateau,  there  may  be  another  factor  or  factors  contributing  to 
the  flatness. 

Since  there  are  pores  against  the  anvil  and  since  the  full  pressure 
in  a  wave  cannot  be  passed  into  the  anvil  until  these  are  closed,  rise 
time  may  be  simply  this  closing  time.  Pores  in  0.67  g/cc  polyurethane 
measured  microscopically  are  0.3  to  0.6  mm  in  diameter;  since  locking 
wave  rise  time  in  Shot  9155  is  1.5  ^.scc,  a  0.2  mm/Vsec  particle  speed 
would  be  needed  to  close  all  the  pores  during  rise  time.  This  is  indeed 
the  particle  speed  expected  behind  the  main  wave  at  the  time  of  its  ar¬ 
rival  at  the  anvil,  see  Fig.  9  for  example. 

As  the  main  wave  becomes  stronger  its  plateau  character  tends  to  be¬ 
come  less  pronounced  (Shots  9894,  9809,  9800,  and  9932  in  Fig.  3).  This 
is  reasonable  because  the  lack  in  uniformity  in  porous  material  becomes 
less  important  as  shock  and  particle  speeds  become  more  nearly  alike. 
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By  comparing  the  results  of  a  few  very  simple  calculations  with 
observed  peak  pressures,  it  can  be  shown  that  the  so-called  "locked 
states"  in  aluminum  and  beryllium,  as  in  polyuret  liane ,  do  noi  have  the 
same  shock  impedances  as  do  the  same  materials  in  undistended  form  Th'* 
disparity  appears  to  be  greater  for  light-weight  foams  than  for  the  more 
dense  If  the  impedance  of  both  solid  aluminum  and  beryllium  is  assumed 
to  be  160  x  106  g  cm  2  sec  ‘  and  the  impedance  of  steel  is  2.25  x  10h 
g  cm  2  sec  1 ,  and  if  it  is  noted  that  all  the  material  (l nil uding  flyer) 
behind  the  locking  front  moves  at  nearly  the  same  speed,  we  ran  calculate 
a  lower  bound  to  the  pressure  expected  upon  meeting  of  steel  wall  and  a 
main  wave  front  behind  which  the  aluminum  or  beryllium  is  in  the  normal 
undistended  state.  In  cases  where  there  is  a  forerunner  ahead  of  the 
main  front  we  subtract  from  the  total  input  momentum  the  momentum  passing 
into  the  steel  before  main  wave  arrival  and  in  computing  the  mass  of 
locked  material  we  neglect  the  mass  piled  up  at  the  wall  by  the  fore¬ 
runner.  To  estimate  the  lower  bound  we  also  neglect  the  pressure  behind 
the  main  front  just  before  impact.  Comparisons  of  these  calculated  lower 
bounds  with  experimental  observations  are  shown  below: 
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The  fourth  column  above  contains  values  of  shock  impedance  computed  from 
the  observed  interface  pressure:*  tin  the  third  «olumn)  and  the  foregoing 
assumptions  Actual  impedances  in  the  lork.-d  foam  should  be  lower  than 
entries  in  the  fourth  col umn 

Graphif,  especially  the  ATJ .  appears  to  be  more  nrailv  ttulv 
"locked”  behind  th*  ma  i  wave  than  am  of  the  above  materials 
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MATERIAL 

ESTIMATED 
LOWER  BOUND 
ON  PRESSURE 

(kbar ) 

OBSERVED 

PEAK 

PRESSURE 
<  kbar ) 

UPPER  BOUND  ON 
OBSERVED  IMPEDANCE 

(10®  g  c»  2  aee  *) 

SHOT  NO. 

1.1  g/cc  C-PT0114 

6.7 

3.9 

0.29 

10.035 

1.1  g/cc  C-PT0114 

4.3 

2.0 

0.22 

10,036 

1.1  g/cc  C-PT0114 

2.7 

1.3 

0.25 

10,037 

1.7  g/cc  C-ATJ 

5.75 

5.7 

1.6 

10,084 

1.7  g/cc  C-ATJ 

3.5 

2.6 

0.39 

10,082 

Estimates  above  are  based  on  an  assumed  undistended  impedance  for  carbon 
of  0.55  *  10®  g  cm”2  sec-1.  This  better  agreement  between  estimated  and 
observed  pressures  may  be  connected  with  the  more  conventional  wave  shapes 
seen  in  graphite  and  noted  earlier  in  this  report. 

Choice  of  an  assumed  "locked”  impedance  for  1.1  g/cc  silica  is  not 
so  clear  but  if  we  take  that  of  a  quartzite,  1.63  x  10®  g  cm”2  sec”1,  we 
find: 


ESTIMATED 
LOWER  BOUND 
ON  PRESSURE 
(kbar) 

OBSERVED 

PEAK 

PRESSURE 

(kbar) 

UPPER  BOUND  ON 
OBSERVED  IMPEDANCE 

(10®  |  e«"J  aae”1) 

SHOT  NO. 

1.8 

7.0 

0.445 

9,325 

8.8 

2.4 

0.295 

9,893 

7.5 

2.2 

0.33 

9,345 

C.  THE  BEST  FOAM 

In  our  present  terms,  the  best  foam  holds  the  peak  pressure  on  a 
structure  stemming  from  an  impact  carrying  a  certain  momentum  density 
below  a  certain  limit  by  adding  the  least  mass  per  unit  area  to  the 
structure.  The  choice  of  the  best  foam  may  also  depend  on  environmental 
conditions,  the  properties  of  the  protected  structure,  and  other  factors 
not  considered  here. 

When  momentum  is  deposited  in  the  foam  layer  instantaneously,  the 
full  forerunner  stress  (greater  than  the  characteristic  pressure  Pe  be¬ 
cause  of  reflection)  will  be  felt  in  the  structure.  In  other  cases  it 
may  be  possible  to  hold  the  peak  stress  below  this  level  but  it  is  not 
likely  these  cases  will  be  of  practical  interest  when  protection  by  the 
use  of  elastic-rigid  foams  is  sought.  The  least  stress  that  can  be  ex¬ 
pected  then  is  what  we  have  called  forerunner  pressure  in  the  anvil. 
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Of  the  foams  studied,  polyurethane  is  the  best  in  the  sense  of  the 
above  definition.  Its  ability  to  stop  the  main  wave  is  the  highest,  and 
the  rate  of  increase  of  peak  pressure  with  increasing  I0/p0l  lowest 
(Fig.  5).  Judging  from  the  relative  independence  of  I0/P0f  on  foam  dis¬ 
tention  seen  in  Fig.  4,  the  demands  of  a  structural  tolerance  lower  than 
the  forerunner  stress  of  a  given  density  of  polyurethane  would  better  be 
met  by  increasing  the  distention  of  polyurethane  instead  of  substituting 
a  material  such  as  graphite  with  a  lower  forerunner  pressure.  The  effect 
of  distention  on  foam  effectiveness  is  less  than  the  effect  of  the  kind 
of  material  (within  our  range  of  materials).  Hie  single  data  point  from 
0.33  g/cc  polyurethane  (Shot  9228 ) 2  indicates  that  the  effect  of  disten¬ 
tion  may  be  greater  in  polyurethane  than  in  aluminum,  but  in  our  exper¬ 
ience  it  is  true  that  the  influence  of  distention  on  foam  quality  is  less 
than  the  effect  of  kind  of  material.  Polyurethane,  however,  has  a  greater 
tendency  to  rebound  (thus  increasing  momentum  delivered  to  structure) 
than  have  materials  of  higher  shock  impedance. 

We  have  had  fair  success  roughly  judging  peak  pressures  in  widely 
different  kinds  of  foams  from  the  calculated  values  of  the  quantity  5. 
This  success  may  possibly  become  greater  as  more  accurate  values  of  the 
parameters  upon  which  S  depends  become  known.  In  Table  5  appear  the  cal¬ 
culated  values  of  5  and  the  observed  values  of  f0/P0I  *  5'  separating  1- 
and  2-wave  shock  systems.  Also  listed  are  the  parameter  values  used  in 
the  computation  of  S.  This  scheme  is  useful  only  in  locating  the  general 
area  in  a  plot  such  as  Fig.  5  in  which  to  expect  observed  peak  pressures. 


Table  5 

CALCULATED  AND  OBSERVED  STOTPING  POWER  OF  FOAMS 


MATERIAL 

U. 

(■■/(in e  ) 

p. 

(kb.r) 

vo 

(ee/f  ) 

V1 

(ee/f) 

S 

(10*  cm/aac) 

S' 

(10*  ea/aae) 

Polyurethane 

m 

1.515 

1.00 

2.2 

1.6  -  2.4 

Aluminum 

(open-cel  1 ) 

i 

1.05 

0.37 

1. 15 

*-1.6 

Aluminum 

( c  losed-ce 1 1 ) 

2.3 

0.715 

0.37 

2. 1 

1.95  -  2.7 

Silica 

1.35 

0.65 

0.95 

0.50 

2.0 

1.25  -  1.35 

Graphite  (ATJ) 

1.3 

0.71 

0.58 

0.50 

0.33 

0.75(’> 

Graphi te 
(PT0114) 

1.0 

0.70 

0.87 

0.50 

0.85 

<0.70 

Beryllium 

2.0 

1.5 

0.83 

0.55 

1.3 

<0.8 
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D.  QI'ASI  -  STATIC  COMPRESSION  OF  FOAMS 


So  v  oral  cylinders  of  eac  li  of  the  foams  studied  were  carefully  dimen¬ 
sioned  and  weighed,  encased  in  hollow  steel  cylinders,  and  compressed 
hydraulically  at  a  rate  of  6,000  psi/min  At  ordinary  temperatures  stress 
is  linearly  related  to  strain  below  a  certain  limit  of  proportionality. 

At  higher  stress  the  rate  of  change  of  stress  with  strain  becomes  much 
smaller  until  a  point  is  reached  at  which  this  rate  becomes  zero  (point 
of  inflection).  As  stress  rises  still  higher,  the  rate  increases  again. 
The  compression  curve  in  Fig.  10  shows  results  for  polyurethane. 

Data  reduced  from  quasi-static  testing  appears  in  Table  6-  The  pre¬ 
dicted  elastic  wave  speed  is  the  square  root  of  the  ratio  of  the  moduius 
of  linear  compression  to  density.  Some  foam  specimens  were  heated  and 
compressed  at  temperatures  above  ambient. 

E.  EFFECT  OF  HEAT 

Although  our  experiments  with  preheated  foams  have  been  too  few  for 
close  calculation  of  the  amount  of  change,  it  is  clearly  established  that 
heat  tends  to  reduce  elastic  modulus,  elastic  yield,  and  forerunner 
strength.  One  experiment  with  preheated  closed. cell  aluminum  (1.4  g/cc) 
shows  a  considerably  larger  lowering  of  this  strength  than  might  be  fore¬ 
cast  from  the  weakening  of  the  static  yield  strength;  the  forerunner 
stress  in  the  anvil  sterns  to  be  reduced  more  than  50  percent  under  pre¬ 
heating  from  ambient  to  318°C,  and  the  static  proportional  limit  falls 
only  30  percent  under  heating  to  426°C.16  In  fact,  the  first  wave  stress 
in  the  anvil  is  quite  near  the  static  yield  of  the  heated  specimen.  As 
a  result,  the  value  of  S  for  this  foam  is  reduced  by  preheating.  Since 
the  measured  peak  pressure  is  higher  as  a  result  of  preheating  it  is  also 
likely  the  experimentally  defined  stopping  power  S'  of  the  foam  is  lowered 
that  is,  if  the  more  or  less  linear  relation  between  peak  pressure  in  the 
anvil  and  I0/IPq  continues  to  hold  in  the  preheated  material,  preheating 
displaces  this  line  leftward  toward  lower  values  of  /0/V0I. 

Tli e  preheated  foam  also  shows  stress  relaxation  in  the  forerunner 
similar  to  the  relaxation  seen  in  the  un heated  samples.  Observation  of 
first  wave  speed  in  the  preheated  aluminum  failed. 

Preheating  polyurethane  does  not  seem  to  have  so  strong  an  effect 
on  shock  behavior  as  expected  from  static  observations;  the  decline  in 
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Table  6 


QUASI- STATIC  ONE-DIMENSIONAL  COMPRESSION  OF  FOAMS 


MATERIAL 

INITIAL 

DENSITY 

(g/cc) 

TEST 

TEkfERATURE 

(°C> 

ELASTIC 

MODULUS 

(kbar) 

PREDICTED 

El.  A  STIC 

WAVE  SPEED 
(mm/fiaec  ) 

YIELD 

STRESS 

0, 2%  OFFSET 
(kbar) 

COMPRESSION 
AT  1.5  kb 

(%) 

MODULUS 

AT  1.5  kb 
(kbar) 

Polyurethane 

0.660 

Ambient 

11.9 

1.34  ±  0.07 

0.331  ±  0.007 

41.8 

12.3 

(closed-cel  1, 

G.  664 

Ambient 

12.9 

1.40  t  0.08 

0.304  ±  0.007 

41 

Not  measured 

Po  1  y  t  ron ) 

0.659 

Ambient 

16.4 

1.58  i  0.08 

0.310  ±  0.007 

41 

Not  measured 

0.661 

Ambient 

12.3 

1.37  i  0.07 

0.296  ±  0.007 

41 

Not  measured 

0.  654 

Ambient 

12.0 

1.35  ±  0.07 

0.331  ±  0.007 

43 

Not  measured 

0.900 

Ambient 

21.6 

1.  58  +  0.08 

0.807  ±  0.007 

17.4 

0.711 

0.667 

120 

4.57 

nifltLmji'im 

47.5 

32.4 

0.668 

157 

-  Negligible  ► 

S 1  *  5 

30.8 

0.893 

120 

KSSfl 

27.3 

12.5 

0.880 

200 

31.4 

37.6 

Polyurethane 

mm 

.Ambient 

1.50  ±  0,08 

0.551  ±  0.007 

27.8 

— 

(closed  cell. 

120 

HESS 

0.905  i  0.045 

0.151  ±  0.007 

45.2 

CPP) 

HU 

mm 

HSU 

n  cc-3 

.  1  • 

XI  11 

0. 01 72  +  0.004 

68.0 

20.8 

A 1  uni n  un 
(open  cell) 

Ambient 

6.61 

0.936  ±  0.047 

0.0276  i  0.007 

62.0 

17.9 

0.946 

Amb i en  t 

6.51 

0.  326  ±  0.041 

0.0690  ±  C. 007 

55.0 

8.96 

1.30 

Ambient 

56.5 

2.08  ±  0.10 

0.  345  ±  0.007 

34.5 

8*  68 

1.30 

200 

23.6 

1.37  i  0.07 

0.248  ±  0.007 

36.8 

10.9 

1.30 

426 

16.9 

1.14  ±  0.06 

0. 186  ±  0.007 

39.0 

8.82 

Aluninun 

1.41 

Ambi ent 

mm 

2.06  ±  0.10 

0.759  i  0.007 

22.2 

83.4 

(closed  cell ) 

1.43 

426 

ESI 

1.11  t  0.05 

0.690  ±  0.007 

27.6 

41.9 

Silica 

1.02 

Ambient 

10.55 

1.01  i  0.05 

0.165  ±  0.007 

36.7 

17.7 

(first  lot) 

Silica 

0.585 

Ambient 

Not 

measured 

0.0552  +  0.007 

60.5 

27.0 

(second  lot) 

0.  931 

Ambient 

15.7 

1.30  ±  0.065 

0.0965  i  0.007 

42.5 

17.65 

0.945 

426 

7.86 

0.  915  i  0.045 

0.206  t  0.007 

28.0 

13.0 

Graphite 

0.975 

Ambient 

2.20 

0.142  ±  0.013 

0.0482  ±  0.007 

50.0 

35.8 

(electrode) 

Graphite 

1.18 

Amb  i  en  t. 

13.1 

1.02  ±  0.092 

0.0551  ±  0.007 

21.2 

14.3 

(PT-0114) 

Graphite 

m 

Ambient 

■na 

1.32  i  0.13 

0.497  i  0.007 

mm 

11.8 

(  A  i J  ) 

gB 

426 

m 

0.905  ±  0.090 

0.504  ±  0.007 

wtam 

9.69 

Beryl  1 ium 

UBS 

Ambient 

56.5 

2.175  ±  0.36 

0.276  i  0.0C7 

■  ' 

EH 

Ambient 

50.9 

2.045  ±  0.04 

0.641  i  0.007 

BUSES1 

NOTES: 

All  specimen*  compressed  into  tight-fitting  steel  dies- 

Uncertsinty  limits  sbove  refer  to  a  particular  measurement  and  do  not  include  the  likely 
variation  between  samples- 

All  densities  measured  at  ambient  temperature- 


static  strength  is  near  70  percent  at  120°C,  but  the  forerunner  is  weak¬ 
ened  by  less  than  40  percent  in  two  experiments  with  preheated  material.1. 
In  these  experiments  main  wave  stress  is  raised  from  about  0.7  to  3.0kbar 
by  preheating,  and  the  peak  pressure  us.  /q/P01  relation  is  displaced  to 
the  left  in  Fig.  5.  Although  the  effect  of  preheating  on  elastic  modulus 
at  120°C  was  strong,  our  investigation  of  the  influence  of  preheating  on 
forerunner  speed  was  inconclusive. 

Preheating  slows  the  forerunner  pressure  rise,  and  the  effect  was 
more  dramatic  in  our  experiment  with  closed-cell  aluminum  than  in  our 
polyurethane  shots. 


* 
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3.  MAJOR  EXPERIMENTS 


A.  DESCRIPTION  OF  TWO-DIMENSIONAL  TECHNIQUES 

In  all  experiments  a  detonation  front  in  sheet  explosive  moving 
normal  to  itself  generated  a  line  impact  traveling  across  one  face  of  a 
thin  slab  of  foam;  the  waves  in  the  foam  stemming  from  this  impact  passed 
into  a  hardened  steel  anvil  in  close  contact  with  the  other  side  of  the 
slab.  The  resultant  angular  displacement  of  images  reflected  in  a  mirror 
free  surface  on  the  anvil  was  photographed  with  a  high-speed  smear  camera. 
Because  the  waves  meet  the  interface  and  the  free  surface  obliquely  the 
technique2  is  called  two-dimensional  (2-D).  Figure  1(a)  illustrates  both 

the  general  arrangement  and  several  special  features  such  as  flyer  tilt, 
anvil  bevel,  standoff,  and  foam  taper  which  were  used  in  various  experi¬ 
ments.  (The  figure  shows  the  barrier  on  flyer  plate  tilted  with  respect 
to  the  foam  free  surface.  When  sound  speed  in  the  foam  is  close  to  or 
above  detonation  speed,  this  tilt  must  be  used  to  achieve  impact  speed 
higher  than  sound  speed.  This  feature  has  not  yet  been  needed  in  our 
work.)  Usually  the  surfaces  and  interfaces  were  parallel  to  each  other. 
When  the  foam  pores  are  connected  with  each  other,  detonation  products  may 
flow  into  the  slab  to  confuse  the  response,  so  in  every  experiment  except 
the  first  two  we  put  a  thin,  solid  aluminum  barrier  (or  flyer)  between  the 
sheet  of  explosive  and  the  impacted  surface  of  the  foam.  When  the  effect 
of  preheating  foams  was  studied,  we  wanted  to  insulate  the  explosive  from 
the  hot  foam  as  much  as  possible,  so  a  3/16-inch-wide  standoff  was  kept 
between  barrier  and  foam  surfaces  by  means  of  a  thin  paper  spacer  on  edge. 
In  certain  experiments  we  did  not  make  the  mirror  surface  of  the  anvil 
parallel  to  the  impacted  surface,  but  beveled  it  to  change  the  collision 
angle  between  wavefront  and  free  surface.  As  a  means  of  measuring  wave 
speed  in  foam,  in  two  shots  we  beveled  or  tapered  the  foam  slab.  Finally, 
as  a  more  general  means  of  measuring  wave  speed  in  foam,  we  let  the  ex¬ 
plosive  shock  in  every  experiment  interrupt  reflection  in  a  glass  mirror 
at  the  end  of  the  detonation  run.  Figure  1(b)  illustrates  this  feature 
and  Table  1  gives  details  of  each  experiment.  Since  the  value  of  the 
rate  of  travel  of  the  initial  impact  along  the  foam  slab  >>  important  in 
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data  reduction,  we  frequently  put  time-of-arrival  switches  on  the  free 
surface  of  the  sheet  explosive. 

To  measure  the  angular  displacement  of  the  mirror  we  used  an  optical 
lever  arm3  generally  about  20  feet  long  with  a  40-inch  focal-length  lens 
placed  about  20  feet  from  the  mirror.  A  typical  smear  camera  record  ap¬ 
pears  in  Fig.  2  (Shot  10,331).  The  lines  are  images  of  points  on  trans¬ 
lucent  grid  lines  ruled  at  known  separations  on  opaque  glass  covering  one 
end  of  an  explosive  argon  light  source.5  The  vertical  scale  in  Fig.  2  is 
computed  from  the  grid  line  spacings  and  the  optical  distances  in  the  ex¬ 
periment,  the  horizontal  or  time  scale  comes  from  the  known  rotation 
speed  of  a  mirror  in  the  smear  camera.  Two  long,  narrow,  mutually  per¬ 
pendicular  apertures  arc  used  in  tho  optical  system:  the  larger  is 
usually  5/16  inch  wide  and  stands  between  lens  and  mirror  to  limit  the 
reflection  zone  of  a  point  in  the  source  to  a  strip  about  4  mm  wide  in 
the  direction  perpendicular  to  the  run  of  the  detonation;  the  smaller 
aperture,  0.05  mm  wide,  lies  in  the  focal  plane  to  pass  on  to  the  film 
only  light  from  essentially  one  point  in  each  grid  line.  The  mirror  or 
polished  area  of  the  anvil  face  is  itself  only  3/4  inch  wide  in  a  direc¬ 
tion  normal  to  the  detonation  travel  but  extends  the  whole  length  of  the 
anvil  in  the  direction  parallel  to  the  detonation  run. 

Thus  we  see  in  the  lines  of  Fig.  2  the  result  of  a  train  of  wave- 
fronts  sweeping  across  the  steel  mirror  and  moving  successively  into  re¬ 
flection  zones  for  successive  grid  lines.  The  speed  of  this  motion  in 
terms  of  distance  on  the  mirror  we  call  in  this  report  “apparent  wave 
speed.’’  To  each  effective  light-source  point  and  for  each  wavefront  at 
least  two  images  correspond,  one  from  the  undisturbed  reflection  zone 
and  the  second  from  the  tilted  zone  at  the  end  of  the  angular  motion. 

If  the  rise  time  in  the  front  were  zero,  there  would  be  only  two  images. 
Because  rise  time  is  usually  quite  an  appreciable  part  of  a  microsecond, 
the  sharp  separation  between  these  two  images  is  not  often  seen.  Judging 
from  the  duration*  of  intermediate  images  in  Fig.  2,  the  rise  time  of  the 


*  Tin*  interval*  between  event*  in  *  etve  trein  «u*t  be  aeaaured  froa  horiiontel  diatancea  in  Fig.  2  be¬ 
tween  eloping  line*  dr*«n  through  iaage*  a*rking  thoae  event*.  After  th*  airror  »urf*ce  i*  tilted  the 
reflection  tone  of  •  given  line  i*  noved  i-  the  direction  oppoiite  to  the  detonation  run;  becauae  of 
the  2-D  geometry  thia  dieplaceaent  iapliea  >  change  in  origin  of  the  tiae  acale.  By  drawing  atraight 
linea  through  iaagea  of  correaponding  event*  falling  in  auiceeaive  grid  point  loci  we  can  correct  for 
thia  ahift  in  origin.  When  aeaaured  in  the  coordinate  ayatea  ahown  in  Fig.  2,  th*  alope  of  all  line* 
joining  correaponding  evant-iaage*  ia  th*  apparent  wav*  apeed. 

i 

'  Th*  tiae  aeparation  between  th*  firat  appearance  of  on*  and  the  laat  of  th*  other  in  thia  ideal  caae 
depend*  on  th*  width  of  a  reflection  ion*,  th*  aeparetion  of  the  two  reflection  tone*,  the  apparent 
apeed  of  th*  wev*  along  th*  airror,  fila  aenaitivity,  and  aaount  of  light. 
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FIG.  11  SMEAR  CAMERA  RECORD  SHOWING  FAST-RISING  FORERUNNER  — 

SHOT  9832 

forerunner  or  A-wave  is  at  least  0.67  y sec,  and  of  the  main  or  B- wave  at 
least  1.5  y. sec.  Presumably  the  forerunner  rise  time  is  much  less  in  Fig.  11 
(Shot  9,832)  where  the  A-wave  images  are  much  more  clearly  separated. 

Most  likely  the  C-wave  in  Fig.  2  stems  from  the  double  reflection 
of  the  main  wave  off  the  anvil  and  again  off  the  barrier.  The  following 
disturbances,  A'  and  B  ,  are  reverberations  of  the  first  and  second  waves 
entirely  within  the  anvil.  One  part  of  the  light-source  cover  is  not 
ruled  but  is  completely  translucent  Light  passing  through  there  is  re¬ 
flected  in  the  glass  mirror  abutting  the  head  of  the  anvil  and,  as  is 
seen  above  the  grid  lines  in  Fig.  2,  exposes  the  top  of  each  film  until 
the  detonation  of  the  explosive  lying  against  the  glass  mirror  destroys 
its  silver  surface.  The  distance  D  on  the  anvil  is  read  off  the  cali¬ 
brated  film  record  and  is  used  to  calculate  average  wave  speed  in  the 
foam  as  illustrated  in  Fig.  1(b)  We  have  measured  sound  speed  at  sev¬ 
eral  temperatures  in  our  hardened  steel  anvils  (5.88  ±  0.01  mm/Msec),  and 
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in  this  calculation  we  assume  that  all  waves  in  the  anvils  move  at  that 
speed.  In  all  experiments  using  parallel  surfaces,  t.e..  those  not  using 
tapered  slabs,  beveled  anvils,  or  tilted  flyers,  the  apparent  wave  speed 
is  assumed  to  be  detonation  speed,  i.e.,  7.35  ±  0.05  ram'^ser  at  ambient 
temperatures,  unless  the  apparent  speed  seen  in  the  record  does  not  agree 
with  this  value  within  experimental  error,  in  which  case  the  apparent 
speed  seen  in  the  camera  record  is  used. 

The  interval  4-4'(Fig.  2)  generally  marks  the  observation  time 
available  in  this  method  and  depends  on  the  anvil  thickness.  Useful  ex¬ 
periment  time  is,  of  course,  also  limited  by  the  incoming  rarefactions 
from  the  sides  of  the  foam  slab  and  of  the  anvil  which  parallel  the  det¬ 
onation  run;  in  the  experiment  represented  by  Fig.  2,  A-A'  is  about 
8.7  /xsec  and  the  limit  to  useful  time  after  ,4-wave  arrival  set  by  side 
rarefactions  is  about  10.5  M sec. 

If  rarefactions  from  the  ends  of  the  anvil  or  slab  are  important, 
in  principle  the  camera  records  a  change  in  either  or  both  apparent  wave 
speed  and  stress  in  the  same  wave  while  detonation  or  impact  is  going  on, 
although  if  the  effect  is  strong  at  every  point  within  camera  view  this 
visible  change  may  be  slight.  There  are  a  number  of  stages  in  designing 
an  experiment  where  these  ending  and  starting  edge  effects  can  be  avoided 
To  establish  the  stability  of  the  first  wave  in  the  anvil  is  simple:  no 
ending  edge  effect  is  possible  and  the  anvil  is  made  at  least  a  cot 
long,  where  a  is  the  anvil  thickness,  cos  -J  =  c  D,  c  is  sound  speed  in 
the  anvil  and  D  is  impact  speed  on  the  anvil.  A  wave  moving  at  steady 
speed  D  along  the  length  of  the  anvil  will  be  stable  everywhere  in  the 
anvil  mirror  beyond  a  distance  a  cot  P  from  the  starting  edge  There  is 
no  ending  edge  effect  for  the  elastic  wave  in  the  foam  and  the  starting 
edge  release  waves  are  outrun  in  a  distance  a1  rot  P'  where  a'  is  foam 
thickness,  cos  O'  -  Uf/D‘ ,  and  D ‘  is  the  impact  speed  along  the  length 
of  the  foam.  After  stability,  of  course,  D-D  ■  Thus  in  our  experi¬ 
ments  the  forerunner  should  be  stable  after  a  run  equal  to 
(a  cot  P  *  a'  cot  P' )  along  the  mirror.  In  an  anvil  2  inches  thick  and 
with  0.67  g/cc  polyurethane  5  mm  thick  impacted  by  explosive  lying 
against  it,  this  distance  is  2.7  inches 

To  discuss  stability  in  waves  trailing  the  forerunner,  we  will 
imagine  an  impact  moving  at  speed  D  along  the  length  of  a  slab  of  thick¬ 
ness  a  ,  wave  speed  of  interest  will  be  l1  but  release  wave  speed  will 
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be  in  general  B" 


Outrunning  begins  at  a  time  t  () 


alter  star'  o I 


l  ;;.[>a<  t 


Also  proportional  to  a"  is  the  corresponding  distance  h'  along  the  unim- 
paeted  surface: 


Figure  12  illustrates  these  quantities. 


FIG.  12  OUTRUNNING  IN  A  SLAB 


A  point  a  distance  Z"  From  the  starting  edge  of  the  mirror  will  be 
free  of  relief  waves  aroused  bv  this  wave  at  the  starting  edge  for  a 
time  interval  *' 


D  2 
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Unless  ft"  >  i  /[l  -  t/" 2 ,7)" 2  ]  3,  end  edge  effect  cannot  i nte -  fere 
with  observations  of  the  wave  causing  it.  If  such  interference  begins 
at  a  distance  g  from  the  end  edge  at  time  tf  after  the  wave  first 
reaches  edge,  then 


«*  *  -•■*]* 


When  ft  <  U  / [ 1  -  U  2  D  2]  J,  the  end  edge  disturbance  begins  a  time 
[a  /V  ][l  -  U'2/D'2]  3  after  first  arrival  of  the  wave  at  the  edge,  and 
moves  from  the  end  edge  of  the  slab  along  the  interface  or  the  free  sur¬ 
face  toward  the  start  edge  at  speed  ft'[l  -  ft'2  L  2  +  ft'  2  D  2]~‘j  until  the 
d i s  t  ance 


has  been  covered,  when  the  disturbance  slows  A  point  on  the  interface 
or  free-surface  distant  ^  from  the  end  edge  will  be  free  of  end  edge 
effect  for  a  linu*  after  the  appearance  there  of  the  wave  causing  the 
effect: 


Kach  of  these  quantities  has  a  meaning  both  in  foam  slab  and  in  the  anvil, 
also  each  ma>  have  a  distinct  meaning  depending  or.  which  wave  is  consid¬ 
ered  We  will  leave  terms  referring  to  the  anvil  unprimed,  corresponding 
foam  quantities  will  be  singly  primed  If  the  wave  considered  to  be 
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causing  the  edge  effect  is  the  forerunner  or  the  main  wave,  terms  will 
carry  the  subscript  1  or  2,  respectively.  Immediately  we  can  write 

~  ^2  ~  =  U j  and  flj*  =  I/j*  .  If  we  are  to  make  observations  on  the 

main  wave  in  the  mirror  at  a  distance  Z  from  the  starting  edge,  the  mir¬ 
ror  surface  must  be  “clear"  of  forerunner  edge  effect  at  Z  until  the 
main  wave  reaches  Z,  and  the  foam-anvil  interface  must  be  clear  of  fore¬ 
runner  edge  effect  at 

Z'  *  Z  “  a  cot  6 

until  the  main  wave  reaches  Z' .  However,  in  all  practical  cases,  obser¬ 
vation  of  the  main  wave  will  be  limited  by  those  release  waves  in  the 
foam  behind  the  main  wave,  i.e.,  the  condition  for  stability  of  the  main 
wave  against  its  own  rarefactions  will  dominate  the  condition  for  stabil¬ 
ity  against  forerunner  relief  waves  in  the  foam.  Thus  we  require: 

Z'  >  b'2,  or  Z  >  fcj  +  a  cot  6, 

and 


a ' ( tan  &2  ~  tan  ) 


Both  and  U'2  are  average  speeds  for  the  particular  experimental  condi¬ 
tions  considered. 

As  an  example  we-  will  find  the  length  of  a  steel  anvil  for  use  with 
a  5-mm  slab  of  0.67  g/cc  polyurethane.  If  the  anvil  is  1  inch  thick, 
a  =  25.4  mm;  other  parameters  are  U  *  R  «  5.88  mm/fisec,  l/j  •  *  1 . 65  mm^sec . 

R'2  2  mm/^sec,  a'  *  5  mm,  and  U'2  1  mm/nsec.  If  the  barrier  or  explo¬ 

sive  is  in  contact  with  the  foam,  D  *  7.35  mm/fisec.  Thus  forerunner 
stability  is  reached  (a  cos  6  +  a'  cos  )  «  35.05  mm  *  1.38  inches  from 
the  starting  edge.  Equilibrium  of  the  main  wave  at  the  interface  comes 
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4> 


after  the  start,  or  a  distance 


\ 


b '  =  /?  2 2  "  a '  2J  =  12.3  mm  =  0.485  inch 


along  the  interface.  The  stable  B-front  will  then  first  appear  in  the 
anvil  mirror 


Z  =  b 2  +  a  cot  6  =  46.2  mm  =  1.82  inches 


from  the  starting  edge,  but  the  forerunner  stress  ahead  of  it  in  this 
region  of  the  anvil  will  have  been  weakened  by  starting  edge  effect. 
This  will  not  destroy  the  main  wave  data  but  can  be  avoided  by  using 
parts  of  the  mirror  yet  farther  from  the  starting  edge.  Such  a  region 
must  be  clear  for 


a '  (tan  ~  tan  ) 

D 


D 


if 

v  9 

-  1  4  -  ] 
tan  cos  —  _  tan  cos 


v: 


D 


i 

D  J 


2 . 0  ftsec 


after  the  passage  of  the  first  wave;  hence  it  will  be  found  at  distances 
Z  along  the  mirror  satisfying  the  inequality: 


Z  > 


rh 


in  which 


+  t 


o  i 


D 


Ti  +  (oi 


*)■  ■  -  (S  ■ ■  ■) 


Vi 


1  -  21 
D2 


01 


U 


(  u2' 

V  ~  D2, 


% 


.2  H sec 


=  {t 20i  U2  ~  a2}  1  ~  34  mm 
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Thus  Z  ^  120  mm  =  4.7  inches.  For  a  2- inch-thick  anvil  this  clear  area 
lies  beyond  a  distance 


Z  >  172  mm  =  6.8  inches 

from  the  starting  edge.  Since  reverberations  in  the  anvil  may  be  more 
annoying  than  starting  edge  effect  ahead  of  the  main  wave,  frequently 
this  lower  limit  on  anvil  length  has  been  ignored. 

End  edge  effect,  in  locked  foam  influences  the  main  wave  at  the  inter¬ 
face  in  the  region  within 


y, 

g‘2  =  [(^2**2  )  2  ~  a  2  ]  =  6  4  mrn  =  0.25  inch 

of  the  end  edge.  Since  a  cot  6  >  g'2  for  both  the  1-  and  2-inch  anvils, 
this  influence  does  not  reach  the  main-wave  jump  in  the  mirror.  In  the 
1-inch  anvil  the  fast  end  edge  effect  extends  a  distance 


h 


1 


aU  1 


1 - 


D2 


33 . 9  mm  =  1.33  inches 


from  the  end  edge;  and  =  67.8  mm  =  2.66  inches  in  the  2-inch  anvil. 
In  this  region  the  interval  between  a  wave  and  its  end  edge  reflection 
i  s 

2C 

—  =  0 . 272C  <  0.272 h. 

D  ~  1 


that  is,  9.22  Msec  when  a  =  1  inch  and  18.4  Msec  when  a  ■■  2  inches.  Out¬ 
side  this  region  the  interval  is  of  course  longer.  Since  the  main  wave 
follows  about  2.0  Msec  behind  the  first,  end  edge  effect  will  appear 
ahead  of  the  main  front  only  for 

£  £  D  =  7.35  mm  =  0.29  inch 

Since  an  ar^a  of  the  steel  mirror  near  the  end  edge  was  always 
masked,  fig.  2  show*-  no  end  edge  influence  on  the  first  two  wavefronts. 
The  anvil  used  in  the  experiment  of  which  Fig.  2  is  the  record  was  2  inche 
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thick  and  8%  inches  long  in  the  direction  of  the  detonation  run.  Under 
these  conditions  only  the  last  nine  or  ten  grid  images,  between  the  first 
and  second  fronts,  should  escape  rarefactions  stemming  from  the  start¬ 
ing  edge,  but  the  record  does  not  seem  to  show  any  effect  of  that  kind 
anywhere.  We  do  not  think  this  lateral  release  within  the  anvil  is  im¬ 
portant  in  any  of  our  experiments,  and  we  have  tentatively  interpreted 
anv  such  apparent  relaxation  of  stress  in  the  4-wave  that  we  have  seen 
as  arising  from  stress  relaxation  behind  the  forerunner  in  the  foam18  — 
due  either  to  an  internal  mechanism  in  the  elastically  strained  material,9 
or  to  sound  waves  moving  from  a  free  edge  at  speeds  higher  than  forerun- 
ne-  speed. 

For  reasons  not  understood,  4  - wave  apparent  speed  failed  to  coin¬ 
cide  within  the  uncertainty  limits  with  detonation  or  impact  speed  in 
three  experiments,  Shots  9,893  and  9,931  in  silica,  Shot  9,933  in  open¬ 
cell  aluminum.  In  all  three  the  first  wave  rose  extremely  slowly,  so 
perhaps  either  the  experimental  error  was  larger  than  expected  or  com¬ 
ponents  of  the  forerunner  were  subject  to  unknown  edge  effects.  These 
apparent  wave  speeds  were  lower  than  impact  speed;  a  wave  in  the  mirror 
surface  outrunning  starting  edge  effect  should  be  moving  faster  than  its 
stable  speed. 

If  for  1.3  g/cc  aluminum  we  estimate  that  \J'2  =  2.5  mm/gsec  and 
B'2  =6.0  mm//isec,  we  find  for  a  specimen  4  mm  thick  (Shots  9,930  and 
9,932)  that  the  main  front  will  be  stable  at  the  interface  a  distance 
46  mm  or  1.8  inches  from  the  starting  edge.  In  a  2-inch  anvil  this  in¬ 
stability  should  disappear  about  4.5  inches  from  the  starting  edge,  and 
thus  would  not  be  photographed  in  our  experiments.  We  cannot  explain  the 
seeming  contradiction  between  4-wave  behavior  in  Shots  9,930  and  9,932'- 
in  the  first  the  6-front  overtook  the  4-front,  yet  in  the  second,  when 
more  explosive  was  used,  a  stable  4-front  appeared  throughout  the  record. 

1  The  lighter  aluminum,  though,  may  be  heavily  affected  by  release 
waves  in  the  locked  foam.  Here  we  estimate  U'2  =  0.8  mm/fisec  and 
R'2  =  6.0  im/^sec.  In  a  typical  experiment  a1  =  6.0  mm,  hence 
b'2  =  493  mm  =  19.4  (inches;  that  is,  the  main  wave  may  never  outrun  re¬ 
lease  waves  behind  it  in  any  of  our  experiments  with  this  material.  Yet 
the  points  in  Fig.  4  stemming  from  these  experiments  do  not  seem  far  out 
of  agreement  with  points  from  heavier  material. 
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We  do  not  know  what  part  of  the  relaxation  behind  t  lie  B- front  seen 
in  the  smear  records  should  be  attributed  to  starting  edge  effect 

Because  of  the  toxicity  of  pulverized  beryllium  we  impacted  foamed 
specimens  of  this  metal  within  an  airtight  steel  Lank,  a  photograph  of 
which  appears  in  Fig.  13-  To  help  the  containment,  air  pressure  in  the 
tank  was  held  below  500  pi  during  the  explosion.  Both  beryllium  shots 
were  on  beveled  anvils,  but  the  foam  slabs  were  not  tapered  and  the  fly¬ 
ers  no",  tilted  because  the  expected  high  forerunner  speed  was  not  indi¬ 
cated  in  preliminary  static  testing.  The  experiments  were  mounted  on 
the  tank  cover  plate  as  sketched  in  Fig.  14. 

B.  QUASI -STATIC  MEASUREMENTS 

To  make  the  quasi-static  measurements  reported  in  Table  6,  we  hy- 
draulically  compressed  machine-cut  cylinders  of  foam  encased  in  close- 
fitting  steel  dies.  (Since  minimum  wall  thickness  was  3/8  inch,  die 
expansion  was  completely  negligible  in  our  tests.)  Displacement  and 
pressure- sensing  gages  connected  to  remote  writing  equipment  made  a 
simultaneous  plot  of  stress  in  the  hydraulic  chamber  vs.  engineering 
strain  in  the  cylinder. 

Without  exception  all  the  st ress - st ra i n  relations  in  quasi - st at i c 
(0.172  kbar/minute )  ,  1-D  compression  took  the  form  shown  for  polyurethane 
in  Fig.  15:  a  seemingly  straight  or  "elastic”  region  from  zero  to  the 
limit  of  proportionality,  followed  by  a  region  of  gradually  falling  sec¬ 
ond  derivative  of  stress  with  respect  to  strain.  There  is  a  point  of 
inflection  when  this  derivative  becomes  zero,  then  a  region  where  it 
rises  until  the  end  of  the  compression.  Most  of  our  tests  ended  at 
1.72-kbar  stress,  but  in  a  few  we  went  to  over  4. 14  kbar.  Except  in  the 
case  of  1.4  g/cc  ciosed-cell  aluminum,  density  at  1.72  kbar  is  in  the 
neighborhood  of  that  of  the  undistended  solid;  density  of  closed-cell 
aluminum  is  only  2.39  g/cc  at  4.14  kbar  (crystal  density  2.70  g/cc). 
Points  of  inflection  always  fall  below  1.72  kbar.  To  show  the  "locking” 
nature  of  the  compression,  relaxation  behavior  is  included  in  Fig.  15 
for  another  sample  of  the  same  kind  and  approximately  the  same  density 
of  foam.  Preheating  the  specimen  of  foam  generally  lowers  both  the 
clastic  modulus  and  yield,  as  shown  in  Fig.  16  for  polyurethane.  The 


* 

Baldwin  Press  Model  MA1B  opersted  st  6,000  psi/min.on  the  specimen. 
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FIG.  14(a)  EXPERIMENTS  MOUNTED  ON  COVER  PLATE 


0.67  g/cc  POLYURETHANE  (POLYTRON) 
2  diam,  2H  LENGTH 
LOAD  RATE  6,000  p*i/min 


FIG.  15  QUASI-STATIC  COMPRESSION  OF  POLYURETHANE 
INCLUDING  RELAXATION 


STRESS  -  Kbor 


FIG.  16  EFFECT  OF  PREHEATING  ON  COMPRESSION  OF  POLYURETHANE 
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difference  among  the  compression  curves  depending  on  the  manufacturer  of 
the  samples  in  Fig.  16  is  typical  of  this  material. 

Most  of  the  quasi-static  properties  we  found  in  our  foams  appear  in 
Table  6.  Predicted  elastic  wave  speeds  are  simply  the  square  roots  of 
the  ratio  of  modulus  to  starting  density.  Yield  stress  is  defined  as 
the  stress  at  which  the  departure  from  a  linear  relation  between  stress 
and  strain  amounts  to  0.2  percent  of  stress  or  strain.  Also  appearing 
in  the  table  are  compression,  the  slope  of  the  stress-strain  curve  at 
1.5-kbar  stress,  and  our  measurements  of  sound  speeds.  Uncertainty  limits 
refer  to  the  single  measurements  and  not  to  the  variation  that  is  likely 
among  different  specimens  of  the  same  kind  and  density  of  material.  From 
our  values  it  appears  that  yield  in  silica  is  raised  by  preheating  and 
is  little  affected  in  closed-cell  aluminum.  Otherwise  results  are  as 
expressed  in  Fig.  16. 

To  assess  the  proportion  of  apparent  stress, arising  in  friction  be¬ 
tween  the  specimen  foam  and  the  will  of  the  die,  we  did  two  things:  (1) 
compressed  various  samples  of  the  same  polyurethane  foam  but  of  different 
length-to-diameter  ratios  and  (2)  during  very  slow,  stepwise  compression 
of  aluminum  in  a  special  apparatus  held  lateral  expansion  to  zero  by  in¬ 
creasing  hydrostatic  pressure  on  the  sides  of  the  specimen.  Compression 
under  various  length-to-diameter  ratios  in  polyurethane  is  reported  in 
Fig.  17;  behavior  of  closed-cell  aluminum  under  1-D  strain  has  been 
shown  in  Fig.  18  by  points  lying  very  Cj.ose  to  the  elastic  portion  of  a 
stress-strain  curve  found  from  a  die-held  sample  of  the  same  material. 
Quasi-static  behavior  at  425°C  is  included  in  the  same  figure.  Differ¬ 
ences  between  curves  for  various  length-to-diameter  ratios  in  Fig.  17 
are  statistically  significant  but  are  seemingly  randomly  related  to  any 
possible  influence  of  die  friction,  which  should  increase  with  increasing 
ratio.  We  think  the  effect  of  variation  in  sample  properties  probably 
masks  the  friction  effect  in  Fig.  17. 

By  making  a  die  of  smaller  diameter  we  were  able  to  reach  higher 
stresses  during  quasi-static  testing  of  one  specimen  each  of  polyurethane 
and  closed-cell  aluminum.  Results  showed  the  usual  divergences  from  tests 
to  lower  stress  limits  due  to  differences  in  the  physical  properties  of 
samples;  but  final  strain  reached  in  aluminum  was  considerably  higher 
than  in  earlier  tests.  Maximum  stress  was  4.35  kbar  at  which  strain  in 
1.432  g/cc  aluminum  was  40.5  percent  (89.2  percent  of  crystal  density) 
and  in  0.698  g/cc  CPB  polyurethane.  48.  2  percent  (final  density  1.34  g/er). 
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0  66  Q/CC  POLYURETHANE  (POLYTRON) 
LOAD  RATE,  6,000  p»i/m.n 


FIG.  17  COMPRESSION  OF  POLYURETHANE  UMDER  VARIOUS 
LENGTH-TO-DIAMETER  RATIOS 
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CLOSED-CELL  ALUMINUM 
(MD-AK) 


STRAIN 


FIG.  18  EFFECT  OF  PREHEATING  ON  COMPRESSION  AND 

ONE-DIMENSIONAL  STRAIN  OF  CLOSED-CELL  ALUMINUM 


Figures  1  4  and  2 0 ,  showing  the  temperature  effect  on  compressive 
behavior  in  open-cell  aluminum  and  graphite,  respect icclv,  are  similar 
to  figs.  16  and  18  The  curves  above  1.0-kbar  stress  show  that  thete  is 
considerable  structural  rigidity  loft  in  the  foam  even  though  the  propor¬ 
tional  limit  has  been  greatly  exceeded,  and  that  this  rigidity  is  reduced 
both  below  and  above  the  limit  by  heat.  The  pre s su re  - vo 1 ume  equation  of 
state  of  >.n  ordinary  solid  or  liquid  generally  shows,  of  course,  an  op¬ 
posite  pressure  dependence  on  temperature  than  that  seen  in  these  figures 
Figure  21.  taken  from  compression  data  in  silica,  is  then  remarkable 
above  the  static  yield  point  in  showing  the  direction  of  temperature- 
pressure  relation  expected  of  an  ordinary  material.  Since  the  apparent 
“Griineisen  constar*  '  irrespond i ng  to  the  temperature  dependence  of 

pressure  in  t\-s  foam  is  of  the  order  of  0.2.  i.e.,  approximately  one- 

tenth  that  ol  ordinary  sc  1 ’.ds  and  liquids,  there  must  be  some  structural 
weakening  due  to  heating.  (Pure  silica  fuses  above  1,000°C.  )  Except  for 
temperature  of  sample  and  die  and  slight  difference  in  sample  densities, 
che  conditions  of  the  tests  leading  to  the  two  curves  in  Fig.  23  were 
the  same.  The  samples  were  taken  from  neighboring  locations  in  the  same 
1 u rge  block 

C.  EXPERIMENTS  IN  PREHEATED  FOAM 

To  explore  the  effect  of  preheating  on  the  shock  response  of  foams, 
the  usual  2-D  experiments  were  enclosed  within  asbestos  insulating  boxes 
fitted  with  electrical  heating  wire.  Power  was  slowly  applied  until 
thermocouples  on  the  free  surfaces  of  both  foam  and  anvil  read  the  same 
temperature.  Figure  22  shows  the  arrangement  of  this  apparatus.  During 
the  heating,  the  flyer  and  explosive  drawn  in  Fig  22,  of  course,  were 
replaced  by  an  insulating  cover;  when  a  temperature  slightly  above  the 
wanted  level  had  been  reached  the  flyer  and  explosive  were  put  in  place 
along  with  a  cold,  removable  heat  shield.  With  this  precaution  we  never 
found  the  temperature  in  the  explosive  just  before  initiation  above  60  C 
However,  because  the  explosive  was  hotter  than  usual,  its  detonation 
speed  was  measured  in  every  experiment  involving  preheating.  (Results 
are  entered  in  Table  1  in  the  Apparent  Wave  Speed  column.)  To  reach 
foam  temperatures  near  300°C.  we  heated  for  approximately  1  hour  with 
800  watts. 

Preliminary  static  work  had  shown  which  temperatures  should  have  a 
marked  effect  of  foam  behavior  without  melting  or  decomposing  the  material 
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FIG.  19  EFFECT  OF  PREHEATING  ON  COMPRESSION  OF 
OPEN-CELL  ALUMINUM 
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FIG.  20  EFFECT  OF  PREHEATING  ON  COMPRESSION  OF 
ATJ  GRAPHITE 
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FIG.  21  EFFECT  OF  PREHEATING  ON  COMPRESSION  OF  SILICA 
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FIG.  22  VIEW  OF  EXPERIMENT  WITH  THERMAL  ISOLATION 


or  wholly  destroying  its  elasticity;  as  a  result  we  sought  to  preheat 
0.67  g/cc  polyurethane  to  approximately  120°C,  and  1.4  g/cc  aluminum  to 
400°C.  Several  technical  problems  arose  in  reaching  the  higher  tempera¬ 
ture  level:  to  find  a  glue  that  would  hold  the  curling  and  cracking 
foam  slabs  on  to  the  anvil  at  high  temperature  took  some  time;  to  avoid 
deterioration  of  the  quality  of  the  anvil  mirror  it  was  vacuum -coated 
with  aluminum.  (Gold  did  not  withstand  heating  so  well.)  Also  the 
smear  camera  record  from  Shot  10,360  (318°C)  is  blurred  by  convection 
currents  in  front  of  the  anvil  mirror;  any  more  experiments  at  that  tem¬ 
perature  level  should  be  done  with  some  sort  of  transparent  heat  isola¬ 
tion  between  air  in  the  heating  box  and  air  outside. 

The  preheated  foam  shots  were  Nos.  10,331  and  10,332  in  polyurethane  (CPH) 
and  10,360  in  aluminum.  Pertinent  data  and  results  are  given  in  Table  1. 

Polyurethane  at  about  lli°C  transmit  a  first  wave  reduced  in  strength 
by  about  40  percent,  although  the  static  yield  at  120°C  is  lessened  by 
70  percent.  Closed-cell  aluminum,  1.4  g/cc,  at  318°C  shows  a  forerunner 
stress  in  the  anvil  lowered  more  han  50  percent,  while  the  static  propor¬ 
tional  limit  falls  only  30  percent  under  heating  to  426°C.  Measurements 
of  average  first-wave  speeds  in  preheated  material  were  inconclusive. 

D.  AFTER- SHOT  OBSERVATIONS 

Even  though  the  smear  record  has  a  clear  B- wave,  specimens  of  the 
polyurethane  foam  recovered  from  the  experiment  after  impact  show  less 
than  5  percent  residual  strain  due  to  the  explosively  driven  compression. 
When  the  barrier  or  flyer  is  used,  there  is  even  no  scoring  of  the  struck 
surface.  Beryllium  also  shows  little  outward  change  or  residual  strain 
after  shocking.  Aluminum  foam  is  found  compressed  to  near  its  unfonmed 
density;  silica  and  graphite  (PT0114)  have  not  been  recovered.  Samples 
machine-cut  from  recovered  specimens  of  polyurethane  and  aluminum  foams 
show  the  densities  and  sound  speeds  listed  in  Table  7.  All  specimens 
in  the  experiments  mentioned  were  subjected  to  a  main  or  8- wave  except 
perhaps  one,  that  from  Shot  9,892.  In  the  smear  record  from  Shot  9,892 
there  is  a  pressure  rise  of  considerable  strength  coming  into  the  anvil 
between  the  first  and  second  reverberations  of  the  4-front  in  the  anvil; 
this  second  wave  is  probably  a  locking  wave.  Aluminum  is  recovered  with 

Epoiylitt,  C4B644,  trada  dim  of  Epoiylita  Corp.,  South  El  Nonta,  Californio. 
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Table  7 

AFTERSHOT  DENSITIES  AND  SOUND  SPEEDS  OF  FOAMS 


SHOT 

NO. 

MATERIAL 

STARTING 
DENSITY 
(g/cc ) 

PEAK 

PRESSURE 

(kbar) 

AFTERSHOT 

DENSITY 

(g/cc, 

±  1 %) 

AFTERSHOT 

SOUND 

SPEED 

(naZ/rsec, 

±  1%) 

TYPICAL 
BEFORE 
SHOT  SOUND 
SPEED 
(»m//isec  ) 

9808 

P-C  (Polytron) 

^0.67 

1.1  -  4.3 

0.672 

9832 

P-C  (Polytron) 

0.645 

5.0 

0.673 

3  | 

Average  0.673 

1.30 

1.73 

10260 

P-C  (CPR) 

~0.714 

0.50 

0.716 

1.83 

Average  0.716 

1.83 

2.11 

9834 

MD-AO 

0.925 

0.81 

1.94 

3.52 

9835 

MD-AO 

0.927 

not  measured 

3.02 

3.62 

9810 

MD-AO 

0.909 

not  measured 

2.16 

3.44 

9803 

MD-AO 

0.960 

1.6 

2.06 

3.03 

9802 

MD-AO 

0.939 

1.9 

2.16 

3.37 

9801 

MD-AO 

0.95 

5.0 

2.30 

3.20 

9800 

MD-AO 

0.96 

12.0 

2.06 

-- 

Avera 

ge  2.  10 

3.20 

2.54 

0930 

MD-AO 

1.38 

1.4 

1.98 

3.96 

9933 

MD-AO 

1.28 

8.1 

2.32 

2.36 

Averi 

ge  2.15 

3. 16 

4.75 

9892 

MD-AK 

l.  445 

1.49 

1.90 

3.48 

9894 

MD-AK 

1.46 

4.5 

2.19 

2.49 

Average  2.04 

3.00 

4.60 

9765 

MD-AO 

~0.76 

not  measured 

2.21 

-  - 

9792 

MD-AO 

0.76 

20.0 

2.04 

-• 

9809 

MD-AO 

0.688 

6.6 

2. 10 

-- 

Average  2. 17 

-- 

10082 

C-ATJ 

1.73 

2.6 

1.76 

1.99 

10084 

C-ATJ 

1.73 

5.7 

1.78 

2.00 

Averi 

ge  1.77 

2.00 

2.32 
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densities  between  2.0  and  2.2  g/cc  or  between  72  and  82  percent  of  the 
unfoamed  density;  during  the  experiments,  of  course,  higher  densities 
may  be  reached. 

After-shot  sound  speed  (1.3  mm/Vsec)  in  polyurethane  is  markedly 
lower  than  sound  speed  in  unshocked  material  (1.7  mm/^sec),  but  aluminum 
is  not  consistent  in  this  regard.  Shocking  seems  to  raise  the  final 
sound  speed  in  0.9  g/cc  aluminum  but  to  lower  it  in  1.4  g/cc  aluminum. 

.  Thin  sections  made  from  both  impacted  and  unused  polyurethane  and 
studied  und®r  a  petrographic  microscope  show  no  differences  except  that 
the  section  from  the  struck  specimen  appears  to  be  bi re f r ingent .  This 
suggests  the  possibility  of  the  existence  of  residual  strain. 

From  flash  X-ray  shadowgrams  we  know  that  0.67  g/cc  polyurethane 
does  compress  40  to  50  percent  under  explosive  impact.  These  pictures 
show  also  that  the  compression  lasts  at  least  17  ^ sec  after  detonation. 

E.  VALUE  OF  EXPLOSIVE  MOMENTUM 


In  those  experiments  described  in  Table  1  as  having  standoff,  that 
is,  a  space  between  the  parallel  flyer  and  foam  surfaces,  we  can  later 
assign  an  accurate  value  to  the  starting  momentum  in  the  impact  on  the 
foam  by  subsidiary  measurements  of  the  angle  between  the  moving  and  still 
portion  of  the  flyer  during  like  detonations.  Without  this  measurement 
momentum  can  be  related  to  weights  per  unit  area  of  explosive  and  accel¬ 
erated  plate  by  the  Gurney 19  the  ory  *  and  G.  R.  Abrahamson's  observations20 

•  i 

See  Fig.  10,  Foulea  and  Curran.  Since  damnation  apeed  »ea  about  7.3  M/^aac  and  tha  apeciaen  *aa 
i  inchea  long,  tha  tiae  interval  bataaan  tha  aarty  and  lata  atagaa  of  coapraaaion  la  17  »iatc.  Micro- 
photodenaitoaetar  atudiaa  of  tha  negative  of  thtu  abadoagraai  indicate  a  coapraaaion  bahind  tha  locking 
•ava  in  tha  ranga  aontionad  above.  Thin  ia  tha  ragion  of  final  atrain  aaan  in  atatic  coapraaaion  to  v  ikbar. 

^  Thia  ralataa  final  plata  apoad  V  to  tha  auparficial  charge  or  eaploaive  danaity  "c"  and  tha  auparficial 
plata  denait<  V*  aa  folloara: 


0.6C£ 


1  H 


0.8  ♦  j  ♦  0.2 


WJ 


In  oar  »ork  Fa  ■  / q.  input  aoaantaa  danaity  to  tha  foaa  and  »a  can  rewrite  tha  above  aapraeaioa  aa 


•hich  often  caa  be  approaiaated  by 


’  j'  I, 

l  ♦  O.h  I  ♦  $.33(  4-  ) 

i.  .  _ _ —A  ,  a  J—l 

•[’(^)  ■] 
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2  /„ 
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of  the  Gurney  constant,  GQ;  but  since  this  constant  strongly  depends  upon 
the  nature  of  the  material  accelerated  by  the  explosive,  the  amount  of 
explosive  momentum  delivered  to  a  foam  lying  against  the  aluminum  barrier 
at  the  start  of  detonation  is  uncertain  within  limits. 

For  instance,  a  ma te r i a  1  o f dens  i  t y  0.  67  g/cc  impacted  by  EL-506D  exp  1 o- 
sive,  Abrahamson  reports,  has  a  Gurney  constant,  G0,  0f  1.88  ±  0.19  x  10s  cm  sec; 
for  aluminum  (»q  =  2.31  i  0.23  x  10^cm/sec.  To  put  1  x  10*  taps  of  impulse 
into  an  aluminum  flyer  0.012  inch  thick  as  in  our  2-0  experiments  with 
standoff,  Gurney’s  formula  requires  a  layer  of  EL-506D  explosive  0.021  inch 
thick  (density  1.4  g/cc).  For  the  experiments  without  standoff  we  can  use 
a  weighted  average  value  ol  G0,  i.e., 

c,  =  - L_ 

0  "ai  + 

where  m  refers  to  mass  per  unit  area,  subscript  Al  means  aluminum,  and 
subscript  /,  foam.  Thus  when  a  0 . 01 2- i nch - th i ck  flyer  is  used  against 
5  mm  of  polyurethane  foam  of  0,67  g  cc  density,  C,  1  97  and  the  explo 
sive  thickness  to  give  1  *  104  dyne  sec  cm  2  is  calculated  as  0.0180  inch. 
When  the  flyer  alone  is  impelled,  the  thickness  of  explosive  for  the  same 
momentum  density  in  the  same  geometry  is  0.022  inch 

Values  of  the  abscissa  in  Figs  4  and  5  are  usual  lv  uncertain  on  two 
scores:  the  foregoing  theory  does  not  give  unchallengeable  values  of 

input  momentum  density  in  2-0  experiments,  and  the  foam  thickness  in 
the  (rhanging)  direction  normal  to  the  shock  front  is  unknown  In  plot 
ting  points  in  those  two  figures  from  the  2*0  experiments.  w»  have  used 
one  constant  ratio  between  explosive  thickness  and  /„  for  experiments 
without  standoff  and  another  in  experiments  with  standoff,  and  we  have 
overlooked  the  greater  effective  foam  thickness  due  to  obliquitv  of  i he 
fronts;  the  horizontal  bars  through  the  points  in  the  figures  are  est  i 
mates  of  our  uncertainties  in  these  respects. 

F  FOAM  UNIFORMITY 

The  rather  wide  scatter  of  the  reduced  data  shows,  we  think  (hr 
importance  of  variation  in  physical  properties  among  even  neighboring 
sites  in  one  foam  block  In  the  absence  of  evidence  «r  see  no  reason  to 
regard  density  as  the  only  controlling  physical  characteristic  but  densitv 
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variations  are  widely  found  in  foam  specimens.  We  have  made  both  \-rav 
and  light  shadowgrams  through  foam  layers,  and  we  have  measured  average 
densities  from  small  samples  cut  from  neighboring  sites.  For  example, 
all  specimens  of  the  same  kind  and  of  approximate  densities  listed  in 
Table  1  were  supplied  in  single  pieces,  there  is  no  reason  to  believe 
that  the  variation  in  average  densities  seen  there  does  not  extend  to 
each  single  specimen  as  well.  But  if  we  traverse  with  a  m 1 c rophot odens i - 
tometer  the  X-ray  shadowgram,  we  should  see  important  trends  Figure  23 
'resulted  from  two  traverses  along  different  axes  of  a  single  shadowgram 
through  a  1-inch-thick  slab  of  0.9  g/cc  aluminum.  The  density  variation 


distance  —  inehtt  ***.t.} 

FIG.  23  MICROPHOTODENSITOMETER  SURVEY  OF  X-RAY  SHADOWGRAMS  OF 
A  THIN  FOAM  LAYER 

fro*  place  to  place  way  be  as  murh  as  20  percent  XII  aluminum  foam 
showed  obvious  outward  marks  of  irregularity  in  density.  especially  the 
lighter  density  foam  Other  foams  showed  less  variation 

We  probed  our  foamed  beryllium  .specimens  u  1 1  nsnn  i  r  a  I  1  y  .  and  here 
•gam  considerable  itregularitv  was  found  in  s"und  speed  within  one  sample 
Of  four  measurements  at  different  sites  in  one  plate  >  *  8'j  ‘  0  23b  . 
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the  average  was  6.71  mm/^sec;  the  standard  deviation,  0.40  mm  i-ser.  and 
the  greatest  deviation  from  the  mean,  0.54  mm  ,xse<  .  Of  four  'vlsnders 
1  1/8  inch  diameter,  1  inch  long,  the  ratio  of  average  density  to  sound 
speed  through  the  length  of  the  sample  was  1.52  *  10  t  O.ligser  cm 
the  greatest  deviation,  0.  16  x  10  6  g  sec/  on  *,  hence  the  correlation  be¬ 
tween  average  foam  density  and  sound  speed  was  not  better  than  10  percent. 


4.  THEORY 


A.  GENERAL 

la  the  three  models  used  in  calculations  so  far  we  have  assumed  a 
uniform  material  into  which  the  conventional  1-D  shock  discontinuities 
run  according  to  the  conservation  laws.  In  all  we  have  taken  an  equation 
of  state  (in  effect,  a  combined  Hugoniot  and  isentropic  equation  of  state) 
as  a  simple  relation  between  pressure  and  specific  volume,  and,  through  a 
linear  “elastic”  region  of  small  compression  in  this  relation,  in  all 
have  achieved  breakup  of  a  single  strong  shock  into:  (1)  the  fast  forerunner 
behind  which  strain  is  small  and  (2)  the  main  shock  behind  which  strain  is 
great.  Entropy  change,  always  concentrated  in  the  main  front,  is  not 
important  in  calculation  because  all  isentropes  containing  any  state 
behind  the  main  front  are  assumed  to  coincide.  In  the  simple  model, 
completely  set  forth  by  Fowlcs  and  Curran1  and  Rempel,*  the  main  front 
is  like  a  snowball,  gathering  more  and  mere  condensed  foam  behind  itself 
and  moving  slower  and  slower,  with  ail  parts  of  the  gathered  mass  moving 
at  the  same  speed  and  having  the  same  density.  Both  the  other  models 
are  slightly  more  sophisticated — but  not  more  successful  —  in  allowing 
waves  of  small  changes  to  move  through  the  condensed  mass. 

In  the  following  two  sections  we  present  the  basic  assumptions  and 
some  results  of  1*D  flow  calculations  by  the  method  of  characteristics 
and  by  the  method  of  artificial  viscosity.  The  first  yielded  pressure 
and  particle  speed  distributions  2.5  //sec  after  impact  in  a  foam  half¬ 
space  struck  by  a  thin  metal  flyer.  The  same  distributions  were  derived 
by  the  second  method.  Since  the  results  compared  closely  and  since  a 
code  for  machine  computation  by  artificial  viscosity  was  available,*  the 
method  of  artificial  viscosity  was  used  in  an  attempt  to  predict  the 
pressure  history  in  the  anvil  during  one  of  our  1-D  experiments. 

In  the  characteristics  method  the  Rank  ine-Hugon  iot  jump  conditions 
are  satisfied  at  the  locking  shock  front  behind  which  release  waves 


* 

See  J.  F.rkmnn,  Phase  2  of  this  report. 
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reverberate.  Interface  states  are  found  by  standard  impedance  matching 
techniques  in  the  pressure-particle  speed  (P  -  u)  plane.  Locations  of 
fronts  and  interfaces  are  followed  in  the  di stance- time  (x  -  t)  plane. 

In  the  problem  considered  the  forerunner  appears  almost  immediately  upon 
impact. 


B.  CALCULATION  BY  THE  METHOD  OF  CHARACTERISTICS 

Figures  8  and  9  result  from  a  graphical  solution  of  the  flow  equations 
carried  out  in  the  pressure-particle  speed  (P  -  u)  and  di stance- time  ( x  -  t) 
planes  with  the  help  of  simplified  equations  of  state.  We  tried  to  simulate 
conditions  in  a  1  ~  D  experiment,  that  is,  a  solid  aluminum  flyer  plate 
0.020  inch  thick  struck  a  polyurethane  foam  half-space  at  a  speed  of 
0.08  mm//i.sec.  The  foam's  original  specific  volume  FQ  was  1.515  cc/g,  locked 
volume  Fj  at  the  elastic  limit  Pe  was  1.00  cc/g,  P(  was  1.00  kbar,  and 
elastic  speed  V  was  1.65  mm/^sec.  The  locked  foam  and  the  solid  aluminum 
had  sound  speeds,  c,  at  zero  pressure  of  2.12*  and  6.23  mm//isec,  respectively, 
and  both  were  given  equations  of  state  which  were  single  straight  lines  in 
the  P  ~  V  plane, *  implying  that  rarefactions  resulting  from  reflections  of 
shocks  at  free  surfaces  were  discontinuities. 


Thus  if  m  is  the  slope  of  the  equation  of  state,  i.e., 


,2  =  v-2 


1  V'  -  V 
*1  M 


and 


m  = 


V'  -  V 
yl  *1 


(u  is  the  particle  speed,  P,  the  pressure,  and  F,  the  volume),  then  waves 
behind  the  locking  front  move  at  speeds 


\U\  =  ±u  +  Vnfi 


y; 

±U  +  C  -  -  P 

c 


Slopes  of  characteristics  for  locked  foam  and  aluminum  in  the  P  “  u  plane 
are  independent  of  the  pressure: 

Chosen  to  (free  with  esrlier  work  using  the  method  of  ertificisl  viscosity  in  which  locked  fesa  wss 
trested  by  sn  spproxiaste  equstion  of  stete  for  lucite. 

■  This  equstion  of  state  for  ths  locked  foem  is  shown  dsshed  in  Fig.  7. 
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A P 

Au 


\l’  ±  u 


Pressure  and  particle  speed  in  the  foam  behind  the  forerunner  are  related 
bv 

P  =  —  u  if  P  <  P 

Vo  "  ' 

In  the  absence  of  a  forerunner  the  relation  between  pressure  and  particle 
speed  just  behind  the  locking  shock  in  the  foam  is: 


2  V  >  2 
r  l 

if  P  >  P  where 

—  C 


When  the  locking  shock  in  the  foam  is  running  into  the  elastically  strained 
region,  the  jump  conditions  and  the  equation  of  state  imply  immediately 
behind  the  shock  front: 


if  Pf  <  P  <  Pe. 

Figure  24  shows  the  relation  written  above  between  P  and  a  and  a  few 
of  the  early  characteristics  used  in  the  solution  of  the  problem  stated. 
This  solution  has  not  been  taken  beyond  2.5  //sec  of  problem  lime. 


U.  CALCULATION  BY  THE  METHOD  OF  AHTIFIC! \L  VISCOSITY 

In  Figs.  8  and  Q  (except  in  the  forerunner  front  t  the  flow  snlul 
bv  the  method  of  artificial  viscosity  ( Q-  met  hods  )*  '*  agree  iloselv  with 

*  The  coding  for  thi»  cnlcuUtion  follnm  in  •ml  particular*  that  found  in  H*r»ann  o  at.1' 
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FIG.  24  RELATION  BETWEEN  PRESSURE  AND  PARTICLE  SPEED  JUST  BEHIND 
THE  LOCKED  FRONT  IN  POLYURETHANE  FOAM  HALF-SPACE  STRUCK 
BY  ALUMINUM  FLYER 

those  found  by  the  method  of  characteristics  outlined  in  the  foregoing 
section,  although  because  of  the  short  time  covered  the  test  is  not  a 
severe  one.  We  think,  however,  our  use  of  the  Q-method  to  explore  rapidly 
the  effects  upon  pressure  history  of  certain  changes  in  the  equation  of 
state  is  valid.  The  only  important  difference  in  the  forecasts  of  pressure 
history  between  the  two  calculations  lies  in  the  slow  pressure  rise  in  the 
forerunner  front  expected  from  Q-method.  The  artificial  parameter  0  is 
also  high  in  value  in  that  region.  We  look  upon  it  as  a  shortcoming  of 
the  (^-method  that  results  in  regions  of  high  Q  are  wrong  to  some  unknown 
extent.  Q  is  defined  as 

Q  •  ( |  At*  i  4  4c*)lAu|  ~ 

where  q  and  4  are  free  numeric  constants,  Au  is  particle  speed  difference 
between  adjacent  cells  of  constant  mass,  c  is  local  sound  speed,  and  V  is 
local  specific  volume. 
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As  seen  in  Fig.  7,  although  the  equation  of  state  used  for  locked 
polyurethane  in  the  Q-method  is  more  complicated  than  the  simple  relation 
used  in  the  characteristics  method,  the  differences  in  the  pressure  range 
of  interest  are  negligible.* 

Although  the  distributions  shown  in  Figs.  8  and  9  are  uninfluenced 
by  it,  a  rigid  wall  at  x  =  5.0  mm  bounded  the  foam  during  the  calcula¬ 
tions  by  the  Q-method,  and  Fig.  6  shows  toe  pressure  history  predicted  in 
the  wall  by  the  same  computation  which  gave  rise  in  part  to  Figs.  8  and  9. 
There  are  a  number  of  obvious  points  of  strong  disagreement  between  the 
predicted  and  observed  histories  shown  in  Fig.  6.  We  would  expect  only 
a  10  to  15  percent  lowering  of  peak  stress  in  the  predictior  if  we  gave 
to  the  wall  the  shock  impedance  of  steel,  we  have  sought  to  lower  this 
peak  further  by  changing  the  equation  of  state  assumed  for  the  locked 
foam  in  such  a  way  as  to  lower  its  shock  impedance,  but  this  leads  to  a 
longer  time  interval  between  B-  and  C-fronts  (main  and  reverberation) 
than  observed.  The  rise  time  in  the  fl-front  is  increased  slightly  in  two 
different  ways:  (1)  by  moving  the  point  (  F  ,  P  )  down  along  the  elastic 
line  halfway  to  (V0,PQ)  and  connecting  ( l'# , Pf  )  with  ( lj , P, )  by  a  straight 
line,  and  (2)  by  moving  (Vf,Pt)  along  the  elastic  line  to  0.85  kbar  and 
connecting  to  the  existing  equation  for  locked  foam  at  P  -  2  kbar  through 
two  second- degree  arcs,  meeting  at  (P  -  1  kbar,  T  =  ( 1 #  +  l()  2l  where 
derivatives  dP/HV  on  both  arcs  vanish.  Figure  25  illustrates  these  modi¬ 
fications.  In  both  (1)  and  (2).  the  pressure  histors  it;  the  rigid  wall 
begins  a  slow  rise  after  an  interval  of  constant  pressure,  but  t  lie  arrival 
of  the  main  front  is  still  quite  distinct,  the  B-C  interval  is  2  gser  too 
long  in  case  (2)  but  approximately  as  seen  experimentally  in  case  tl). 

We  can  successfully  duplicate  the  failure  of  the  anvil  pressure  to  return 
to  *ero  within  8  i*sec  after  .4- front  arrival  bv  •  boost  tig  a  relaxation  path 
for  locked  foam  along  the  compression  path,  as  for  example  in  <  asr  (2*. 
None  of  the  many  P  “  I  constitutive  relations  tried  so  far  begins  to  re¬ 
move  the  clear  pressure  relaxation  between  the  B  and  f.'-praks.  although 
th?  two  peaks  can  easily  be  made  less  distinct  but  farther  apart  in  lime. 


*  tl  »k«»l4  k*  »m*4  «!>•»  »»*<  *»  lk*  l>  MtM  lk*  r»pr*«*»t»<i  h»  T  o»«  t  k*  •vppl*. 

vtlH  b»  Ik*  intlniiia  of  tl itu  tlnftf  Of  In**  P  P  kM»**»  tk*  i  lf.  *M  iVj.I’j  >  -•'•'I 

Ik*  pro*i*to«  I»tk*  (rna  lk^»*  ml*r«*di*tr  •!•«*•  fr  l-l«  |>i*»»wr*.  »k|«*i  •»  k**»  t»k»o  f»r*ll*l 
«o  tk*  r**k  (ro*  O'  Ifj.fjl.  So**  of  lk*»*  >«  p.«»»ibl»  10  •  ppl  ir«l  ««o  I  Sr  ■*|K.«| 

of  rk»r*rl*rt*l it*. 
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FIG.  25  MODIFICATIONS  OF  EQUATION  OF  STATE  FOR 
Q-CALCULATIONS 


D.  INTERPRETATION  OF  EXPERIMENTAL  RESULTS  BY  THE 
SIMPLE  THEORY 

From  calculated  impact  speed  and  from  the  peak  stresses  in  the  main 
waves  seen  in  three  experiments  with  1  ~  D  symmetry,21  we  have  tried  deducing 
consistent  values  for  average  shock  speed  in  the  locked  foam.  If  our  inter¬ 
pretation  of  the  experiments  has  been  correct  in  locked  0.67  g/cc  polyure¬ 
thane,  the  calculations  show  an  average  shock  speed  increase  from  0.8  to 
1.5  mm//Li.sec  with  stress  in  the  range  2  to  4  kbar.  The  results  of  the  cal¬ 
culation  can  be  seen  in  Table  8.  The  computation  of  those  results  is  ex¬ 
plained  in  the  following. 

Using  the  simple  theory  of  e 1  as t i c- r i gi d  locking  behavior  (Fowlesand 
Curran,1  and  Rempei2),  the  measured  flyer  momentum  of  the  1  -  D  experiments 
in  polyurethane,  and  the  accurately  measured  forerunner  speed  in  polyure¬ 
thane,  we  explore  the  effect  of  assuming  different  values  of  locked  density 
until  we  find  a  range  of  values  in  each  experiment  that  leads  to  calculated 
time  intervals  between  arrivals  of  first  and  second  waves  (A  ~  B  intervals) 
falling  within  the  range  actually  seen  in  the  same  1  -  D  experiments.  Thus 
in  column  6  of  Table  8  appear  the  measured  least  and  greatest  values  of  the 
A  -  B  intervals,  and  in  column,  the  calculated  intervals  corresponding  to 
each  of  the  values  of  specific  volume  in  column  8.  We  see  that  a  locked 
specific  volume  between  0.90  and  1.05  cc/g  can  always  be  found  to  produce 
a  calculated  A  ~  B  interval  to  match  any  observed  value.  In  fact  1.00  cc/g 
is  a  good  compromise;  so  we  then  calculate  particle  speed  and  stress  behind 
the  locking  wave  at  the  moment  of  collision  of  the  two  locking  waves  by  as¬ 
suming  the  locked  volume  is  1.00  cc/g.  Next,  assuming  the  shock  speed  and 
density  stay  constant  within  the  locked  mass  after  .this  collision  and  as¬ 
suming  shock  impedance  of  steel  is  4.62  x  106  g  cm-2  sec-1,  we  compute  the 
value  of  foam  shock  impedance  corresponding  to  the  observed  peak  pressure 
of  the  main  wave  (column  7).  From  this  we  get  the  shock  speed  in  column  11. 
(Sound  speeds  in  two  samples  of  this  foam  after  being  subjected  to  locking 
shocks  were  1.32  and  1.29  min/gsec,  Shots  9808  and  9832,  respectively.  As 
noted  elsewhere  [Table  7]  the  recovered  samples  appear  to  hnve  their  original 
average  density.!  Our  experiments  give  us  a  verification  of  these  inferred 
speeds  in  locked  foam  because  a  third  front  is  reflected  back  to  the  anvil 
from  the  sharp  impedance  mismatch  at  the  foam- flyer  interface. 

Entries  in  the  last  column  of  Table  8  show  thut  the  average  reverber¬ 
ation  speed  based  on  B  ~  C  intervals  is  much  higher  than  speeds  inferred 
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from  the  computed  speed  of  locked  mass  and  the  pressure  induced  in  the 
anvil  by  the  main  wave.  The  disagreement  is  so  great  we  are  led  to  doubt 
whether  the  third  wave  originates  in  the  meeting  of  main  wave  and  anvil, 
but  wonder  if  it  stems  instead  from  a  significantly  earlier  collision  of 
a  wal 1- reflected  wave  and  the  main  wave.  If  this  is  true  there  should  be 
a  fourth  wave  in  our  record  of  anvil  motion,  and  there  may  have  been  in 
Shot  9217;  the  observation  was  broken  off  too  soon.  However,  it  is  not 
likely  we  would  have  missed  the  fourth  wave  in  Shot  0175.  Kurt  lie  rmo  re 
i f  we  argue  that  the  fourth  wave  is  not  seen  because  the  ft  am  has  bounced 
off  the  anvil  before  the  wave  crosses  the  foam,  we  must  dispute  the  ac¬ 
curacy  or  the  pertinence  of  the  record  of  Shot  0155  which  shows  pressure 
held  on  the  anvil  for  a  long  time  after  the  third  wave  enters  it  (Fig.  p )  , 
(It  is  theoretically  possible  for  the  locked  foam  to  bounce  off  the  steel 
wall  before  the  fourth  wave  reaches  it,  provided  there  is  substantial  in¬ 
crease  in  the  shock  impedance  o f  locked  material  with  higher  pressure.  ) 

Another  way  to  soften  the  contrast  between  columns  11  and  11  of 
Table  8  is  to  assume  that  the  pressure  behind  the  main  wave  when  it 
collides  with  tin-  wall  or  the  wall  reflection  of  the  forerunner  is  no 
greater  than  forerunner  stress  P f.  I  rider  this  assumption  the  shock  im¬ 
pedance  of  the  locked  foam  must  he  greater  to  give  rise  to  the  observed 
peak  /{-pressure;  but  the  resulting  increase  in  values  of  f  in  column  II 
for  Shots  9153  and  9217  amounts  to  only  about  13  percent.  Kven  taking 
pressure  behind  the  B  front  as  zero  does  not  wholly  reconcile  columns  11 
and  11. 


K.  FOAM  KFFHCTl YKNKSS  AS  A  UNCTION  OF  PAR  AMKTF.RS 
FROM  Till-  SIMPI.F.  TIIKORY 

Mien  the  main  wave  front  strikes  the  wall  reflection  of  the  fore¬ 
runner  just  as  the  front  disappears,  according  to  the  simple  theory,  t  lie- 
peak  pressure  in  the  anvil  is 


P  -  P, 


(  2  > 


where  Pf  ”  stress  in  forerunner  and  l\  is  the  shock  speed  in  t  he-  locked 
foam.  If  the  beginning  foam  thickness  is  great  enough  in  make 
lQ  lt0  S'*  .  the  reflected  wave  will  pass  through  the  rarefaction  fan 

Kqu«t"in  fll,  Sectic-n  2B. 
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issuing  toward  the  wall  from  the  locked  mass  at  time  t*.  and  the  second 
wave  on  the  wall  will  be  elastic  and  will  follow  a  relaxation  of  pressure. 
In  other  words,  for  l  slightly  higher  than  that  given  above,  the  locked 
mass  is  brought  to  rest  with  the  exertion  of  pressure  at  forerunner  stress 
or  below  on  the  wall.  The  complicated  calculation  needed  to  express  this 
value  of  l  exactly  is  not  justified  by  the  accuracy  of  the  model,  and  we 
consider  S  as  giving  the  approximate  "stopping  power”  of  the  foam  for  the 
locking  wave.  Furthermore,  if  we  write  as  a  first  approximation 


V 


t 


♦ 


where  K  probably  depends  somewhat  on  the  kind  of  material  making  up  the 
foam,  then  P  in  Eq.  (2)  becomes  the  same  for  all  materials  which  have 
the  same  proportionality  between  distention  and  reduction  of  elastic 
speed  due  to  distention.  On  this  basis,  then,  S  is  the  superficial  mass 
density  of  foam  needed  to  keep  pressure  on  the  structure  below  the  value 
given  by 

P  *  Pt(  1  ►  K) 


It  is  also  the  superficial  density  of  foam  at  which  increments  of  foam 
thickness  abruptly  become  more  effective  than  at  lower  densities.  Letting 


and 


we  write  stopping  power  5  as: 


5 


V 


\  *  d  -  -L 
%  /« 
1 


<  3) 


(The  asterisk  means  the  quantity  is  evaluated  under  conditions  such  that 
the  collision  of  the  two  locked  masses  takes  place  at  time  t*.)  Wc  ran 


Although  th#y  art  not  ro*(U*i*r  it  nur  remlti  in  rtprawH’  »Unnus  sj.-  nm  nakr  ih»«  a  arts 

fooi*  appro* tail  ion. 
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thus  express  the  stopping  power  as  a  function  of  three  independent  parameters 
from  the  simpl  i  fi  ed  equation  o  f  state,  namely  t'ff  Otand/fl.  If  Vf  and  V'0  are 
held  constant,  a  increases  wi  th  increasing  P r:  f 3,  the  only  variable  depending 
on  locked  volume  Vl ,  ri  ses  wi  th  growing  distentior.  l't  itself,  of  course, 
depends  both  on  the  slope  of  the  elastic  locus  and  the  initial  volume: 

■  * 

As  either  one  of  the  two  variables  a  and  ft  is  varied  through  the  range 
0  <  Ct  <  ft  and  V  held  constant,  the  right  side  of  Eq.  (3)  has  no  maximum; 
stopping  power  grows  without  limit  wi th  both  increasing  distention  and  in¬ 
creasing  elastic  pressure.  If  we  take  a'  *  (V0  -  and 

ft  =  ( V0  -  vt/v^  the  result  is  similar,  i.e., 

S  ■  r,  7-f-Z  • 

a'  ft' 


Other  things  being  equal,  stopping  power  increases  with  final  specific 
volume,  slope  of  the  elastic  segment  of  the  stress-strain  curve,  disten¬ 
tion  (l',,  -  V t ) / V j,  and  elastic  pressure. 

In  most  foams  increasing  distentions  generally  implies  falling  P  and 
modulus,  so  there  may  be  an  optimum  distention  which  must  be  found  from 
empirically  determined  relations  between  elastic  properties  and  distention. 
As  an  illustration  of  this  procedure  but  not  as  a  demonstration  of  actual 
optimization  in  polyurethane,  we  will  write  empirical  relations  between 
yield  stress  and  initial  density  and  between  modulus  and  density  based  on 
the  static  measurements  of  Lindberg,  Gates,  and  Baer22  in  lightweight 
polyurethane: 

P,  K  , 

—————  »  —  +  A 


where 


14.15  *  l:09cm2,'sec2 


~1.65  *  109dyne/cm2 


4.74  *  I0*cm2/sec2 


“0.674  H  10®Hvne/cm2 
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Equation  (3)  becomes: 


S 


We  seek  maxima  in  this  function  of  V’Q  so  that  P e  =  / V Q  +  B  <  Pfm,  or 

V'o  >  /(Pe B  -  B).  We  are,  of  course,  limited  to:  >  V’j. 

PeB  is  the  most  Pf  can  be,  if  the  second  pressure  jump  due  to  the  col¬ 

lision  at  t*  is  to  bring  the  peak  stress  in  the  structure  just  up  to  some 
previously  specified  limit.  If  ^  =  1.0  cm3/gandPe>(  =  1.0  kbar,  the  second 
inequality  controls,  and  in  the  range  1  <  V0  <  7  the  function  S  has  a  broad 
maximum  between  V  =2  and  ^.=  4.  Row  I  in  the  tabulation  shows  the  behavior. 


K0(c»3/g) 

1.0 

1.03 

1.5 

..0 

3.0 

4.0 

5.0 

6.0 

7.04 

I 

10  4S(cm/sec) 

0 

0.922 

1.105 

1.16 

1.08 

0.922 

0.680 

0 

II 

10  4S(cm/sec) 

00 

1.31 

1.57 

1.79 

1.90 

1.96 

1.99 

2.02 

III 

10  ^S( cm/sec) 

1.79 

1.77 

1.48 

Imaginary 

Setting  A  -  B  =  0  may  make  the  relations  between  distention  and  the  forerunner 
parameters  more  nearly  correct  and  has  the  e f f ec t  on  S  seen  i n  Row  II.  (The 
theory  and  constitutive  rel  ations  break  down  as  -*  1, 0.  )  Row  III,  which  also 
shows  little  influence  of  distention  on  S,  i  s ca 1 cul a  ted  f rom  a  different  set 
of  linear  relations  between  modulus  and  yield  stress  on  the  one  hand  and  K0  on 
the  other.  Here  quasi-dynamic  data  for  0. 67  and  0. 335  g/cc  polyurethane  is 
used,  specifically  P g  =  1.0  and  0.  15  kbar,  respectively,  and.  {/  *1.6 

and  1.3  mm//isec,  respectively.  These  val ues  o f  P(  ard  the  values  of  V  for 
0.67  g/cc  material  are  reported  by  Rempel 2  from  shock  measurements.  The  value 
of  U t  for  the  lighter  polyurethane  i  s  rel ated  to  quasi - stat i c  compression  i n for¬ 
mation  also  reported  by  Rempel .  Thur  4  =  “5.7  x 1 0 9 ,  B  -  "0.7  * 109  dyne  cm-2, 

K j  =  33.7  x  1 0 9 ,  and  Kz  =  2.54  *  109  cm  2  sec“2  for  Row  III. 


It  is  possible  foam  effectiveness  is  not  highly  sensitive  to  distention. 
Figure  4  seems  to  confirm  that  distention  over  the  range  used  in  our  experi¬ 
ments  has  relatively  small  bearing  on  pressure  attenuation  effectiveness  of 
aluminum  foam.  The  single  point  in  Fig.  5  from  0.335  g/cc  polyurethane 
(Shot  9228)  does  not  appear  to  be  in  good  agreement  with  this  conclusion. 
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CONCLUSIONS  AND  FUTURE  DIRECTIONS 


o . 


t 


Good  methods*  exist  to  estimate  shock  response  of  certain  elastic- 
rigid  foams  at  very  high  pressures  when  no  forerunner  is  present,  and  at 
very  low  pressures  when  no  locking  wave  is  present.  For  the  intermediate 
region  of  pressure  (approximately  0.  I  to  6  kbar),  when  both  waves  are 
present  and  of  the  same  order  of  magnitude,  predictions  are  less  reliable 
and,  in  our  experience,  values  forecast  of  pressure  in  this  region  based 
on  general  theories  are  too  high.  From  considerations  of  pressure  wave 
shape  we  believe  this  anomaly  may  be  related  to  the  porous  nature  of  the 
foam.  We  have,  however,  experimentally  measured  peak  pressures  trans¬ 
mitted  to  a  steel  wall  to  about  6  kbar,  found  conditions  for  double-wave 
structure,  and  observed  wave  shapes  produced  by  in. pacts  carrying  momentum 
densities  in  the  range  0.6  to  3  x  104  taps  in  0.7  to  1.4  g/cc  aluminum, 
0.67  g  fee  polyurethane,  1.1  g/cc  beryllium,  1.1  to  1.7  g/cc  graphite,  and 
1.0  g/cc  silica  foam. 


Predictions  of  pressure  histories  in  the  general  theories  are  based 
on  statically  measured  densities  and  dynamically  measured  forerunner 
speeds  and  strengths.  In  the  materials  studies,  forerunner  speeds  can 
be  related  to  statically  measured  moduli  within  ±  15  percent;  forerunner 
strengths  are  two  to  four  times  static  yields  and  have  not  been  related 
to  static  yields  by  theory. 


In  the  sense  of  transmitting  the  least  peak  pressure  for  the  least 
added  weight  under  impact  of  a  given  amount  of  momentum,  polyurethane 
appears  to  be  the  best  of  the  foams  studied  (polyurethane,  graphite, 
silica,  beryllium,  and  aluminum).  However,  the  momentum  gain  due  to  the 
bounce-off  of  a  polyurethane  laver  from  a  steel  wall  is  higher  than  that 
of  other  foams,  particularly  aluminum,  where  the  difference  in  total 
momentum  passed  to  the  structure  by  the  two  foams  may  be  as  large  as  50 
percent.  Moreover,  polyurethane  is  subject,  to  a  considerable  loss  of 


Ba  ftd  on  the  experimentally  established  Hugoniots  of  Fowles  and  Curran. 
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effectiveness  as  a  pressure  attenuator  at  100°C,  whereas  aluminum  per¬ 
formance  does  not  deteriorate  to  the  same  extent  until  300°C  is  reached. 
Static  data  in  heated  silica  point  to  a  possible  improvement  in  high 
temperature  performance. 

Improvements  in  the  prediction  system  will  most  probably  come  from 
the  following  studies:  (1)  establishment  of  the  pressure- volume  equation 
of  state  of  foams  in  the  neighborhood  of  the  elastic  limit  by  observing 
forerunner  impacts  on  a  series  of  anvils  of  rigidity  less  than  that  of 
steel,  (2)  removal  of  the  doubt  introduced  in  all  static-dynamic  correla¬ 
tions  by  variations  in  physical  properties  of  foam  specimens,  (3)  cataloging 
of  pressure  histories  in  one  material  over  a  wide  range  of  input  momenta 
and  foam  thickness  through  the  use  of  1-D  experiments,  and  finally  (4)  an 
attempt  to  interpret  these  histories  as  revealing  the  interaction  of  a 
wall-reflected  wave  known  from  (1)  above  with  the  main  wave  and  as  reveal¬ 
ing  the  nature  of  shock  flow  in  porous  material  (as  contrasted  to  a  material 
of  uniform  density).  Although  it  is  not  likely  that  standard  methods  of 
analysis,  such  as  the  method  of  artificial  viscosity,  will  be  useful  in 
revealing  the  nature  of  locking  shock  motion  at  low  pressures  in  porous 
materials,  exploration  by  the  method  of  artificial  viscosity  of  the 
effects  on  forecasted  main  wave  shapes  of  more  drastic  departures  from 
the  simple  equation  of  state  than  have  been  studied  so  far  could  be 
informative . 

Foam  forerunners  could  provide  fruitful  fields  for  investigation  of 
strain-rate  effects  and  the  phenomenon  of  stress  relaxation  discussed  by 
Duvall. 9 
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1 .  INTRODUCTION 


The  effects  produced  in  solids  by  events  such  as  the  impact  of  pro¬ 
jectiles  and  the  detonation  ol  high  explosives  and  nuclear  devices  ran  be 
studied  directly  by  performing  experiments.  in  such  studies  the  object 
may  be  to  observe  the  physical  response  of  materials  at  high  pressure. 
Experiments  for  collecting  data  on  the  equation  of  state  of  a  material 
are  an  important  example  of  such  work  and  are  usually  on  a  small  scale. 
Other  studies  are  of  less  fundamental  nature  and  are  done  to  improve  or 
to  test  the  design  of  a  structure  or  a  device.  for  example,  a  missile 
may  be  damaged  or  cestroyed  if  a  nuclear  device  is  detonated  in  its  vi¬ 
cinity.  In  such  :  case,  the  experimental  approach  is  at  least  expensive 
and  time  consuming,  if  not  virtually  forbidden.  Because  of  these  reasons, 
much  of  In?  study  of  the  vulnerability  of  missiles,  underground  structures, 
and  other  installations  must  be  performed  by  the  application  of  theory  to 
predict  shock  attenuation.  The  work  reported  here  is  a  continuation  of  an 
effort  to  verify  the  validity  of  the  application  of  hydrodynamic  theory  in 
such  studies1’'  and,  where  discrepancies  are  observed,  to  formulate  im¬ 
proved  models  that  will  provide  more  accurate  predictions.  The  term  hydro- 
dynamic  theory  is  taken  to  mean  a  theory  in  which  material  rigidity  is 
neglected  and  the  velocities  of  waves  in  the  media  are  hydrodynamic. 

At  very  high  pressures,  material  rigidity  is  expected  to  be  unimpor¬ 
tant  so  that  the  principal  stresses  are  approximately  equal,  and  the  hydro- 
dynamic  theory  is  a  good  approximation.  For  phenomena  associated  with  the 
shork  front  itself,  hydrodynamic  theory  seems  to  be  valid  at  much  lower 
pressures;  e.g..  in  determining  the  llugomot  equation  of  state.  As  a  shock 
front  propagates  from  its  source  it  is  attenuated.  In  this  report,  it  ts 
assumed  that  attenuation  is  a  wave  phenomenon  in  which  rarefaction  waves 
overtake  and  interact  w u h  the  shock  front.  If  this  attenuation  is  to  be 
calculated,  the  speed  of  these  relief  waves  must  be  known.  It  rannnt  be 
assumed  that  their  speed  is  always  independent  of  rigidity,  even  if  rtgidilv 
can  be  ignored  for  the  shock  transition.  An  objective  of  the  work  reported 
here  is  to  determine  if  rigidity  is  important  in  the  flow  behind  shocks  in 
some  representative  solids. 
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The  effects  of  rigidity  on  the  speed  of  relief  waves  can  be  studied 
by  measuring  directly  the  velocity  of  the  waves.  What  amounts  to  the  same 
thing  is  the  study  of  the  attenuation  of  shock  waves.  An  attempt  is  made 
in  this  work  to  study  both  the  velocity  of  the  waves  and  the  attenuation 
caused  by  them.  In  previous  work,  attenuation  of  shocks  was  studied  in 
aluminum,  copper,  rock  salt,  and  other  rocks.1,2  In  this  work,  more  obser¬ 
vations  have  been  made  for  both  aluminum  and  copper.  A  plastic 
(Armstrong  C-7  epoxy)  and  gold  have  also  been  studied.  Targets  made  of 
these  materials  are  struck  by  aluminum  plates  moving  with  a  velocity  of 
about  0.125  cm/Vtsec.  This  arrangement,  which  induces  an  essentially  plane 
wave  in  the  target,  is  one  of  the  si  iplest  for  which  flow  calculations  can 
be  readily  performed.  Comparison  of  the  results  of  calculations  and  ex¬ 
periments  determine  the  extent  of  agreement  with  the  hydrodynamic  theory 
in  each  of  the  several  different  situations. 


90 


2.  SUMMARY 


Attenuation  of  shock  waves  was  studied  in  1060  and  2024  aluminum, 

OFHC  copper,  gold,  and  Armstrong  C-7  epoxy.  Shock  waves  were  induced 
in  samples  of  these  materials  by  striking  them  with  aluminum  plates 
traveling  about  0.125  cm/jusec.  Because  the  strains  induced  in  the  speci¬ 
mens  were  virtually  one- dimensional ,  the  experiments  could  be  simulated 
with  a  computer  code.  Comparison  of  experimental  and  calculated  results 
indicated  that  an  elastopl astic  relation  between  stress  and  strain  should 
be  used  for  aluminum,  copper,  and  epoxy.  More  data  are  needed  before  a 
s t ress - s t rai n  relation  can  be  suggested  for  gold. 

The  amplitude  of  the  elastic  relief  wave  in  2024  aluminum  was  approxi¬ 
mately  20  kbar  when  the  maximum  stress  was  120  kbar.  On  the  assumption 
that  Poisson’s  ratio  is  constant,  a  yield  strength  value  of  5  kbar  was 
inferred  (see  Sections  4  and  6).  This  is  the  yield  strength  in  simple 
tension,  or  twice  the  critical  resolved  shear  stress.  For  10b0  alum¬ 
inum,  the  amplitude  of  the  elastic  relief  wave  was  16  kbar,  so  that  the 
yield  was  4  kbar  when  the  maximum  stress  was  120  kbar.  The  relief  wave 
in  copper  had  an  amplitude  of  27  kbar,  so  that  the  yield  was  app roxi ma t e 1 v 
7  kbar  for  a  stress  of  160  kbar.  These  values  of  the  yield  strengths  are 
all  greater  than  those  observed  at  engineering  stress  levels  and  imply  that 
the  yield  strength  is  a  function  of  the  stress. 

It  was  impossible  to  obtain  u  value  of  the  yield  strength  directly 
from  the  experimental  data  for  either  gold  or  the  epoxy.  For  the  epoxy, 
the  calculated  attenuation  matched  the  experimental  observations  fairly 
well  when  the  yield  was  assumed  to  be  1.0  kbar. 

Two  methods  were  used  for  recording  the  response  of  materials  to 
shock  waves.  One  employed  the  smear  cameru  to  observe  the  f ree- surface 
motion  of  specimens  shaped  like  wedges.1  This  technique  has  been  used  for 
a  number  of  years  and  gives  reliable  results.  As  ordinarily  used,  the 
technique  records  only  the  velocity  of  the  free  surface  at  the  instant  the 
shock  interacts  with  the  surface.  Its  use  raises  the  question  us  to  the 
relation  of  the  free-surface  velocity  to  the  particle  velocity  induced 
by  the  shock  front. 
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In  an  attempt  to  circumvent  some  of  the  disadvantages  of  the  optical 
technique,  a  second  method  of  recording  response  by  means  of  a  transducer 
was  used  in  several  experiments.  The  device  consisted  of  a  Manganin  wire 
potted  in  Armstrong  C-7  epoxy. 4  It  was  calibrated  to  measure  the  pressure 
induced  in  the  epoxy  by  the  transmission  of  a  shock  from  a  specimen  into 
the  epoxy.  Although  the  transducer  is  still  in  the  developmental  stage,  it 
was  hoped  that  it  would  permit  the  recording  not  only  of  peak  pressure  but 
also  of  pressure  as  a  function  of  time  behind  the  shock  front.  Measure¬ 
ment  of  the  peak  pressure  was  as  good  as  can  be  expected,  although  some 
transducers  gave  obviously  erroneous  results  for  reasons  that  are  not 
understood.  The  transducers  do  not  give  the  expected  response  after  the 
passage  of  the  shock  front,  perhaps  because  of  a  hysteresis  effect  which 
is  being  investigated. 5  Presently,  the  transducers  appear  incapable  of 
giving  the  amplitude  and  arrival  time  of  the  elastic  relief  waves  which 
are  implied  to  exist  in  the  observations  obtained  with  the  smear  camera. 

Detailed  behavior  of  an  explosive  system  such  as  is  used  to  throw 
aluminum  plates  is  not  understood.  Some  studies  of  the  system  are  re¬ 
ported  here,  but  they  ^re  far  from  being  exhaustive  stidies  even  on  the 
one  system  which  throws  aluminum  plates  at  a  velocity  of  0.125  cm/gsec. 
Enough  data  are  given  to  show  that  the  flyer  plates  are  probably  not 
spalling  and  that  their  shapes  are  tolerable.  During  preliminary  studies, 
reproducibility  of  the  velocities  of  flyer  plates  did  not  appear  to  be  a 
problem.  It  became  a  problem  later  in  the  attenuation  experiments  in 
which  Baratol*  pads  were  obtained  on  a  new  order.  Much  more  work  should 
be  done  on  methods  of  using  explosives  for  throwing  flat,  strain-free, 
unspalled  aluminum  plates  not  only  at  0.125  cm/Vsec,  but  at  greater 
velocities;  methods  of  throwing  and  using  plates  made  of  other  materials 
should  also  be  investigated. 


*  Btrttol  it  67  ptrctnt  btrma  turtle  ttd  33  parettt  TNT. 
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3.  EXPERIMENTAL  ARRANGEMENT  FOR  STUDYING 
ATTENUATION  OF  SHOCKS  IN  SOLIDS 


A.  FLYER  PLATE  ASSEMBLIES 

All  the  attenuation  experiments  were  performed  by  throwing  aluminum 
plates  with  charges  of  high  explosives.  The  irrangement  used  in  the  ex¬ 
periments  is  represented  in  Fig.  1,  which  shows  a  plane-wave  generator 
(PWG)  and  an  explosive  pad  in  contact  with  a  brass  plate.  For  most  of 
the  work  the  explosive  pad  was  Baratol.  A  thermosetting  plastic  cement 
about  0.001  inch  thick  affixed  the  aluminum  plate  to  the  brass  plate. 

The  bra  ss  plate  was  used  partly  because  it  held  back  the  smoke  from  the 
explosive  and  permitted  a  clear  view  of  the  experiment,  and  because  it 
also  served  as  a  base  for  the  vacuum  chamber.  The  curvature  induced  in 
the  plate  when  the  chamber  was  evacuated  was  negligible.  The  functions 
mentioned  above  could  also  be  achieved  by  using  plates  of  materials  other 
than  brass — for  example,  aluminum.  However,  there  had  to  be  some  separa¬ 
tion  between  the  aluminum  flyer  plate  and  the  brass  buffer  plate  so  that, 
during  the  time  of  observation,  the  target  responded  only  to  the  impact 
of  the  aluminum  flyer  plate.  The  impedance  mismatch  between  aluminum  and 
brass  is  such  that  the  aluminum  separates  from  a  brass  buffer  plate  and 
travels  with  a  considerably  higher  velocity  than  the  brass. 

The  1  '8- inch-thick  aluminum  flyer  plates  attained  a  velocity  of 
0.  125  cm  i. sec  when  the  explosive  pad  was  Baratol.  Greater  velocities  can 
be  attained  by  the  use  of  more  energetic  explosives,  such  as  Composition  B- 3 
or  KMX .  Use  of  these  explosives  may  introduce  other  problems,  however. 

One  of  the  desirable  features  of  explosively  driven  plates  is  that  the  plate 
have  a  large  useful  diameter.  fchen  an  8- i ndi-d i amrt er  plane-wave  generator 
is  employed,  the  useful  diameter  of  the  flyrr  plate  may  be  as  large  as 
I  inches. 

If  the  strain  induced  in  a  target  by  the  impact  of  a  high-speed  plate 
is  to  be  one-dimensional ,  the  plate  should  be  flat  and  oriented  parallel  to 
the  surface  of  I  he  target  at  the  time  of  impact.  This  is  difficult  to 
achieve,  especially  when  the  fiver  plate  is  as  much  as  4  inches  in  diameter. 
Several  shots  were  fired  to  determine  an  acceptable  method  of  driving  plates 
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CAP  AND  BOOSTER 


FIG.  1  FLYER  PLATE  ASSEMBLY 


for  tb*  attenuation  studies.  In  these  shots,  a  flat  glass  mirror  was 
held  about  1/2  inch  from  and  parallel  to  the  initial  position  of  the 
aluminum  flyer  plate.  The  shot  was  suitably  illuminated  and  the  record 
of  the  collision  was  taken  by  the  smear  camera. 

Results  of  six  of  these  shots,  where  the  deformation  is  given  as  a 
function  of  the  diameter  of  the  plate,  are  given  in  Fig.  2.  The  scale 
for  the  deformation  is  enlarged  with  respect  to  the  ordinate.  Each  plate 
should  be  visualised  as  moving  from  right  to  left  so  that  the  renter  of 
the  plate  is  the  last  part  to  reach  a  targrt.  In  Shots  9886,  9887,  and 
9902,  the  original  diameter  of  the  aluminum  plates  is  4‘y  inches. 

Shot  9902  shows  a  fairlv  good  plate,  the  center  cf  which  lags  the  most 
forward  part  by  0.015  inch.  Berau.ie  the  velocity  of  the  piste  was 
0.125  rmVeer,  the  plate  first  struck  the  edges  of  the  target  and  then, 
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PLATE  CENTER  - inch** 


approximates  0.3  >.*vr  later,  struck  t  lie  renter.  I’lse  uselul  lii.imeter  <>l 
the  plate  was  under  4  inches.  Hesults  i>f  Shots  988')  and  9887  ate  eusc 
in  that  the  delay  of  the  center  was  ([renter  and  the  uselul  diameter  was 
I  e  s  s  . 

In  orde  to  obtain  plates  with  larger  useful  diameters.  t»  -  i  n  >  h  - 
diameter  plates  were  used  in  the  following  three  experiments.  Ihe  lag 
of  the  renters  in  these  shots  was  determined  h\  drawing  a  cloud  1  mi  lies 
long  across  t  fie  profile.  This  chord  represented  the  maximum  dimension 
of  the  target  for  most  experiments.  In  Shot  9938,  the  center  lag  was 
about  0.011  inch  or  0.2  /.sec;  in  Shot  10,017,  0.013  in<h  or  0.3  ..sec;  and 
in  Shot  10,018,  0.015  inch  or  0.3  /.sec.  hath  of  these  shots  gave  delay 


DEFORMATION  Of  PLATE  -  iftClt** 


FIG.  2  RESULTS  OF  FLYER  PLATE  PLANARITY  TESTS 
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times  for  the  center  less  than  or  as  good  as  the  times  for  Shot  9902,  the 
best  of  the  shots  using  4%-  inch-diameter  plates.  As  expected,  the  useful 
diameter  of  the  plates  was  improved  also.  No  great  advantage  was  given  in 
Shot  10,017  by  using  2-inch-thick  Baratol.  The  combination  of  an  8-inch- 
diameter  PWG  and  a  1-inch-thick  pad  of  Baratol  was  used  as  a  reasonable 
compromise  in  the  attenuat ion  experiments  which  followed.  Heproduc i bi 1 i ty 
of  such  a  system  is  inferred  by  comparing  the  duplicate  Shots  9938  and 
10,018,  the  lag  times  of  which  were  0.2  and  0.3  /xsec,  respectively. 

Previous  studies  showed  that  aluminum  plates  projected  by  a  driving 
assembly  containing  a  pad  of  Composition  B- 3  did  spall.  Those  studies 
involved  the  recovery  of  flyer  plates  as  well  as  the  use  of  the  flash 
X-ray  camera.  A  description  of  the  mechanism  which  produces  the  spall  is 
given  in  Appendix  I.  The  event  has  been  simulated  by  the  use  of  a  com¬ 
puter  code,  and  the  results  based  on  a  fracture  strength  of  30  kbar  indi¬ 
cate  that  the  aluminum  plate  should  spall. 

At  the  beginning  of  this  project  it  was  not  known  if  the  use  of 
Baratol  between  the  PWG  and  the  brass  plate,  Fig.  1,  would  cause  the 
aluminum  flyer  plate  to  spall.  Previous  experience  showed  that  the  prob¬ 
lems  associated  with  the  recovery  of  the  flyer  plates  after  they  had  been 
accelerated  to  a  high  velocity  by  an  explosive  charge  were  many,  and  that 
misleading  results  could  be  obtained.  Three  methods  were  used  for  deter¬ 
mining  if  the  flyer  plate  had  spalled 

The  first  method  used  the  high-speed  framing  camera.  The  flyer  plate 
was  thrown  across  the  field  of  view  of  the  camera,  and  the  event  was  back¬ 
lighted.  Although  great  attention  was  given  to  the  control  of  the  deto¬ 
nation  produce  gases,  no  useful  records  were  obtained  because  of  obscuration 
by  these  gases. 

Another  method  for  determining  spallation  employed  the  flash  X-ray 
camera.  A  double  flash  unit  was  used,  and  eacli  unit  was  flashed  at  pre¬ 
determined  times  after  the  flyer  plate  was  put  into  motion.  These  records 
(see  Fig.  3)  show  that  the  flyer  plates  remained  intact  when  Baratol  was 
used  as  the  explosive  (or  a  bare  plane-wave  generator  which  used  Baratol 
as  the  low-velocity  explosive).  As  a  test,  a  shot  was  also  fired  using 
Comp  B-3,  and  a  clear  indication  of  spalling  was  obtained,  as  shown  in 
Fig.  4. 

The  third  method  used  a  pressure  transducer  which  is  being  developed 
at  Stanford  Research  Institute.4  It  was  reasoned  that  if  the  plate  spalled 
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FIG.  3  X-RAY  PICTURES  OF  FLYER  PLATES  IN  MOTION 
(a)  Shot  9813  —  6-inch  diam  PWG 
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FIG.  3  X-RAY  PICTURES  OF  FLYER  PLATES  IN  MOTION 
(b)  Shot  9767  —  6-inch  diam  PWG  +  1-inch  Baratol 
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FIG.  4  X-RAY  PICTURE  OF  SPALLED  FLYER  PLATE  (Shot  9768) 


before  hitting  the  transducer,  the  transducer  record  should  show  a  second 
shock.  When  Baratol  was  used  in  the  explosive  train,  the  transducer 
record  showed  only  one  pulse.  For  the  shot  in  which  Comp  B- 3  was  used, 
the  record  did  show  a  second  increase  in  pressure.  Hence  the  transducer 
confirmed  the  X-ray  record,  because  both  instruments  were  used  simulta¬ 
neously,  as  indicated  in  Fig.  3(a). 

The  plate  thrown  by  Baratol  had  a  velocity  of  approximately 
0.125  cm/ftsec,  which  is  adequate  for  many  experiments.  Use  of  Comp  B- 3 
resulted  in  flyer  plate  velocities  of  0.19  cm/ftsec.  Such  flyer  plates 
would  be  useful  if  some  way  can  be  developed  to  prevent  the  plates  from 
spalling.  A  possible  method  is  to  introduce  a  layer  of  plastic  between 
the  brass  plate  and  the  aluminum  flyer  plate.  This  system  was  simulated 
by  a  computer  code  with  encouraging  results.  The  calculated  results  also 
indicate  that  the  plastic  will  not  seriously  interfere  with  the  experiment. 
The  computed  results  should  be  confirmed  by  experimental  results  before 
the  system  is  used. 

B.  METHODS  OF  RECORDING  FOH  ATTENUATION  EXPERIMENTS 

The  smear  camera  was  used  in  previous  work  to  record  the  free-surface 
motion  of  specimens  in  which  attenuating  shock  waves  were  being  studied. 

A  simple  way  of  recording  the  free-surface  velocity  was  to  mount  a  small, 
partially  reflecting  mirror  close  to  the  surface  being  watched.  Both  the 
mirror  and  the  surface  beneath  it  were  illuminated.  When  the  shock  ar¬ 
rived  at  the  surface  of  the  specimen,  the  reflectance  of  the  surface 
changed,  and  the  arrival  time  of  the  shock  was  thereby  recorded  on  the 
film  in  the  camera.  Because  the  surface  was  accelerated  by  the  shock,  it 
soon  crossed  the  gap  and  collided  with  the  mirror,  which  caused  another 
change  in  the  intensity  of  the  light  reaching  the  smear  camera.  Thus  the 
flight  time  of  the  free  surface  of  the  specimen  across  the  gap  was  re¬ 
corded.  This  method  of  obtaining  the  velocity  of  a  free  surface  is  not 
satisfactory  when  the  shock  wave  is  nonuniform.  In  such  a  case,  the  free 
surface  does  not  stay  in  uniform  motion,  so  that  some  average  velocity  is 
obtained  by  this  technique.  The  gap  technique  was  employed  with  specimens 
cut  in  the  shape  of  wedges  so  that  any  variation  of  the  free-surface  ve¬ 
locity  with  the  distance  of  travel  of  the  shock  in  the  specimen  could  be 
observed.  The  results  of  such  experiments  are  still  regarded  as  being 
fairly  reliable. 
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More  reliable  results  should  be  obtainable,  however,  with  the  optical 
lever  arm  technique,  in  which  targets  were  in  the  shape  of  wedges  and  the 
surface  facing  the  camera  was  polished  to  a  mirror  finish.3  This  require¬ 
ment  restricts  the  technique  to  homogeneous  materials  whose  surface  takes 
a  polish.  The  camera  was  not  focused  on  the  surface  of  the  specimen  but 
was  focused  on  the  light  source  reflected  in  the  mirror  surface  of  the 
speiimen.  The  light  source  was  covered  with  a  mask  with  alternate  trans¬ 
parent  and  opaque  spaces  (Fig.  5).  As  the  shock  wave  hit  the  free  surface 
of  the  specimen,  the  surface  was  turned,  and  the  image  of  each  grid  line 
was  displaced.  As  shown  in  Appendix  II,  the  displacement  of  the  grid  line 
and  the  velocity  of  the  shock  along  the  inclined  surface  of  the  specimen 
give  the  particle  velocity  induced  in  the  specimen  by  the  shock  front. 

The  advantage  of  this  method  is  that  the  velocity  of  the  surface  is  re¬ 
corded  at  the  appropriate  time — at  the  time  of  breakout  of  the  shock  front, 
and  thus  there  is  no  averaging,  as  was  mentioned  above  for  the  gap  method. 

A  disadvantage  is  the  necessity  for  differentiating  the  curve  defined  by 
the  termini  of  the  straight  lines  on  the  record,  see  Fig.  6.  If  the  flyer 
plate  were  perfectly  flat  and  were  oriented  parallel  to  the  target  at  the 
time  it  hits,  the  curve  described  above  would  have  no  inflection  points; 
that  is,  the  shock  front  would  travel  along  the  inclined  surface  of  the 
wedge  (see  Fig.  5)  at  constant  velocity,  or  at  a  continuously  decreasing 
velocity.  Curvature  of  the  flyer  plate  can  cause  an  inflection  in  the 
curve,  so  that  the  apparent  velocity  increases  near  the  thick  end  of  the 
wedge,  even  though  the  pressure  carried  by  the  shock  is  being  attenuated. 
This  inflection  makes  it  difficult  to  differentiate  the  curve  and  is  prob¬ 
ably  the  chief  source  of  error  in  the  analysis  of  the  records.  Otherwise, 
the  analysis  requires  knowledge  of  the  camera  writing  speed,  which  is 
known  to  a  precision  of  about  0.5  percent;  of  the  camera  magnification, 
which  is  determined  from  each  record  to  a  precision  of  about  1  percent; 
and  of  the  displacement  of  each  trace,  which  is  measured  to  about  5  percent. 

Both  the  gap  method  and  the  optical  lever  arm  method  give  information 
about  the  free  surface  of  the  target.  It  is  usually  assumed  that  the  free- 
surface  velocity  is  twice  the  velocity  of  the  particles  immediately  behind 
the  shock  front  prior  to  the  time  of  incidence  of  the  shock  on  the  surface. 
This  assumption  was  used  in  the  analysis  of  the  records  because  it  simpli¬ 
fied  the  calculations.  The  error  introduced  was  of  the  same  order  as  the 
experimental  error. 
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FIG.  5  WEDGE  AND  LIGHT-SOURCE  COVER  MOUNTED  ON  VACUUM  CHAMBER 


FIG.  6  SMEAR  CAMERA  RECORD  OF  AN  OPTICAL  LEVER  SHOT 


A  transducer  was  also  used  in  studying  th*>  attenuation  of  shock  waves. 
The  device,  which  was  developed  at  Stanford  Research  Institute,  consisted 
of  a  Manganin  wire  cast  in  C- 7  epoxy  manufactured  by  Armstrong  Products 
Company.4  Prior  to  the  recording  time,  a  constant  current  was  passed 
through  the  Manganin  wire.  As  the  wire  was  compressed  by  the  shock,  its 
resistance  changed.  Signal  leads  were  connected  to  the  wire,  Fig.  7,  so 
that  the  resistance  change  could  he  monitored  by  an  oscilloscope  as  a 
change  in  voltage.  The  device  was  calibrated  to  a  pressure  of  about 
200  kbar,  and  the  relation  between  pressure  and  the  resistance  was  found 
to  be  linear.  The  calibration  experiments  depended  on  the  use  of  the 
smear  camera  to  determine  the  input  to  the  transducer  and  on  knowledge  of 
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(b)  PICTORIAL  VIEW 


FIG.  7  CONSTRUCTION  OF  MANGANIN  WIRE  TRANSDUCER 

the  equation  of  state  of  the  epoxy  and  of  the  material  used  to  transmit 
the  shock  into  the  transducer. 

The  response  of  the  transducer  was  adequate  to  record  accurately  the 
arrival  of  the  shock  wave.  The  transducer  also  appeared  to  respond  to 
changes  in  pressure  after  the  passage  of  the  shock  front;  this  is  the 
chief  attraction  in  the  use  of  transducers  in  the  study  of  attenuating 
shock  waves.  Accurate  records  showing  pressure  as  a  function  of  time 
would  be  helpful  in  determining  the  velocity  of  the  initial  relief  wave 
behind  a  shock  front.  If  the  relief  wave  is  elastic,  transducer  records 
should  be  useful  in  measuring  the  strength  of  the  wave.  These  data  would 
be  helpful  in  formulating  models  for  the  relation  between  stress  and 
strain  for  the  solids  being  studied. 

Application  of  the  transducer  to  the  study  of  attenuation  of  shock 
waves  in  solids  was  accomplished  by  placing  the  transducer  on  the  surface 
of  the  specimen  being  studied.  The  Manganin  wire  was  separated  from  the 
interface  by  about  0.1  cm  of  the  epoxy.  In  order  to  expedite  the  work, 
J.hree  transducers  were  mounted  together,  'as  shown  in  Fig.  8(a).  One 
transducer  was  placed  at  the  surface  of  the  specimen.  The  other  two  were 
placed  in  holes  machined  into  the  specimen.  This  was  relatively  easy  to 
accomplish  because  the  wires  could  be  positioned,  and  the  epoxy  could 
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FIG.  8  ARRANGEMENT  OF  THREE  TRANSDUCERS  FOR  A  SINGLE  SHOT 
(o)  Top  Viow 
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FIG.  8  ARRANGEMENT  OF  THREE  TRANSDUCERS  FOR  A  SINGLE  SHOT 
(b)  Bottom  Vi*w  —  Assembly  Roody  for  Mounting  of  Explosive 

then  be  poured  to  the  dean  red  tlut  knr.u.  The  hardened  epoxv  was  mi  f  f  i  * 
rientlv  rapid  mi  that  it  and  the  sjirt  narn  formed  the  lop  of  the  vaiuum 
chamber.  Figute  8(b)  shown  an  assembly  ready  for  loading  with  explosive. 

Each  transducer  was  connected  to  two  oscilloscopes.  The  primary 
oscilloscope  was  set  up  to  record  the  signal  induced  in  the  transducer  bv 
the  application  of  a  sourer  of  constant  current  to  the  transducer.  This 
current  was  turned  nn  a  lew  mirroseconds  before  the  shock  front  reached 
the  transducer.  The  primary  scope  also  showed  the  signal  due  to  the  change 
of  resistance  of  the  Manganin  wire  as  it  *a»  compressed.  This  stratagem 
made  it  unnecessary  to  measure  either  the  resistance  of  the  transducer  or 
the  current  induced  in  it  b>  the  power  supply;  it  was  only  necessarv  to 
assume  that  the  power  supply  gave  a  constant  current.  Any  change  in  the 
voltage  output  of  the  gage  due  to  the  application  of  pressure  to  the  wire 
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is  proportional  to  t  In*  <  h  ailin'  in  tr^isiain  r  of  t  li  <■  gage .  !!••  to  »•  ,i  i  n  mil 

such  as  that  shown  in  hig.  9  was  analv/ed  l>\  fir.-t  him  -•  u  >  t  ng  tin-  del  lei  - 
t  i  on  4  to  B,  a  nil  then  a  punt  i  ou  I  ii  lx-  tea  if  oil  t  li  ••  s  i  un.i  I  i  i  >  u  I  t  i  ng  from 
compressing  the  transducer,  >m  li  as  point  .  The  Mtm  of  IW.  to  \B  g  a  v  e 
the  relative  change  of  resistance  of  the  I  i  iiiimIiii  ct  .  \  .1  i  1 1>  1  a  t  1  on  <  u  r  v  e 

was  then  used  to  ofitain  the  pressure. 

Because  the  record  of  the  signal  given  t>\  the  priitiarv  scope  was 
relatively  small,  a  set  ondu  r  \  scope  was  used  to  d  1  sp  i  .1  \  the  signal  after 
it  had  been  amplified.  This  gave  .1  reiord  which  could  he  ana)  wed  more 

accurately. 


FIG.  9  TYPICAL  TRANSDUCER  RECORD 
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4.  RESULTS  OF  WEDGE  EXPERIMENTS 


A.  C- 7  EPOXY  WEDGE-SHAPED  SPECIMEN 

The  optical  lever  arm  technique  which  requires  the  use  of  a  specimen 
in  the  shape  of  a  wedge  is  described  in  Section  3B.  A  record  for  such  a 
shot  is  shown  in  Fig.  6  and  the  complications  introduced  in  the  analysis 
due  to  the  curvature  of  the  flyer  plate  are  also  discussed  in  Section  3B. 
This  record  was  further  complicated  because  of  a  slight  bend  in  the  wedge 
near  the  apex.  The  specimen,  made  of  C-7  epoxy,  was  mounted  in  the  cover 
of  the  vacuum  chamber  so  that  it  had  to  support  the  force  exerted  by 
normal  atmosphere.  This  force  produced  in  the  thin  portion  of  the  wedge 
some  deformation  which  made  the  record  somewhat  unreliable  near  the  symbol 
A  in  Fig.  6-  Otherwise  the  data  obtained  from  the  record  were  readily 
reducible. 

The  curve  from  A  to  B  in  Fig.  6  is  a  record  of  the  travel  of  the 
shock  front  along  the  inclined  surface  of  the  specimen.  In  order  to  ob¬ 
tain  the  velocity  of  the  shock  along  this  surface,  the  curve  must  be 
differentiated.  This  was  done  by  reading  the  time  and  distance  coordinates 
of  the  end  of  each  undeflected  trace  in  a  Vanguard  Motion  Analyzer  (VMA). 
Where  lines  were  omitted  on  the  grid,  synthetic  data  were  generated  by 
averaging  the  results  of  the  two  adjacent  lines.  The  data  were  then 
differenced  in  order  to  determine  if  any  errors  are  detectable,  and  any 
obvious  errors  were  corrected.  Next,  the  data  were  smoothed  by  use  of  a 
5-point  smoothing  formula,  following  which  a  5-point  differentiation 
formula  was  used  to  obtain  the  slope.  Application  of  the  appropriate 
reduction  factors  then  gave  the  velocity  of  the  shock  along  the  surface 
of  the  wedge. 

The  displacement  of  the  traces  were  also  read  on  the  VMA,  and  these 
data  were  also  differenced  and  smoothed  in  order  to  minimize  errors.  The 
method  for  obtaining  the  shock - i nduc ed  particle  velocity  is  discussed  in 
Appendix  II.  Peak  pressure  could  then  be  determined  from  Hugoniot  data 
for  the  material  concerned. 
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PRESSURE -  kbor 


FIG.  10  PRESSURE  IN  C-7  EPOXY  vs.  x/xQ  (xQ  is  thickness  of  flyer  plate) 

Results  for  the  C-7  epoxy  wedge  are  given  in  Fig.  10,  which  shows 
the  shock  -  induced  particle  velocity  as  a  function  of  the  distance  into 
the  C-7.  The  record  was  unreliable  near  the  thin  end  of  the  wedge,  as 
noted  above,  so  no  reliable  values  of  the  particle  velocity  were  obtained 
for  values  of  x/x 0  less  than  approximately  3.*  An  impedance  "  mi smatch ” 
solution  gave  a  pressure  of  about  46  kbar  in  the  absence  of  any  attenuation. 


*  It  wee  convenient  to  divide  the  diatence  the  ahock  had  traveled  into  the  apecimen,  x,  by  the  thickneas 
of  the  flyer  plate,  x^,  ehich  waa  0.322  c»  for  moat  experiaenta  reported  here. 


110 


Hence  attenuation  of  the  shock  commenced  at  a  distance  of  about  2.5  times 
the  thickness  of  the  aluminum  flyer  plate.  The  hydrodynamic  calculation, 
which  will  be  described  later,  predicted  that  attenuation  started  at  4.9. 

B.  ALUMINUM  WEDGE-SHAPED  SPECIMEN 

Specimens  were  fabricated  from  both  2024-T4  and  1060  aluminum.  After 
being  machined,  the  2024  specimens  were  heat  treated  in  an  air  atmosphere 
at  775°F  for  three  hours,  and  were  then  cooled  at  500°F  at  a  rate  of  50°F 
per  hour.  Specimens  cut  from  1060  aluminum  were  held  at  a  temperature  of 
650°F  for  one-half  hour. 

Smear  camera  records  showing  the  free-surface  motion  of  aluminum  and 
copper  are  shown  in  Figs.  11(a)  and  (b).  The  undeflected  traces  did  not 
terminate  along  a  monotonic  curve,  as  was  observed  for  the  C-7  wedge. 

That  the  record  for  copper  showed  less  curvature  than  the  record  for 
aluminum  was  accounted  for  on  the  basis  that  shock  velocity  is  less  in 
copper  than  in  aluminum.  When  a  curved  flyer  plate  hits  a  target,  a 
phenomenon  occurs  which  is  analogous  to  refraction  in  optics.  The  shock 
is  curved  more  than  the  flyer  plate  if  the  shock  velocity  is  greater  than 
the  velocity  of  the  flyer  plate.  Hence  analyses  of  wedge  experiments  be¬ 
come  easier  for  materials  having  low  shock  velocities.  Conversely,  greater 
flyer  plate  velocities  give  records  which  are  more  easily  analyzed 
(assuming  that  the  shock  velocity  does  not  increase  as  fast  or  faster  than 
the  flyer  plate  velocity). 

It  should  be  noted  that  all  records  are  not  as  clear  as  those  shown 
in  Fig.  11.  For  some  experiments  the  displaced  traces  were  so  smeared 
that  individual  traces  could  not  be  recognized,  and  measuring  the  trace 
displacement  was  impossible.  The  cause  of  such  smearing  is  not  clear. 

It  has  been  attributed  to  a  poor  polish  on  the  face  of  the  specimen  or  to 
residual  polishing  material  on  the  surface.  Errors  in  focusing  the  camera 
may  be  responsible,  or  curvature  of  the  displaced  surface  may  depend  in  a 
critical  manner  on  the  wedge  angle. 

Results  were  obtained  from  the  record  for  the  metal  wedges  in  the 
same  way  as  described  for  the  C-7  epoxy  wedge  above.  Results  for  two 
shots  using  annealed  2024  aluminum  target  wedges  are  shown  in  Fig.  12(a). 
Because  what  appeared  to  be  very  early  attenuation  was  observed  in  earlier 
work,  Shot  10,227  was  designed  to  give  data  at  a  few  plate  thicknesses 
into  the  target.  This  was  accomplished  by  using  a  small  wedge  angle. 
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F!G.  11  RECORD  OF  MOTION  OF  FREE  SURFACES  OF  ALUMINUM 
AND  COPPER  WEDGES 
(b)  Copper  (Shot  10,366) 
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This  shot  shows  what  can  be  called  attenuation  at  about  5  plate  thick¬ 
nesses;  at  3  plate  thicknesses  the  results  compare  favorably  with  the 
hydrodynamic  solution  which  will  be  described  later  (in  this  region  it 
amounts  to  an  impedance  mismatch  solution  using  only  Hugoniot  data). 

The  apparent  attenuation  exhibited  by  the  data  from  Shot  10,227  did 
not  appear  in  the  earlier  analysis.7  There,  differentiation  was  performed 
graphically,  and  apparently  that  method  smoothed  the  data  to  an  extent 
greater  than  the  numerical  method  employed  later. 

Because  no  attenuation  was  observed  in  the  preliminary  analysis  of 
Shot  10,227,  another  experiment,  Shot  10,354,  was  designed  to  yield  data 
for  greater  values  of  */*0-  The  results  are  given  in  Fig.  12(a);  the 
particle  velocity  is  surprisingly  high  at  values  of  x/x 0  near  5.  All 
shots  are  designed  on  the  assumption  that  the  explosive  systems  give  re¬ 
producible  flyer  plate  velocities,  which  are  measured  in  separate 
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FIG.  12  PEAK  "ARTICLE  VELOCITY  IN  ALUMINUM  v«.  x/*c 
(i)  2024  Aluminum 


experiments.  Plate  velocities  have  been  observed  in  attenuation  experi¬ 
ments  at  a  cost  of  more  difficulty  in  design  of  the  experiment  and  in 
less  satisfactory  records  of  the  motion  of  the  free  surface  of  the  speci¬ 
men.  The  present  results  show  that  reproducibility  is  not  good  when 
Baratol  is  used.  The  average  of  the  four  points  near  5  plate  thicknesses 
gives  a  particle  velocity  of  close  to  0.068  cm/gisec.  If  no  attenuation 
at  this  thickness  is  assumed,  a  flyer  plate  velocity  of  0.136  cm/gsec  is 
implied.  In  a  later  experiment,  it  was  observed  that  an  0  040- i nch - th i ck 
aluminum  plate  achieved  a  velocity  of  0.133  cm/gtsec.  Hence  the  unex¬ 
pectedly  high  values  of  particle  velocities  in  Shot  10,354  may  be  due  to 
the  inadvertent  use  of  a  charge  of  Baratol,  which  was  capable  of  giving 
greater  velocities. 

In  Shot  10,354  attenuation  commenced  between  5  and  6  plate  thick¬ 
nesses,  in  fair  agreement  with  Shot  10,227.  Averaging  the  first  four 
points  for  Shot  10,354  gives  a  particle  velocity  of  0.068  cm/gxsec,  which 
corresponds  to  a  stress  of  120  kbar  on  the  Hugoniot  for  aluminum.  At 
x/*0  approximately  8,  the  particle  velocity  is  approximately  0.058,  and 
the  stress  is  100  kbar.  If  this  attenuation  is  entirely  due  to  the 
elastic  relief  wave,  the  value  of  the  yield  strength,  Y,  is  approximately 
5  kbar,  or  about  twice  that  expected  for  2024  aluminum  under  more  ordinary 
conditions.  The  relation  between  the  yield  strength  and  the  amplitude  of 
the  elastic  relief  wave  is  given  in  Section  6. 

Besults  of  two  shots  using  annealed  1060  aluminum  are  shown  in 
Fig.  12(b).  Attenuation  in  Shot  10,226  commenced  at  5  or  6  plate  thick¬ 
nesses.  The  results  agree  well  with  the  hydrodynamic  solution  for  small 
values  of  x/xQ  also.  Results  of  Shot  10,353  show  unexpectedly  high  values 
of  the  particle  velocity  at  about  5  plate  thicknesses.  As  before,  this 
must  indicate  that  the  flyer  piate  had  a  velocity  greater  than  expected — 
0.125  cm/gsec.  If  this  is  true,  attenuation  in  1060  aluminum  commences 
at  about  5  or  6  plate  thicknesses,  just  as  it  does  in  2024  aluminum.  The 
data  are  not  precise  enough  to  determine  to  a  fractional  part  of  a  fiver 
plate  thickness  the  point  at  which  attenuation  commences. 

In  Shot  10,353,  the  particle  velocity  drops  from  about  0.067  to 
0.059  cm/gisec  in  the  vicinity  of  x  xQ  of  about  7.  By  using  the  same 
Hugoniot  as  was  used  above  for  2024  aluminum,  the  stress  dropped  from 
118  kbar  to  102  kbar,  so  that  the  value  of  )  is  about  4  kbar.  In  both 
cases,  the  value  of  Y  is  approximate,  and  it  is  not  certain  that  its 
value  for  1060  aluminum  is  really  smaller  than  it  is  for  2024  aluminum. 
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FIG.  12  PEAK  PARTICLE  VELOCITY  IN  ALUMINUM  vs.  x/xQ 
(b)  1060  Aluminum 


C.  COPPER  WEDGE-SHAPED  SPECIMEN 

Copper  specimens  were  cut  from  forgings  of  OFHC  copper  having  a 
density  of  8.94  g/cc.  Principal  impurities  in  the  copper  were  0.001  per¬ 
cent  silver  and  possibly  0.002  percent  nickel  (from  semiquant i tat i ve 
analyses).  After  machining  and  polishing,  the  samples  were  heat  treated 
at  a  temperature  of  1100°F  for  an  hour  in  an  atmosphere  of  dry  hydrogen. 

Results  of  two  experiments  in  which  aluminum  flyer  plates  hit  copper 
targets  are  shown  in  Fig.  13.  As  pointed  out  previously,  the  lesser  value 
of  the  shock  velocity  in  copper  made  it  easier  to  analyze  records  from 
copper  experiments  than  those  from  aluminum  experiments.  Results  of 
Shot  10,366  are  in  reasonable  agreement  with  the  hydrodynamic  solution 
out  to  6  plate  thicknesses.  There  is  an  indication  that  the  flyer  plate 
velocity  may  have  been  greater  than  0.125  cm/fiscc.  Shot  10,3b5  indicated 
tlat  some  attenuation  occurs  at  5.5  to  6  plate  thicknesses.  Hence  this 
pair  of  copper  shots  gave  results  which  are  more  consistent  than  are  the 
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results  from  either  pair  of  aluminum  shots  referred  to  above, 
reproducibility  in  this  case,  as  compared  to  the  aluminum  shots,  mav  be 
because  the  two  copper  shots  used  Baratol  which  was  from  the  same  order. 


FIG.  13  PEAK  PARTICLE  VELOCITY  IN  COPPER  v».  */*0 


The  peak  particle  velocity  for  copper  attenuates  from  approximut el \ 
0.04  to  approximately  0.034  cm  ,.sec  between  3  and  7  plate  thicknesses. 

By  using  the  Hugonint  for  copper,  the  stress  attenuates  from  |63  kbar  to 
1.36  kbar.  or  27  kbar.  If  this  is  due  entirelv  to  the  elastic  wave  m 
copper,  the  reverse  yield  strength  is  about  7  kbar.  The  precision  of  th 
data  mav  be  inferred  from  the  scatter  in  the  results  shown  in  K»g.  11. 
The  scatter  is  about  0.002  cm  >soc  near  »  t0  8.  or  about  *3  percent. 


D.  GOLD  WEDGE-SHAPED  SPECIMEN 


The  gold  wedge  was  made  from  99.99  percent  gold  by  Western  Gold  jr.d 
Platinum  Company,  Belmont,  California.  After  casting,  the  gold  was  re¬ 
duced  50  percent  in  thickness  by  rolling  and  then  was  machined  into  the 
required  shape.  After  the  inclined  surface  was  poli-hed,  the  material 
was  softened  by  heating  to  650°C  for  ore  hour.  Because  the  sample  was 
small  (3  inches  long),  it  was  placed  inside  a  vacuum  chamber  rather  than 
being  molded  into  the  top  cover  of  the  chamber,  as  wer  the  other  wedge- 
shaped  samples.  Because  of  the  size  of  the  '■-ample,  the  flyer  pL-jte  was 
also  reduced  to  0.040  inch  in  thickness.  This  flyer  plate  was  observed 
to  have  a  velocity  of  0.133  cm/^ser,  for  which  the  probable  error  is  less 
than  1  percent. 

Results  for  the  gold  shot  are  shown  in  Fig.  14  along  with  the  results 
of  the  hydrodynamic  calculations.  At  x / x Q  =  3.3,  the  experimental  results 
give  a  particle  velocity  of  0.028  cm/^sec  which  agrees  well  with  predicted 
value,  0.0288  cm//zsec.  Some  attenuation  takes  place  between  values  of 
x/x0  of  3.3  and  6.9,  the  particle  velocity  being  reduced  from  0.028  to 
0.024  cm/^sec.  These  values  of  particle  velocities  correspond  to  Hugoniot 
pressures  of  190  and  159  kbar,  respectively.  This  is  interpreted  as  being 
the  amplitude  of  the  elastir  relief  wave  for  gold  at  the  pressures  referred 
to  above.  As  was  done  for  aluminum  and  copper,  a  value  of  the  yield 
strength,  Y,  is  inferred.  For  gold,  Poisson’s  ratio  v  is  0.42.  so  that 
Y  is  about  4  kbar.  (See  Section  6  for  the  method  for  calculating  the  con¬ 
version  factor.  ) 
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FIG.  14  PEAK  PARTICLE  VELOCITY  IN  GOLD  v*.  *  *0 
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5.  RESULTS  OF  EXPERIMENTS  USING  GAGES 


A.  C- 7  EPOXY  SPECIMENS 

A  typical  gage  record  is  shown  in  Fig.  9  and  a  brief  description  of 
the  method  of  analysis  was  given  in  Section  3B.  Three  gages  were  included 
in  each  shot,  also  as  described  in  Section  3B.  Results  of  experiments  in 
which  “bare  gages”  were  hit  by  aluminum  plates  are  given  in  Fig.  15.  The 
term  bare  gage  is  applied  to  the  arrangements  in  which  no  material  other 
than  C-7  separates  the  flyer  plate  from  the  Manganin  wire.  Because  these 
wires  were  located  at  different  depths  in  the  C-7,  attenuation  of  peak 
pressure  could  be  observed.  In  the  experiments  the  oscilloscopes  were 
triggered  individually,  so  that  the  relation  in  time  of  one  record  to 
another  is  not  known  accurately.  Placement  of  the  records  in  Fig,  15 
(and  in  subsequent  figures)  is  therefore  arbitrary,  being  based  on  calcu¬ 
lations  which  will  be  described  later.  Some  indication  of  the  reproduc¬ 
ibility  of  results  can  be  obtained  by  comparing  the  curves  for  which 
values  of  *A0  are  0.92  and  0.95.  These  differ  about  6  percent  in  peak 
pressure  and  have  considerably  different  shapes.  This  latter  disagreement 
is  particularly  distressing  because  it  indicates  that  these  gages  cannot 
be  relied  upon  to  give  the  shape  of  the  pulse. 

The  envelope  of  the  curves  shown  in  Fig.  15  shows  that  some  attenu¬ 
ation  takes  place  between  */x0  “  1.9  and  x./xq  =  2.9.  This  is  .  fair 
agreement  with  the  results  of  the  wedge  experiment  shown  in  Fig.  16  where 
attenuation  is  observed  at  */x0  =  3.  At  */*0  =  6,  the  wedge  experiment 
gives  a  pressure  of  about  26  kbar  which  agreed  well  with  the  27.5  kbar 
given  by  the  gage  at  x/xQ  =  5.8,  Fig.  16. 

B.  ALUMINUM  SPECIMENS 

When  gages  were  used  on  metal  targets,  the  Manganin  wire  was  separated 
from  the  surface  of  the  specimen  by  about  0 .  1  -  cm- th i c k  C-7  epoxy,  as  noted 
in  Section  3B.  Results  of  gages  mounted  on  targets  of  1060  aluminum  are 
shown  in  Fig.  17(a).  Reproducibility  can  be  checked  by  comparing  the  re¬ 
sults  of  the  gages  located  at  7.82  and  7.83  (three  significant  figures  are 
used  in  order  to  identify  the  curves,  not  to  indicate  the  accuracy  of  the 
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FIG.  16  PRESSURE  IN  C-7  EPOXY  vs.  x/xQ.  RESULTS  FROM  OPTICAL 
AND  TRANSDUCER  SHOTS 
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depth  of  the  wire  in  the  epoxy).  These  two  curves  are  separated  on  the 
time  scale  more  than  warranted  so  that  the  curves  can  be  more  easily 
distinguished  and  reproducibility  is  seen  to  be  good.  Considering  the 
envelope  of  the  curves,  the  gage  at  3.9  appears  to  have  been  erratic. 

If  this  gage  is  ignored,  the  gage  records  suggest  that  attenuation  in 
1060  aluminum  starts  at  values  of  x/xQ  between  3  and  5.  This  is  some  if  hat 
earlier  than  was  indicated  by  the  wedge  experiments,  Figs.  12(a)  and  12(b). 

Results  for  annealed  2024  aluminum  are  presented  in  Fig.  17(b)  where 
attenuation  appears  to  commence  for  a  value  of  x/x0  of  about  4  or  5.  At 
least  the  pulse  top  is  no  longer  flat  for  gages  located  at  3.9  and  4.9, 
which  is  one  indication  that  an  elastic  relief  wave  has  overtaken  the 
shock  front.  The  flat  top  exhibited  by  the  gage  at  7.9  could  be  explained 
by  use  of  the  elastoplastic  model  for  relating  stress  to  strain.  This  is 
discussed  in  a  later  section  where  results  of  calculations  are  compared 
with  experimental  results. 

For  experiments  with  aluminum  targets,  the  peak  pressure  observed 
for  thin  targets  was  observed  to  be  51  kbar.  When  bare  gages  were  struck 
by  flyer  plates,  the  peak  pressure  was  45  kbar.  This  difference  in  re¬ 
sponse  has  not  been  explained,  and  the  records  presented  in  Figs.  17(a) 
and  17(b)  have  normalized  to  45  kbar. 

The  impedance  mismatch  of  aluminum  and  C-7  epoxy  is  relatively  large 
so  that  a  strong  relief  wave  was  reflected  back  into  the  aluminum  at  the 
a  1 uminum- epoxy  interface.  This  relief  wave  interacted  with  the  wave  from 
the  back  side  of  the  flyer  plate  and  put  a  portion  of  the  aluminum  in  a 
state  of  tension.  The  tension  was  apparently  great  enough  to  cause  the 
aluminum  to  fracture  or  spall.  This  is  the  explanation  of  the  relatively 
long  tails  on  the  waves,  like  that  shown  by  the  gage  at  7.9  in  Fig.  1 7 1  b ) . 
That  portion  of  the  record  beyond  5.5  /r sec  is  due  to  the  spalling  of  the 
aluminum  target.  This  phenomenon  is  also  observed  for  copper  targets  and 
will  be  discussed  when  computed  profiles  are  compared  with  experimental 
profi les . 

C.  COPPER  SPECIMENS 

Results  for  gages  on  copper  are  shown  in  Fig.  18.  An  impedance  match 
solution  showed  that  a  pressure  of  30  kbar  should  be  induced  in  the  C-7 
epoxy  when  the  copper  target  was  thin.  The  gage  on  the  thinnest  copper 
target  gave  a  pressure  of  36  kbar  for  reasons  that  are  not  understood. 
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(b)  2024  Muminom 


x/x 


Decause  the  gages  at  2.0  and  3.0  gave  the  expected  pressure,  assuming  nr. 
attenuation,  it  is  probable  that  the  gage  at  1.0  malfunctioned.  Assuming 
that  the  gage  at  •*/■*•',}  =  1  did  malfunction  obviates  the  necessity  for  nor¬ 
malizing  the  remaining  records.  The  gage  records  (Fig.  18)  show  some 
attenuation  of  the  peak  pressure  at  4.9  plate  thicknesses  which  is  a  some¬ 
what  smaller  value  than  that  obtained  from  the  wedge  experiments  (Fig.  13). 

The  impedance  mismatch  between  <■*  pper  and  C-7  is  greater  than  that  for 
aluminum  and  C-7.  Hence  spaMing  should  be  even  more  probable  for  the 
copper  targets.  That  the  copper  did  spall  was  inferred  from  the  long  tail 
on  the  pressure  profiles.  It  was  also  inferred  from  the  reloading  which 
is  shown  by  most  of  the  profiles  in  Fig.  18.  For  example,  the  profile  of 
the  gage  at  2.0  shows  approximately  1  kbar  reloading  of  the  gage  at  a  time 
of  about  2.8  /-* sec.  This  reloading  was  caused  by  the  sudden  fracture  of 
the  copper  which  released  compression  waves  on  each  side  of  the  break. 
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6.  CALCULATION  OF  ATTENUATING  SHOCK  WAVES 


The  ii’.ajor  uncertainty  in  the  calculation  of  attenuating  shock  waves 
in  solids  is  the  lack  of  information  concerning  the  equation  of  state  of 
solids.  Much  work  has  been  done  on  the  assumption  that  solids  behave  as 
fluids  above  some  limiting  pressure  or  stress.  The  work  reported  here 
was  performed  in  an  attempt  to  determine  if  departure  from  fluid  behavior 
could  be  observed.  Ii  order  to  do  this,  the  results  of  the  experiments 
reported  in  Section  5  must  be  compared  with  calculated  results. 

In  some  of  the  calculated  results,  the  material  is  assumed  to  behave 
hydrodynamic  a  1 ly ,  i.e.,  material  rigidity  is  neglected.  Pressure  is 
related  to  density  by 


where  P  is  the  pressure,  p  is  the  density,  and  p0  is  the  density  at  zero 
pressure.  Values  of  the  constants  <4  and  K  are  obtained  by  requiring  the 
equation  to  fit  the  Hugoniot  data  of  the  medium.  No  adiabats  are  derived, 
it  being  assumed  that,  for  the  materials  observed  and  for  the  pressures 
encountered,  the  adiabats  are  not  distinguishable  from  the  Hugoniot  curves. 
Thus  the  term  “hydrodynamic”  is  used  in  a  restricted  sense. 

Table  1  gives  the  values  of  <4,  K,  and  p0  for  the  materials  for  which 
calculations  have  been  performed.  The  source  of  the  data  used  in  deter¬ 
mining  the  values  of  A  and  K  is  also  given  for  each  case.  Data  for  the 
C-7  epoxy  range  from  about  0.04  to  0.08  Mbar,  For  this  material,  the 
pressure  encountered  in  the  experiments  was  0.046  Mbar  or  less. 


Table  1 

VALUES  OF  CONSTANTS  FOR  MURNAGHAN 
EQUATION  OF  STATE 


material 

^0 

(g/cn3 ) 

A 

(Mbar) 

K 

SOURCE 

C-7  Epoxy 

1.18 

wma 

■ 

Reference  4 

Aluminum 

2.785 

Bn 

Reference  1 

Copper 

8.9 

4.71 

Reference  1 

Gold 

19.24 

0.305 

5.79 

Reference  8 
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lie  along  the  curve  ae.  This  curve  is  related  to  the  hydrostatic  curve 
by  the  relation 


where  p is  the  stress  along  the  curve  ae,  and  p  is  the  hydrostatic 
pressure.  The  variable  Y  represents,  what  is  usually  called  the  yield 
stress  in  simple  tension  and  is  twice  the  maximum  resolved  shear  stress. 
The  curve  bf  is  displaced  from  the  hydrostat  by (-2/31  T. 

Comparisons  of  theory  and  experiment  based  on  this  model  with 
constant  yield  strength  show  improvement  in  agreement  over  that  obtained 
with  the  fluid  model;  however,  significant  differences  still  remain. 

Good  agreement  can  be  obtained  by  assuming  that  the  yield  strength  is 
a  function  of  the  hydrostatic  pressure 


Y  =  T0  +  M(p  -  pa ) 


(3) 


where  F0  and  M  are  constants,  p  and  pa  being  defined  in  Fig.  19.  The 
hydrostatic  curve  is  represented  by 

p  -  A/j.  +  Bp:  +  C/i3  ^ 

where  p.  =  p/p0  ~  1(  P  is  the  density  and  p0  is  the  density  at  normal 
pressure.  Combinations  of  values  of  the  constants  for  aluminum,  copper, 
and  C-7  are  given  in  Table  2.  Hugoniot  data  for  aluminum  and  copper  were 
taken  from  Rice  et  al.B  while  values  of  YQ  and  M  were  estimated  from 
attenuation  experiments.  Data  for  C-7  epoxy  were  obtained  from  Keough.4 


Table  2 

VALUES  OF  CONSTANTS  FOR  ELASTOPLASTIC 
EQUATIONS  OF  STATE 


MATERIAL 

y0 

( Mb  a  r  ) 

H 

A 

(Mbir ) 

B 

(Mb.r ) 

C 

(Mb.i) 

^0 

(g/cc ) 

Aluminum 

Aluminum 

ism 

0.055 

0.0 

m 

1.29 

1.74 

2.785 

?.78b 

Aluminum 

0.0 

0.0 

0.765 

1.66 

2.785 

Copper 

0.0007 

0.031 

1.49 

0.546 

11.85 

8.936 

C-7 

0.0 

0.Q 

0.0782 

0.196 

0.221 

1.18 

C-7 

0.0 

0.0585 

0.301 

0.0726 

1.18 

C-7 

0.0 

0.0528 

0.328 

0.0379 

1.18 
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Poisson’s  ratio  was  assumed  to  be  1/3  for  the  metals  and  0.4  for  C.-7,  the 
same  as  for  lucite.^  Both  the  bulk  and  shear  moduli  are  expressed  in  terms 
of  Poisson's  ratio,  v,  and  Young’s  modulus,  E,  and  do  not  appear  explicitly 
in  the  equations  describing  the  e lastoplastic  model.  Poisson’s  ratio  is 
assumed  to  be  independent  of  pressure,  and  E  is  represented  by 


E 


-3(1 


2v) 


V± 

dV 


(5) 


where  V  is  the  specific  volume. 

The  elastic  relief  wave  decreases  the  stress  from  p  to  p  Fig.  19. 

1  k  *  /  ’  15 

This  change  in  stress10  is  given  by 


P 


X  C 


P*f 


i  -  v 

i  -  Tv 


O’ 


*  V 


(6) 


If  Y  is  not  strongly  dependent  on  the  pressure,  and  v  -  1/3,  the  usual 
result  is  obtained  for  which 


P 


x  e 


(7) 


The  assumption  that  Poisson’s  ratio,  v,  is  a  constant  is  not  supported 
by  an  experimental  observation.  The  change  in  stress  carried  by  the 
elastic  relief  wave  is  dependent  on  the  value  v,  see  Kq.  ( 6 ) .  The  elastic 
longitudinal  sound  speed  also  depends  on  v,  as  shown  in  the  following. 

Sound  speed  c  is  defined  as 


c 


(8) 


The  elastic  sound  speed  is  obtained  by  using  the  relation  between  />  and 

p  which  connects  points  e  and  /  in  Fig.  19.  The  relation  10  is 

Pxe  "  P,  =  3(1  “  -  p)/ (1  +  v)  ,  (9) 

where  pie  is  the  stress  at  point  e,  Fig.  19,  and  pe  is  the  hydrostatic 

pressure  at  the  pL  .  Use  of  Fq.  (4)  and  the  definition  of  p  permits 
Eq.  (9)  to  be  written  so  that  p  is  a  function  of  the  density,  p.  The 
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expression  for  the  sound  speed  then  becomes 


Cj  =  [-  3K2  ( 1  -  v)(dp/dV)/(l  +  v)f2 


(10) 


where  V  is  the  specific  volume.  By  use  of  Eqs.  (2)  and  (3),  the  hydro* 
static  pressure  p  can  be  replaced  by  the  stress  p  .  There  results 


c 


1 


{3(1  -  v)(dpJdp)/[{  1  +  v)(l 


+  2/3M)]f2 


(11) 


where  dp^/dp  is  the  derivative  of  the  upper  curve,  Fig.  19,  which  represents 
the  Hugoniot  data.  In  the  computer  code  where  the  difference  scheme  is 
based  on  the  method  of  characteristics,  sound  speed  is  simply  ( dp/dp , 
where  p  vs.  p  represents  the  Hugoniot  data  and  the  derivative  is  evaluated 
by  using  Eq.  (1).  For  v  =  1/3,  the  elastic  sound  speed,  from  Eq.  (11),  is 


(12) 


so  that  the  elastic  sound  speed  is  always  greater  than  the  hydrodynamic 
sound  speed,  chydr0*  The  elastic  sound  speed  is  reduced  a  few  percent 
when  M  is  not  zero  ( M  =  0.055  for  aluminum). 

Any  variation  of  v  with  pressure  would  affect  both  the  amplitude  of 
the  elastic  wave,  Eq.  (6),  and  its  velocity,  Eq.  (11).  As  pointed  out 
above,  it  is  not  known  if  the  value  of  v  depends  on  the  stress.  If  it 
increases  with  stress,  the  speed  and  amplitude  of  the  elastic  relief  wave 
are  both  increased. 

The  computer  code  which  uses  the  elastoplastic  model  for  the  equation 
of  state  employs  the  method  of  integrating  the  flow  equations  given  by 
von  Neumann  and  Richtmyer.11  This  method  uses  an  artificial  viscosity  term 
which  smears  the  shock  front  over  a  few  cells  or  meshes.  Both  linear  and 
quadratic  terms  are  used  in  the  expression  for  artificial  viscosity, 


Q  =  [l .  72A(i  +  0.\c]Ml/V 


(13) 


where  c  is  the  speed  of  sound  (the  code  uses  the  elastic  sound  speed), 

V  is  the  specific  volume  and 

A  U  -  U)+l-U}  ,  (14) 
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is  the  difference  in  particle  velocities  in  two  adjacent  cells.  This 
relation  is  used  only  when  the  material  is  being  compressed,  i.r.,  for 
A U  <  0.  Otherwise,  the  value  of  Q  is  taken  to  be  zero. 

Most  of  the  calculations  presented  in  the  remainder  of  this  report 
were  obtained  by  use  of  the  artificial  viscosity  code  (the  Q-code). 

Each  set  of  results  is  identified  with  the  values  of  1’0  and  A/  used  in 
the  equation  of  state  for  that  set.  Where  both  1'0  and  A/  are  zero  for 
both  projectile  and  target,  the  results  should  agree  with  results  obtained 
with  the  previously  described  characteristics  code.  The  agreement  is  not 
exact  because  Eqs.  (1)  and  (4)  cannot  give  the  same  results,  and  because 
of  the  smearing  of  the  shock  front  by  the  Q-code. 
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7.  COMPARISON  OF  RESULTS  OF  CALCULATIONS 
AND  RESULTS  OF  EXPERIMENTS 


A.  C- 7  EPOXY 

Some  comparisons  of  calculated  results  with  experimental  results  have 
been  made  in  the  figures  discussed  above.  The  calculated  results  were  ob¬ 
tained  by  a  computer  code  using  the  method  of  characteristics  and  a  simple 
hydrodynamic  equation  of  state;  i.e.,  rigidity  was  not  included.  The  com¬ 
parison*  show  that  early  attenuation  occurred  in  all  the  experiments.  In 
this  section,  results  obtained  by  the  artificial  viscosity  code  are  com¬ 
pared  with  the  experimental  results.  The  artificial  viscosity  code  has 
the  advantage  that  rigidity  can  easily  be  included  in  the  equation  of 
state.  However,  rigidity  need  not  be  included  in  the  equation  of  state 
for  the  Q-code.  Results  obtained  using  hydrodynamic  representations  of 
the  equations  of  state  of  both  aluminum  and  C- 7  epoxy  are  given  in  Fig.  20(a) 
for  four  different  locations  in  the  epoxy. 

Experimental  data  from  four  gages  are  superimposed  in  Fig.  20(a)  on 
the  results  of  the  calculations.  Gage  data  showed  greater  attenuation  of 
peak  pressure  than  did  the  results  of  calculations.  This  is  viewed  as 
evidence  that  the  flow  is  not  hydrodynamic.  The  situation  was  improved 
very  little  by  adding  rigidity  to  the  equation  of  state  of  aluminum  [see 
Fig.  20(b)].  These  results  led  to  the  inclusion  of  rigidity  in  the  equa¬ 
tion  of  state  of  C-7  epoxy  also. 

Calculated  peak  pressures  agree  more  closely  with  those  from  the  gages 
when  the  yield  in  C-7  is  assumed  to  be  0.6  kbar,  Fig.  21(a).  Fair  agree¬ 
ment  of  pulse  shapes  is  also  obtained.  Even  closer  agreement  is  obtained 
for  peak  pressure  by  increasing  the  value  of  the  yield  in  C-7  to  1.0  kbar, 
Fig.  21(b).  Pulse  shapes  are  in  somewhat  better  agreement  also.  In  both 
cases  Y  is  assumed  to  be  independent  of  pressure,*  i.e.,  M  -  0. 

Results  for  C-7  epoxy  are  summarized  in  Figs.  22(a)  and  22(b)  which 
show  the  peak  pressure  as  a  function  of  distance  of  penetration  of  the 
shock  front  into  the  epoxy.  Figure  22(a)  shows  the  results  of  using  the 
hydrodynamic  equation  of  state  in  the  Q-code  for  both  aluminum  and  C-7, 
i.e.,  Y0  and  M  were  zero  for  both  materials.  These  results  differ  from 
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FIC.  70  COMPARISON  Of  CALCULATED  WAVE  PROFILES  WITH  RESULTS  OF 
CAGES  FOR  C-7  EPOXY 
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PRESSURE - kbor 


FIG.  22  PRESSURE  IN  C-7  EPOXY  v*.  »/«0  FROM  EXPERIMENTS  AND 
CALCULATIONS 

(a)  Hydrodynamic  Equation  of  Stoto  for  C-7  Epoxy 
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FIG.  22  PRESSURE  IN  C-7  EPOXY  vs.  x/xQ  FROM  EXPERIMENTS  AND 
CALCULATIONS 

(b)  Elastoplastic  Equation  of  State  for  C-7  Epoxy 


those  obtained  by  the  use  of  the  code  using  tlie  method  of  characteristics. 
Part  of  the  disagreement  (about  half  )  is  due  to  the  fact  that  the  equa¬ 
tions  of  state  are  slightly  different,  see  Eqs.  (1)  and  (4).  The  remainder 
of  the  disagreement  must  be  due  to  the  way  in  which  the  artificial  vis¬ 
cosity  code  operates.  Apparently  the  true  velocity  of  the  leading  edge  of 
a  rarefaction  is  not  preserved  by  the  code. 

Including  rigidity  in  the  equation  of  state  of  aluminum  changed  the 
ca'culated  results  very  little,  Fig.  22(a).  Thus  the  difference  between 
the  experimental  and  calculated  results  shown  in  Fig.  22(a)  must  be  due  to 
a  poor  representation  of  the  equation  of  state  of  the  epoxy.  Figure  22(b) 
summarises  the  results  obtained  by  the  inclusion  of  rigidity  in  the  equa¬ 
tion  of  state  of  the  epoxy.  The  calculated  results  using  )  =  1.0  kbar, 

M  =  0  agree  fairly  well  with  the  experimental  data. 

B.  ALUMINUM 


An  attempt  was  made  to  calculate  the  response  of  gages  on  both  a  1  usninum 
and  copper  targets.  As  discussed  above,  the  gages  on  thin  aluminum  targets 
gave  unexpectedly  high  pressures,  so  that  all  the  data  were  normalized  for 
aluminum  targets.  This  unexplained  behavior  of  the  gages  makes  the  com¬ 
parison  of  calculated  and  experimental  results  of  questionable  value. 

Such  comparisons  are  made  however  for  two  gages  in  Fig.  23(a)  which  show 
that  the  calculated  peak  pressures  for  x.Xq  =  3.9  are  in  fair  agreement. 
Agreement  is  not  so  good  when  the  aluminum  specimen  is  7.9  times  as  thick 
as  the  flyer  plate.  Of  some  interest  is  the  reloading  of  the  gage  for  the 
calculated  results,  at  times  of  4.9  usee  for  the  thin  target  and  5.9  /user 
for  the  thick  target.  There  is  no  indication  of  this  reloading  in  the  two 
gage  records.  This  appears  to  be  evidence  that  too  great  a  strength  at 
fracture  (15  kbar)  was  used  in  the  calculation. 

The  calculated  results  shown  in  Fig.  23(b)  were  obtained  by  using 
5  kbar  for  the  strength  of  the  material  at  fracture.  The  shapes  of  the 
calculated  pulses  now  agree  with  the  gage  records  much  better  than  in 
Fig.  23(a) — at  least  as  far  as  the  later  parts  of  the  pulses  are  eoiiccrnod . 
These  results  do  not  prove  conclusively  that  the  spall  strength  of  aluminum 
is  closer  to  5  kbar  than  to  15  kbar.  One  reason  for  this  is  the  necessity 
for  extrapolating  the  equation  of  state  of  aluminum  to  negative  pressure, 
which  raises  questions  which  will  not  be  discussed  tn  this  report.  A  spall 
strength  for  annealed  2024  aluminum  of  10  kbar  is  consistent  with  results 
reported  by  Butcher  et  a/.12 
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GAGES  FOR  ALUMINUM  TARGETS 
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COMPARISON  OF  CALCULATED  WAVE  PROFILES  WITH  RESULTS  OF 
CAGES  FOR  ALUMINUM  TARGETS 


When  the  relation  between  stress  and  strain  is  assumed  to  be  elasto- 
plastic,  the  initial  relief  of  pressure  is  always  elastic.  This  relief 
is  accomplished  by  a  wave  which  is  visualized  as  showing  little,  if  any, 
dispersion.  Its  effect  should  be  virtually  that  of  a  rarefaction  shock 
and  should  carry  a  pressure  drop  of  about  4K.  For  a  pressure  of  about 
100  kbar,  Eq.  (2)  gives  a  value  of  Y  of  7.5  kbar.  Thus  the  shock  in 
aluminum  should  be  overtaken  by  an  elastic  wave  and  the  interaction  should 
reduce  the  pressure  by  about  30  kbar.  Following  the  interaction,  the 
profile  of  the  wave  may  again  be  almost  flat,  as  it  "/as  just  following 
the  impact  of  the  flying  plate.  The  “Q”-coc!e  handles  this  interaction 
fairly  well  as  shown  in  Fig.  23(a).  In  the  first  calculated  profile,  the 
elastic  wave  has  not  overtaken  the  shock  front.  But  its  interaction 
reduces  the  pressure  by  10  kbar  about  0.3  f/sec  after  the  shock  front 
reaches  the  gage.  This  is  the  effect  of  the  30-kbar  wave  in  the  aluminum 
mentioned  above  which,  because  of  the  impedance  mismatch,  propagates  into 
the  epoxy  as  a  10-kbar  wave. 

The  record  of  the  gage  at  x/x0  =  3.9  is  inappropriate  when  compared 
with  the  results  of  the  calculations  referred  to  above.  From  the  slope, 
dp/dt,  of  the  back  side  of  the  wave,  it  appears  that  the  initial  relief 
wave  is  more  dispersed  in  time  than  the  following  wave.  The  effect  de¬ 
scribed  above  as  virtually  being  a  rarefaction  shock  is  not  observed. 

The  decay  is  not  as  rapid  as  the  results  given  by  the  Q-code,  which  is 
expected  to  smear  out  such  shork-like  effects.  If  elastic  relief  waves 
are  involved  in  the  flow,  their  arrivals  are  net  apparent  in  the  records 
from  the  gages. 

None  of  the  gage  records  obtained  indicates  clearly  the  arrival  of 
the  elastic  relief  wave.  This  may  be  due  to  the  inability  of  the  gage 
to  respond  to  a  sudden  release  in  pressure.  Some  ” hysteres  is *'  has  been 
reported  in  other  experiments  in  which  the  pressure  on  the  gage  was 
relieved  in  a  controlled  situation. s 

For  the  second  profile  shown  in  Fig.  23(a),  the  calculated  results 
show  that  the  elastic  relief  wave  has  attenuated  the  shock  to  a  pressure 
of  32  kbar  in  the  epoxy.  The  top  of  the  wave  is  relatively  flat,  as  is 
the  top  of  the  profile  given  by  the  gag'.  Here  the  gage  record  is  con¬ 
sistent  with  the  e  1  astopl  asl  ic  model  and  is  inconsistent  will,  tin  l.vdrn- 
dynamic  model  which  always  gives  a  triangular  profile  after  the  miti.il 
flat  top  has  attenuated. 
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C.  COPPER 


J  „  \ 


Results  of  calculations  are  shown  in  Fig.  24  for  three  thicknesses 
of  copper  targets.  Spall  strength  is  arbitrarily  taken  to  be  5  kbar  and 
the  copper  spalled  in  several  places  for  each  example.  Spalling  started 
in  each  case  about  0.2  ^tsec  after  the  shock  front  hit  the  gage.  This 
means  that  the  fractures  occurred  very  close  to  the  copper-epoxy  interface. 
Calculated  and  experimental  results  agree  reasonably  well  over  most  of 
each  pulse.  This  is  surprising  because  of  the  complications  introduced 
into  the  calculations  by  the  necessity  for  permitting  the  material  to 
spall.  Spall  strengths  of  24  and  28  kbar  are  reported. 12  Use  of  such 
values  for  the  spall  strength  would  increase  the  calculated  reloading  at 
the  gage  position. 

As  in  the  case  for  aluminum  targets,  the  experimental  results  for 
copper  do  not  show  the  characteristics  of  e las t oplas t ic  flow.  Note  the 
shape  of  the  top  of  the  calculated  pulse  shown  by  the  gage  at  3.9.  There 
is  an  initial  spike  and  then  a  flat  portion  from  3.2  /xsec  to  3.5  fisec. 

This  shape  is  expected  from  the  use  of  the  e lastoplastic  model.  The  gage 
at  3.9  does  not  show  this  shape,  or  the  gage  may  not  be  able  to  respond 
to  such  a  sudden  release  of  pressure. 

As  in  the  case  for  aluminum  targets,  the  gage  on  the  thick  copper 
target  does  not  show  the  sharp  peak  characteristic  of  the  attenuation 
predicted  bv  hydrodynamic  theory, 

D.  GOLD 

No  calculations  were  performed  for  gold  in  which  the  e lastoplast ic 
equation  of  state  was  used.  However,  a  value  of  the  yield  strength,  Y , 
was  derived  from  tne  experimental  data  in  Section  4D.  It  was  necessary 
to  use  a  handbook  value,  0.42,  for  Poisson’s  ratio  for  hardened  gold  to 
convert  from  the  amplitude  of  the  supposed  elastic  wave  to  a  value  of  Y. 

The  value  of  Y  was  found  to  be  4  kbar,  so  that  it  is  expected  that  the 
yield  strength  for  gold  is  also  dependent  on  the  stress. 
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APPENDIX  I 


MECHANISM  CAUSING  FLYER  PLATES  TO  SPALL 


It  is  possible  that  the  flyer  plates  may  spall,  i.e.,  fracture  into 
two  or  more  pieces,  because  the  material  in  the  flyer  plate  is  put  into 
a  state  of  tension  at  the  instant  that  it  is  separating  from  the  brass 
buffer  plate.  This  possibility  exists  even  though  the  wave  entering  the 
aluminum  from  the  brass  has  a  constant  state  behind  the  shock  front, 
i.e.,  the  shock  is  a  uniform  shock.  Furthermore,  no  bonding  is  necessary 
at  the  interface  for  the  production  of  a  tensile  stress  in  the  flyer  plate  . 
As  shown  in  Fig.  1-1,  .the  shock  wave  reflects  from  the  free  surface  of  the 
flyer  plate  as  a  rarefaction  fan.  In  the  figure  the  rarefaction  wave  is 
assumed  to  consist  of  a  succession  of  small  waves,  each  being  represented 
by  a  line  which  divides  the  area  in  the  (x,t)  diagram  into  fields.  When 
a  wavelet  reaches  the  interface,  each  in  turn  is  transmitted  into  the 
brass  and  reflected  into  the  aluminum.  The  reflected  wavelets  are  also 
rarefactions.  At  some  time  the  part  of  the  flyer  plate  adjacent  to  the 
brass  will  acquire  a  higher  velocity  than  the  brass,  and  the  interface 
will  open  up.  Thus,  the  field  marked  E 3,  for  example,  will  have  a  stress- 
free  boundary.  The  fields  labeled  F 3,  <»3,  and  //3  are  reached  from  field 
E 3  by  successively  crossing  elements  of  the  original  rarefaction  centered 
at  the  point  V.  Thus  the  fields  named  must  have  successively  lower  pres¬ 
sures  than  field  E 3,  which  is  stress-free.  These  fields  then  represent 
locations  in  the  plate  which  are  in  a  state  of  tension.  Using  the  method 
of  characteristics,  the  basis  for  drawing  Fig.  1*1,  Fowles  and  Curran6 
estimated  the  maximum  tension  to  be  35  khar  when  the  flyer  plate  was 
driven  at  a  velocity  of  0.19  cm'>.scc.  The  spall  threshold  lor  aluminum 
may  be  well  below  35  kbar,  so  that  the  plate  fractures  soon  after  it  starts 
to  move.  The  estimate  may  not  be  ir liable,  however,  because  the  equation 
of  state  of  aluminum  had  to  ue  extrapolated  into  the  negative  pressure 
region,  and  also  because  the  work  has  indicated  that  hydrodynamics  is  not 
reliable  as  a  method  for  calculating  flow  in  solids. 
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\  (K,»)  diagram  of  flyer  plate  driven  by  a  brass  buffer  plate 
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Figure  I  was  drawn  for  illustrative  purposes  which  required  tin- 
fields  to  be  large  enough  for  labeling.  The  angular  width  of  the  rare  - 
faction  fan  is  probably  much  smaller  than  shown,  and  the  triangular 
region  including  fields  Cl,  Cl  and  (>5  is  much  smaller  than  in  the  drawing. 
Thus  the  field  of  maximum  tension,  field  //3,  begins  very  close  to  the 
interface.  If  there  are  no  tine  effects,  the  spall  fracture,  if  any, 
should  occur  very  close  to  the  rear  cf  the  flyer  plate. 

If  the  shock  is  not  uniform,  the  flyer  plate  can  be  caused  to  spall 
by  the  interaction  of  the  rarefaction  fan  centered  at  the  point  M  in 
Fig.  1-1  with  the  rarefaction  following  the  shock  front.  If  the  shock 
wave  is  sufficiently  peaked,  the  spall  fracture  can  form  near  the  front 
of  the  plate,  whereas  the  spall  caused  by  the  mechanism  described  abo\e 
should  occur  near  the  rear  of  the  flyer  plate.  Thus  tor  a  nonuniform 
shock  wave  the  plate  can  fracture  at  almost  any  location. 
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APPENDIX  II 


THE  OPTICAL  LEVER  ARM  TECHNIQUE 


The  optical  lever  arm  technique  as  described  by  Fowles3  is  appli¬ 
cable  when  a  portion  of  a  reflecting  surface  is  abruptly  bent  through  a 
small  angle  by  the  arrival  of  a  shock  wave.  Light  from  the  point  source, 
Point  A  in  Fig.  1 1 - 1 ,  is  reflected  in  the  undisturbed  surface  and  has  an 
image  at  the  point  A ' .  The  shock,  whose  normal  velocity  is  A'0 ,  contacts 
the  free  surface  at  F  and  deflects  the  surface  through  the  angle  a.  Also 
in  Fig.  1 1 - 1 ,  the  image  of  the  light  source  A  is  reflected  in  the  deflected 


FIG.  11-1  GEOMETRY  FOR  OPTICAL  LEVER  ARRANGEMENT 


surface  and  is  located  at  A",  This  image  is  not  fully  illuminated  because 
the  camera  accepts  a  cone  of  light  from  a  point  source  and  point  F  must 
move  upward  before  the  full  cone  is  reflected  in  the  deflected  surface  to 
the  camera.  In  Fig.  1 1  - 1  the  axis  of  the  camera  is  horizontal,  and 


* 


normals  to  the  original  and  deflected  surfaces  are  inclined  at  the  angles 
6  and  B  to  the  horizontal.  The  distance,  d,  between  points  .4  and  B  is 
known,  and  the  displacement  of  the  traces  on  the  film  is  proportional  to 
the  distance  CE.  From  Fig.  1 1 - 1  it  is  evident  that 


BC  =  d  tan  6 

CD  =  BC  cos  9  -  d  sin  6 

BE  =  d  tan  (a  -  B) 

DE  =  BF  cos  6  =  d  tan  (a  -  /3)  cos  6 

CE  =  a  =  CD  +  DE 

6  -  a  +  B 


so  that 


tan  (2a) 


( a/d )  cos  6 
1  -  (a/d)  sin  6 


( II- I ) 


The  value  of  a  must  be  obtained  from  the  displacement,  J,  of  the  trace  on 
the  film,  so  that  the  magnification  of  the  optical  system  is  needed. 

Suppose  another  source  of  light  is  at  the  point  M,  Fig.  1 1 - 1 ,  at  a  dis¬ 
tance  G  from  the  source  4.  The  images  of  sources  4  and  M  as  seen  in  the 
undisturbed  surface  appear  to  be  separated  by  a  distance  G  cos  9.  On 
the  screen  of  the  Vanguard  Motion  Analyzer  the  two  sources  are  separated 
by  a  distance  Z,  so  that  the  magnification  is  G  cos  6/Z.  Thus  the  quantity 
a  is  found  to  be 


a  =  ( JG  cos  1 9)/Z 


Equation  ( 1 1 -  I )  then  becomes 


tan  (2a) 


(JG  cos2  6) /Zd 
1  -  (JG  sin  6  cos  6)/7.d 


( 1 1  -  2  ) 


Where  the  value  of  6  is  zero,  the  result  reduces  to  that  obtained 
by  Fowles.  The  velocity  of  the  free  surface  in  the  direction  of  the 
normal  to  the  original  surface  is 

tan  a  (11-3) 

where  q0  is  defined  in  Fig.  1 1  - 1 . 
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This  is  not  the  true  velocity  of  the  surface,  because  particles  will  move 
perpendicularly  to  the  undisturbed  surface  only  if  the  angle  0  is  zero. 

A  better  approximation  is  given  in  the  following  paragraph. 

As  the  shock  advances,  the  point  F  in  Fig.  I I - 1  moves  past  point  C 
and  the  image  A'  disappears.  Until  that  time,  the  image  A'  is  recorded 
on  the  film  as  a  streak  because  of  the  rotation  of  the  mirror  in  the  smear 
camera.  Many  light  sources  are  used  in  a  shot,  and  the  termini  of  the 
streaks  as  viewed  in  the  undeflected  surface  define  a  curve  on  the  film 
(Figs.  5  and  6  in  the  main  report).  Differentiation  of  this  curve  gives 
the  velocity  q0,  the  "apparent  velocity”  of  the  shock  along  the  free  sur¬ 
face.  Previous  analyses  used  the  nomenclature  of  one-dimensional  shock 
wave  theory;3  oblique  or  two-dimensional  shock  wave  relations  are  pre¬ 
ferred  here.  Referring  to  Fig.  1 1 - 2 ,  and  using  the  nomenclature  of 
Courant  and  Friedrichs,  q  Q  is  the  velocity  of  the  oblique  shock  and  0  is 
the  angle  between  the  vector  q0  and  the  shock  line.13  The  velocity  of 
flow  behind  the  shock  is  represented  by  the  vector  q^,  which  intersects 
q0  at  the  angle  8,  the  angle  of  turning  of  a  streamline  at  the  shock. 

The  velocity,  Uj  =  /VQ  -  is  the  particle  velocity  behind  the  shock  wave. 


FIG.  11-2  DIAGRAM  OF  TWO-DIMENSIONAL  SHOT 


From  the  diagram, 


Uj  cos  0 


tan  8 


q0  -  u j  sin  0 


( II-4) 


f 


from  which 


qQsin  b 
cos  { <p  -  b  ) 


(1 1-5) 


Hydrodynamic  theory  predicts  the  angle  of  turning  of  the  free  surface  to 
be  twice  the  value  of  b,  i.e.,  <X  =  2b.  Thus  the  value  of  b  is  obtained 
from  the  optical  lever  arm  equation.  The  value  of  0  can  be  determined 
for  the  two-dimensional  experiments  by  referring  to  Fig.  1-1,  where  the 
shock  is  represented  by  a  straight  line  for  convenience.  If  the  wedge 
angle  is  represented  by  y,  then 

cr  =  >  +  0 

The  detonation  advances  toward  the  apex  of  the  wedge  with  the  velocity 
D  while  point  F  moves  with  the  velocity  qQ.  Using  the  sine  law 

D/ sin  0  =  q„/sin  (180  ~  cr)  =  q0/ sin  cr  . 

Eliminating  cr  gives 

0  sin  y  ,  . . 

tan  0  =  — — ——  (  1 1  - 6  ) 

q0  ~  D  cos  y 

so  that  F,q.  ( 1 1  -  5 )  can  be  solved  for  the  particle  velocity  immediately 
in  back  of  the  shock  front.  Thus  values  of  Uj  may  be  found  for  many 
points  on  the  face  of  the  wedge. 

The  components  of  velocity,  /Vfl  and  are  related  to  density  and 
pressure  by  the  conservation  of  mass  and  momentum  equations 

Po^'o  y  p,‘v,  o1*-* 

-  5  *  r,  •  .  in-si 

In  these  equations,  p  is  the  density  and  P  is  the  pressure,  the  subscript 
0  refers  to  the  state  in  front  of  the  shock,  and  subscript  1  refers  to 
the  state  behind  the  shock. 

In  addition, 

«*,  »  <V0  -  A',  ,  (II-1)) 


J 
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so  that 


\ 

J 


Pi  •  Po  =  PoNoux 


But 


Nq  =  qQ  sin  0 

so  that 


( 1 1- 10) 


Pl  ~  P0  =  0 u  1  sin  ^  •  ( 1 1  - 1 1 ) 

Thus  the  pressure  of  the  shock  wave  can  be  calculated  for  any  point  on 
the  wedge  face.  The  density  behind  the  shock  is  obtained  from  the  relation 


P  0  U1 

1 - -  -  .  (11-12) 

P  <70  sin  0 


When  the  value  of  the  angle  0  is  nearly  90.0  degrees,  the  quantity 
q0  sin  0  becomes  equal  to  N0  and  the  results  are  the  same  as  those 
derived  by  Fowles. 3  The  use  of  the  oblique  shock  relations  obviates  the 
annoying  question  about  the  direction  of  the  motion  of  the  free  surface. 

In  the  case  ol  the  one -dimens iona 1  experiments,  it  ir  not  possible 
to  calculate  the  value  of  the  angle  0,  that  is,  the  angle  of  incidence 
of  the  shock  on  the  wedge  face  is  not  known,  principally  because  of  the 
curvature  of  the  flyer  plate.  For  this  reason,  those  experiments  do  not 
yield  equation  of  state  data.  The  wedge  angle  is  used  as  an  approximation 
in  the  analysis  so  that  an  approximate  value  of  the  particle  velocity  u{ 
is  obtained.  An  approximate  shock  velocity,  Af 0 .  may  then  be  calculated 
if  the  equation  of  state  of  the  specimen  material  is  known.  A  new  approxi¬ 
mation  to  the  value  of  0  is  then  obtained  from 


sin  0  •  A0/q0  (11-13) 

where  q0  is  known  from  the  camera  record.  Thus,  a  closer  approximation 
for  the  value  of  u(  can  be  obtained. 
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